如何使用Biopython将筛选后的FASTA序列保存至TXT文件？

有两种简单的方法可以实现需求，推荐优先使用Biopython原生的SeqIO.write方法，更规范且不易出错：

方法1：使用SeqIO.write（推荐） #

该方法会自动处理FASTA格式的规范，无需手动拼接序列头和序列：

from Bio import SeqIO

# 收集所有长度≥30的序列记录
filtered_records = []
for seq_record in SeqIO.parse("aminoacid_example.txt", "fasta"):
    if len(seq_record.seq) >= 30:
        filtered_records.append(seq_record)

# 将筛选后的序列写入目标文件，指定格式为fasta
SeqIO.write(filtered_records, "filtered_sequences.txt", "fasta")

方法2：手动写入文件（适配现有打印逻辑） #

如果你想沿用当前的打印逻辑，只需将输出定向到文件即可，用with语句可以自动管理文件的打开和关闭：

from Bio import SeqIO

# 打开输出文件，w表示覆盖写入，若需追加内容可改为a
with open("filtered_sequences.txt", "w") as out_file:
    for seq_record in SeqIO.parse("aminoacid_example.txt", "fasta"):
        if len(seq_record.seq) >= 30:
            # 若要和你的输出示例一致（>后带空格），则写f"> {seq_record.id}"
            print(f">{seq_record.id}", file=out_file)
            print(seq_record.seq, file=out_file)

两种方法最终都会生成符合要求的FASTA格式文件，方法1更符合Biopython的最佳实践，能避免手动处理格式时可能出现的错误。

内容的提问来源于stack exchange，提问作者Vykov