合并基因组文件

将基因组fasta文件合并成一行：

import os
import time

start = time.time()

os.chdir('C:/Users/Administrator/Desktop/')
print(os.getcwd())

res_dict = {}

with open('YLG.sort.fasta','r') as pep:
    for line in pep:
        if line.startswith('>'):
            name = line.strip().split()[0]
            res_dict[name] = ''
        else:
            res_dict[name] += line.replace('\n','')

gen = ''
for cds_id, sequence in res_dict.items():
    #print(sequence)
    gen = gen + sequence

with open('genome.txt', 'a') as file:
        file.write(gen)

end = time.time()
print(end - start)

💌lixiang117423@foxmail.com
💌lixiang117423@gmail.com