使用Snakemake运行Picard/GATK MergeVcfs时出错求助

一、缩进错误（IndentationError） #

Snakemake强制使用4个空格进行缩进，严禁混合制表符（Tab）和空格，否则必然触发缩进错误。排查和修正要点：

• 检查规则内所有层级块（input:/output:/shell:/wrapper:等）的缩进，确保每个子项都比父级多4个空格。
• 编辑器中设置将Tab自动转为4个空格，避免手动混用Tab和空格。
• 示例正确缩进的Picard规则：

rule merge_vcfs:
    input:
        vcf1="data/sample1.vcf",
        vcf2="data/sample2.vcf"
    output:
        "results/merged.vcf"
    shell:
        """
        picard MergeVcfs I={input.vcf1} I={input.vcf2} O={output}
        """

二、subprocess.CalledProcessError（退出状态127） #

退出码127的核心原因是系统无法找到Picard执行命令，对应解决方法：

1. 验证Picard可执行性
   在终端直接执行picard（或java -jar /path/to/picard.jar，若为jar包形式），确认能正常调用。若报错，说明Picard未安装或未添加到系统PATH中。
2. conda环境适配
   若使用conda管理依赖，需确保流程激活了包含Picard的环境，或在规则中显式指定conda环境：

   rule merge_vcfs:
       input:
           vcfs=["data/sample1.vcf", "data/sample2.vcf"]
       output:
           "results/merged.vcf"
       conda:
           "envs/picard_env.yaml"  # 环境文件需包含picard依赖
       wrapper:
           "v3.9.0/bio/picard/mergevcfs"
   

3. wrapper参数匹配
   确认wrapper的输入参数格式正确：多数Picard mergevcfs wrapper要求输入为VCF文件列表，而非单独的输入项。比如将input: vcf1=..., vcf2=...改为input: vcfs=[...]，确保与wrapper的参数定义匹配。
4. 指定Picard完整路径
   若Picard未在PATH中，可在命令中直接指定完整路径：

   rule merge_vcfs:
       input:
           vcf1="data/sample1.vcf",
           vcf2="data/sample2.vcf"
       output:
           "results/merged.vcf"
       shell:
           """
           java -jar /opt/picard/picard.jar MergeVcfs I={input.vcf1} I={input.vcf2} O={output}
           """
   

额外排查技巧 #

• 查看Snakemake输出的错误日志，找到实际执行的命令，复制到终端手动运行，确认是否能复现问题，以此快速定位是环境问题还是命令参数问题。
• 若使用wrapper，可查阅对应版本的wrapper说明，确认输入输出参数的要求。

内容的提问来源于stack exchange，提问作者Wei-Chen Pan