from Bio import SeqIO import os folder_path = 'C:/Users/fang/Desktop/合适的程序/序列/M/MC0' # 文件夹路径 output_file = 'MC0/MC0.fasta' # 合并后的fasta文件名 sequences = [] # 初始化序列列表

遍历文件夹下的所有fasta文件

for file_name in os.listdir(folder_path): if file_name.endswith('.fasta'): file_path = os.path.join(folder_path, file_name) # 构造文件路径 # 读取fasta文件 with open(file_path, 'r') as f: for record in SeqIO.parse(f, 'fasta'): record.id = record.description # 将description作为id record.description = '' # 将description清空 sequences.append(record) # 将读取的序列添加到序列列表

将所有序列写入输出文件中

with open(output_file, 'w') as f: SeqIO.write(sequences, f, 'fasta')

合并FASTA序列文件并保持原始格式

原文地址: https://www.cveoy.top/t/topic/lLai 著作权归作者所有。请勿转载和采集!

免费AI点我,无需注册和登录