博文

fasta转成txt

已有 6034 次阅读 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://wap.sciencenet.cn/blog-1094241-1032900.html

上一篇：MAC清除mackeeper
下一篇：“高通量、单细胞是核心，后面的都是扯淡”总结

收藏 IP: 58.213.93.*| 热度|

当前推荐数：0

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

数据加载中...

返回顶部

马省伟

扫一扫，分享此博文

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

博文

fasta转成txt

当前推荐数：0

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

马省伟

全部作者的精选博文

全部作者的其他最新博文

全部精选博文导读

相关博文

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

博文

fasta转成txt

当前推荐数：0

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

马省伟

全部作者的精选博文

全部作者的其他最新博文

全部精选博文导读

相关博文

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