假如以下这个脚本为Gene_length.py, 测试的fasta文件为test.fasta,输出文件为outfile
代码的运行命令(在linux中):python Gene_length.py test.fasta outfile
* 其中 ID = re.split(r'\.',line)[2] 这一步是因为原始fasta文件里ID名称是这样的>evm.model.Contig62652.1.我们需要的是Contig62652这一部分,这个ID是用.来分割,而我们需要的是第三部分,在列表中是属于【2】
Python | 统计FAST文件里多个基因的长度
最后编辑于 :
?著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事?!?“怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...