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']
解析fasta文件
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成，浏览时请结合常识与多方信息审慎甄别。
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。
【社区内容提示】社区部分内容疑似由AI辅助生成，浏览时请结合常识与多方信息审慎甄别。
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。
推荐阅读更多精彩内容
- 之前了解了一点解析ips文件的方式，觉得每次都要敲命令行老不爽了，于是乎想要一个一键解析的方式。学习了一会Shel...
- Dom4j是一个简单灵活的开源代码的库要使用它，就要下载它所对应的jar文件1，使用Dom4j通过dom解析解析x...
