使用Biopython按坐标删除序列

from Bio import SeqIO
from Bio.SeqRecord import SeqRecord
from Bio.Seq import Seq

input_file = "input.fasta"  # 输入文件名
output_file = "output.fasta"  # 输出文件名
seq_records = SeqIO.parse(input_file, "fasta")  # 解析fasta格式文件

new_records = []
for record in seq_records:
    seq = str(record.seq)  # 获取序列的字符串表示
    start = 10  # 删除序列的起始坐标（示例值）
    end = 30  # 删除序列的结束坐标（示例值）
    new_seq = seq[:start-1] + seq[end:]  # 删除指定坐标范围内的序列
    new_record = SeqRecord(Seq(new_seq), id=record.id, description=record.description)  # 创建新的序列记录
    new_records.append(new_record)

SeqIO.write(new_records, output_file, "fasta")