Snakemake中如何通过配置文件设置最终输出文件名

问题根源 #

你遇到的"Target rules may not contain wildcards"错误，核心原因是all作为默认目标规则，其输出必须是明确的固定路径/文件名，不能包含Snakemake的通配符（比如{sample}这类动态占位符）。你之前的写法可能误将配置变量按通配符逻辑处理，或是规则定义的语法存在问题。

解决方案步骤 #

1. 配置config.yml

在配置文件中明确指定最终输出文件名，同时可维护输入VCF列表等其他参数：

# config.yml
output_filename: "merged_multisample.vcf"
input_vcfs:
  - sample1_raw.vcf
  - sample2_raw.vcf

2. 修改Snakefile

加载配置文件后，直接在规则中引用配置变量作为输出路径，确保all规则的输出为固定值：

# Snakefile
# 加载配置文件
configfile: "config.yml"

# 单个VCF提取变异规则（根据你的实际提取逻辑调整）
rule extract_variants:
    input:
        lambda wildcards: wildcards.raw_vcf
    output:
        "extracted/{sample_id}.vcf"
    shell:
        """
        # 替换为你的变异提取命令，示例用bcftools筛选PASS变异
        bcftools view -f 'PASS' {input} > {output}
        """

# 多样本VCF合并规则
rule merge_vcfs:
    input:
        # 自动生成所有提取后的VCF文件列表
        vcfs = expand(
            "extracted/{sample_id}.vcf",
            sample_id=[vcf.split('_raw')[0] for vcf in config["input_vcfs"]]
        )
    output:
        # 直接引用配置中的输出文件名
        config["output_filename"]
    shell:
        """
        # 替换为你的合并命令，示例用bcftools merge
        bcftools merge -o {output} {input.vcfs}
        """

# 默认目标规则：指向最终输出文件
rule all:
    input:
        config["output_filename"]

3. 运行命令

直接执行指定核心数的命令即可：

snakemake --cores 4

关键说明 #

• all规则的输入直接绑定配置文件中的固定文件名，彻底规避了通配符，解决了报错问题。
• 输出文件名的控制权完全集中在config.yml，修改配置即可切换结果文件名，无需改动Snakefile逻辑。
• 通过expand()函数结合配置中的输入VCF列表，可自动生成样本ID和对应输入文件路径，避免硬编码样本信息。

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