QQ登录

只需一步,快速开始

登录 | 立即注册 | 找回密码

主题

帖子

荣誉

新鱼友

Rank: 1

积分
9
查看: 275|回复: 1

哪位大神知道如何将两个文件中根据相同序列ID将序列合并

[复制链接]
最佳答案
0 
累计签到:4 天
连续签到:1 天
xiaoxiaonao 发表于 2018-4-24 17:13:14 2751 | 显示全部楼层 |阅读模式

马上注册加入鱼C,享用更多服务吧^_^

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
我有A、B两个文件,序列ID按照顺序都一致,现在想将两个文件按照相同的ID将序列以换行符合并,求哪位大神指导
A文件
>1_CHROMOSOME3780765_3785720
tgttgaaagttaaacttgattttgaatcaagtttaattattggatcaattatccaataat
taattatggccaaatccaagttctagagttttctctagaaatatcatcatttccacctcc
>1_CHROMOSOME5691759_5696703
taaatacctaaccggtttgttccgggttaattaggttagatagggaaagtgtatataatg
tcttgtattcactctttaataagataacagaaacattttccactttcggctccttacatt
B文件
>1_CHROMOSOME3780765_3785720
tgttgaaagttaaacttgattttgaatcaagtttaattattggatcaattatccaataat
taatta
>1_CHROMOSOME5691759_5696703
ttgttaatacctaaccggtttgttccgggttaattaggttagatagagaaaatgtatata
atgtcttgtattca
我想得到C文件如下
>1_CHROMOSOME3780765_3785720
tgttgaaagttaaacttgattttgaatcaagtttaattattggatcaattatccaataat
taattatggccaaatccaagttctagagttttctctagaaatatcatcatttccacctcc
tgttgaaagttaaacttgattttgaatcaagtttaattattggatcaattatccaataat
taatta
>1_CHROMOSOME5691759_5696703
taaatacctaaccggtttgttccgggttaattaggttagatagggaaagtgtatataatg
tcttgtattcactctttaataagataacagaaacattttccactttcggctccttacatt
ttgttaatacctaaccggtttgttccgggttaattaggttagatagagaaaatgtatata
atgtcttgtattca
就是以换行符将相同序列ID的序列合并了,求大神指导
楼层
跳转到指定楼层
最佳答案
831 
累计签到:346 天
连续签到:86 天
ba21 发表于 2018-4-24 18:48:20 | 显示全部楼层
  1. #把文件转换成字典
  2. def get_dict(file):
  3.     t_dict = {}
  4.     with open(file, "r") as f:
  5.         txt = f.read()
  6.         txt_list = txt.split('>')
  7.         for item in txt_list:
  8.             if item != '':
  9.                 item_list=item.splitlines()
  10.                 temp_txt = ''
  11.                 for i in range(len(item_list)):
  12.                     if i == 0:
  13.                         id_txt= '<' + item_list[i] + '\n'
  14.                     else:
  15.                         temp_txt += item_list[i] + '\n'
  16.                 t_dict[id_txt] = temp_txt
  17.     return t_dict
  18.             
  19. #获取2两文件的内容,转换成字典保存
  20. a_dict = get_dict("a.txt")
  21. b_dict = get_dict("b.txt")

  22. #把b_dict中的值合并到a_dict中
  23. for k, v in b_dict.items():
  24.     if k in a_dict.keys():
  25.         a_dict[k] += v

  26. #把合并后的内容写入文件         
  27. with open("c.txt", "w") as f:
  28.     for k, v in a_dict.items():
  29.         f.writelines(k)
  30.         f.writelines(v)
复制代码

发表回复

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

关闭

小甲鱼强烈推荐 上一条 /1 下一条

    移动客户端下载(未启用)
    微信公众号

小黑屋|手机版|Archiver|鱼C工作室 ( 粤ICP备11014136号

Copyright 2018 鱼C论坛 版权所有 All Rights Reserved.

Powered by Discuz! X3.1 Copyright
© 2001-2018 Comsenz Inc.    All Rights Reserved.

小黑屋|手机版|Archiver|鱼C工作室 ( 粤公网安备 44051102000370号 | 粤ICP备11014136号

GMT+8, 2018-5-25 07:35

快速回复 返回顶部 返回列表