mashengwei的个人博客分享 http://blog.sciencenet.cn/u/mashengwei

博文

fasta转成txt

已有 6544 次阅读 2017-2-10 23:10 |系统分类:科研笔记| fasta

fasta转成txt

有的时候需要将fasta格式的序列文件转换成以tab键分割的txt文件

写了一个简单的命令。用法如下:

fasta2txt -i input.fa -o out.txt

代码如下

#!/usr/bin/env python
# -*- coding: utf-8 -*-
__author__ = "Sheng-Wei Ma"
import click   #需要你的电脑安装click包
@click.command()
@click.option('-i', type=click.File('rb'), help='fasta file')
@click.option('-o', type=click.File('w'), help='txt file')
def fasta2txt(i, o):
   for line in i:
       if line.startswith(">"):
           lin = line.strip().split()[0][1:]
           o.writelines('n'+lin+'t')
       else:
           o.writelines(line.strip())
if__name__=='__main__':
   fasta2txt()

别忘了上述代码保存为fasta2txt,并修改权限使之可运行,最后可将文件放在系统全局环境变量下。




https://blog.sciencenet.cn/blog-1094241-1032900.html

上一篇:MAC清除mackeeper
下一篇:“高通量、单细胞是核心,后面的都是扯淡”总结
收藏 IP: 58.213.93.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-12-23 19:28

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部