解析fasta文件

sequences={}
ac = ''
seq = ''
for line in open("swissprot.fasta"):
    if line.startswith(">") and seq != '':
        sequences[ac]=seq
        seq = ''
    if line.startswith('>'):
        ac=line.split('|')[1]
    else:
        seq=seq+line.strip()
sequences[ac]=seq
print sequences.keys()
print sequences['P62258']
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容