Snakemake中bwa_map规则的MissingInputException错误修复求助

问题分析 #

你遇到的MissingInputException核心原因是Snakemake找不到指定的输入文件，同时你的Snakefile存在语法和命令格式问题，需要逐一修正。

分步解决方案 #

1. 验证输入文件的实际存在性 #

首先确认报错中提到的FASTQ文件是否真实存在于指定路径，在终端执行：

ls /home/muhammad/MinION_1stRun_fastq_pass/guppy_demult/fastq_all_barcodes/barcode01.fastq_runid_3c11097ede35a7f061046232b44e0228595f2f23_0.fastq

如果命令返回"文件不存在"，检查以下几点：

• 路径中的文件夹名称是否拼写正确（比如guppy_demult、fastq_all_barcodes）
• 文件名是否完全匹配（注意大小写、下划线等细节）
• 确认FASTQ文件是否真的在fastq_all_barcodes目录下

同时验证参考基因组路径的正确性：

ls /home/muhammad/miniconda3/data/G6PD_RefSeqGene_FASTA.fa

2. 修正Snakefile的语法与命令逻辑 #

你的Snakefile存在缩进错误、命令格式错误，以及未创建输出目录的问题，修正后的版本如下：

rule bwa_map:
    input:
        genome = "miniconda3/data/G6PD_RefSeqGene_FASTA.fa",
        samples = "MinION_1stRun_fastq_pass/guppy_demult/fastq_all_barcodes/barcode01.fastq_runid_3c11097ede35a7f061046232b44e0228595f2f23_0.fastq"
    output:
        sample_bams = "mapped_reads_allsamples/barcode01.fastq_runid_3c11097ede35a7f061046232b44e0228595f2f23_0.bam"
    shell:
        """
        mkdir -p mapped_reads_allsamples
        bwa mem {input.genome} {input.samples} | samtools view -Sb - > {output.sample_bams}
        """

关键修改点：

• 对齐input区块内samples的缩进，保证Python语法合规
• 在shell命令中明确调用{input.genome}和{input.samples}，符合bwa mem的标准语法（参考基因组在前，测序reads在后）
• 添加mkdir -p mapped_reads_allsamples，自动创建输出目录（Snakemake不会默认创建目录）
• 使用三重引号包裹多行shell命令，提升可读性

3. 测试与运行 #

修正后先执行干运行测试，确认无路径或语法错误：

snakemake -n

如果输出无报错，再执行实际运行：

snakemake

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