搜索
楼主: Jimmy

详细讲解如何抽取fasta文件里面指定序列名的序列

[复制链接]

17

主题

42

帖子

270

积分

中级会员

Rank: 3Rank: 3

积分
270
发表于 2017-8-24 02:00:46 | 显示全部楼层
银色麦穗 发表于 2016-9-12 20:15
我用的是python,也写过类似的一个,原理是一样的。从命令行里读取要提取的序列,如Chr1,然后就会输出到屏 ...

这个flag这里起什么作用,感觉没看懂,是对每一行的标记吗?为什么flag=0就不输出行,flag不等于0就输出行。楼主能详细解释一下吗。万分感激
回复 支持 反对

使用道具 举报

0

主题

15

帖子

215

积分

中级会员

Rank: 3Rank: 3

积分
215
QQ
发表于 2019-1-15 20:37:23 | 显示全部楼层
learnyoung 发表于 2017-1-17 15:05
对于fasta和gennebank文件的解析可以用biopython这个包,我选取ls_orchid 前十个gene的id存入gene_id 文件 ...

这个运行的速度确实不快,比较慢,但是能用。就是输出结果怎么写到如新文件,还在尝试中
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|手机版|小黑屋|生信技能树 ( 粤ICP备15016384号  

GMT+8, 2019-11-20 21:07 , Processed in 0.027261 second(s), 23 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.