使用Python计算fasta文件的序列长度

最新推荐文章于 2022-09-13 22:52:49 发布

张益凯的生信笔记

最新推荐文章于 2022-09-13 22:52:49 发布

阅读量4.2k

点赞数 1

CC 4.0 BY-SA版权

分类专栏： python 文章标签： python fasta length

本文链接：https://blog.youkuaiyun.com/zhangyikai2017/article/details/84654391

python 专栏收录该内容

2 篇文章

订阅专栏

本文介绍了一种使用Python脚本来计算FASTA文件中各序列长度的方法。通过读取FASTA文件，该脚本能够准确地统计每个序列的长度，并将结果输出到指定文件中。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

使用Python计算fasta文件的序列长度

在这里插入代码片使用Python计算fasta文件的序列长度

#！/usr/bin/python
#-- coding:utf-8 --

import sys

f = open(sys.argv[1],‘r’)
out = open(sys.argv[2],‘w’)

def chr_length(infile,outfile):
f = open(sys.argv[1],‘r’)
out = open(sys.argv[2],‘w’)
dic = {}
name = ‘’
seq = []
for line in f:
if line.startswith(’>’):
name = line.strip(’\n’).strip(’>’)
seq = []
dic[name] = seq
else :
line = line.strip(’\n’)
seq.append(line)
dic[name] = seq
#print “%s sequences in total” % len(dic)
out.write(“Chromosome\tstart\tend\n”)
for (name ,seq) in dic.items():
out.write("%s\t1\t%s\n" % (name,sum(map(len,seq))))
f.close()
out.close()
chr_length(sys.argv[1],sys.argv[2])