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

xiaoxiaonao

我有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的序列合并了，求大神指导

ba21

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. 
    
23. #把b_dict中的值合并到a_dict中
    
24. for k, v in b_dict.items():
    
25.     if k in a_dict.keys():
    
26.         a_dict[k] += v
    
27. 
    
28. #把合并后的内容写入文件         
    
29. with open("c.txt", "w") as f:
    
30.     for k, v in a_dict.items():
    
31.         f.writelines(k)
    
32.         f.writelines(v)
    

复制代码