我有一个contigs文件（contigs名字和序列）和一个paired reads文件，请书写一段程序完成以下要求：根据contigs文件中序列信息，将paired reads文件中的reads回贴（mapping），并将成功回贴的reads导出到一个新的文件

您可以使用BioPython库来处理contigs文件和paired reads文件，并使用Bowtie2或BWA等工具来进行reads的回贴。以下是一个示例程序，演示如何完成您的要求：

from Bio import SeqIO
import subprocess

# 读取contigs文件
contigs_file = "contigs.fasta"
contigs = SeqIO.to_dict(SeqIO.parse(contigs_file, "fasta"))

# 读取paired reads文件
reads_file = "reads.fastq"
reads_output = "mapped_reads.fastq"

# 创建一个空的输出文件
with open(reads_output, "w") as output:
    pass

# 逐个回贴reads并将成功回贴的reads写入输出文件
for contig_id, contig_seq in contigs.items():
    # 创建一个临时SAM文件，用于存储回贴结果
    sam_file = f"{contig_id}.sam"

    # 使用Bowtie2进行reads回贴
    bowtie2_command = f"bowtie2 -x {contig_id} -U {reads_file} -S {sam_file}"
    subprocess.run(bowtie2_command, shell=True, check=True)

    # 解析SAM文件，将成功回贴的reads写入输出文件
    with open(sam_file, "r") as sam:
        with open(reads_output, "a") as output:
            for line in sam:
                if line.startswith("@"):
                    continue
                fields = line.split("\t")
                read_name = fields[0]
                output.write(f"@{read_name}\n")
                output.write(f"{fields[9]}\n")

    # 删除临时SAM文件
    subprocess.run(f"rm {sam_file}", shell=True, check=True)

print("回贴完成！")

请将contigs文件命名为`contigs.fasta`，paired reads文件命名为`reads.fastq`。程序将生成一个名为`mapped_reads.fastq`的输出文件，其中包含成功回贴的reads。

请注意，此示例程序使用Bowtie2进行回贴，您需要先安装Bowtie2并创建一个包含contigs序列索引的Bowtie2数据库。您还可以根据需要更改回贴工具和相关参数。

此示例程序仅供参考，您可能需要根据实际情况进行调整和优化。