Python根据ID从fasta文件提取序列

weixin_38101791 weixin_38101791 墨衍会员 · AI 创作全网分发 墨衍智能分发 通过墨衍一键同步至各平台 可同步平台 CSDN微信公众号微博知乎掘金百家号博客园抖音小红书 你的文章也可以这样分发 墨衍会员 ¥399起/年,写一次发全网 我也要 + 加入社区
根据ID从fasta文件提取序列,为什么结果只有最后一个ID的结果,之前的被覆盖了?以前用这个提取过序列,没有出现问题,这次不知道为什么,只有一个ID的结果 #!/usr/bin/python3 import sys args=sys.argv fr=open(args[1], 'r') fw=open('./out.fasta', 'w') dict={} for line in fr:     if line.startswith('>'):         name1=line.split(">")[1]         name=name1.split(":")[0]         dict[name]=''     else:         dict[name]+=line.replace('\n','') fr.close() ls=open(args[2], 'r') for lin in ls:      for name in dict.keys():          if name ==lin:             fw.write(">"+ name)             fw.write('\n')             fw.write(dict[name])             fw.write('\n') fw.close()
...全文

相关推荐

评论
收藏 已收藏
分享
返回顶部