Snakemake技术咨询：如何在rule all中引入配置参数

问题根源 #

你遇到的两类错误本质是两个核心问题：

• expression cannot contain assignment：rule all的input列表中第二个expand未缩进，且试图引用规则内部的params.xxx变量——expand属于全局作用域，无法识别规则内的局部变量，同时语法格式也不符合要求。
• Missing files for rule all：生成的目标文件名与AdapterRemoval_se规则实际输出的文件名不匹配，比如你在expand里写错了参数名（params.qual应为params.m_qual），且规则输出路径的占位符与params变量名不一致（{params.mlength} vs params里的m_length），导致Snakemake无法关联目标文件与对应规则。

解决方案 #

1. 修正rule all的目标文件生成逻辑 #

直接从config中提取参数，用全局变量在expand里生成目标路径，避免引用规则内的params：

configfile:"config.yaml" 

print(config)

SAMPLES= ["HgSim"]

rule all:
    input:
        expand("01_fastQC/{samples}_fastqc.log", samples=SAMPLES),
        # 直接用config变量生成路径，确保和规则输出完全匹配
        expand(
            "Resume_Stats/AR.mlength{mlength}.mqual{mqual}/{samples}.settings",
            samples=SAMPLES,
            mlength=config["mlength"],
            mqual=config["mqual"]
        )

2. 统一规则输出的占位符与参数名 #

确保AdapterRemoval_se规则的输出路径占位符和params变量名完全一致，消除文件名匹配误差：

rule AdapterRemoval_se:
    input: "reads/{samples}.fq"

    output: 
        # 统一使用params.m_length和params.m_qual，与params定义对齐
        Settings="Resume_Stats/AR.mlength{params.m_length}.mqual{params.m_qual}/{samples}.settings",
        Truncated="02_AdapterRemoval_pe/mlength{params.m_length}.mqual{params.m_qual}/{samples}.truncated",
        Discarded="02_AdapterRemoval_pe/mlength{params.m_length}.mqual{params.m_qual}/{samples}.discarded"
    
    threads: 20   

    params:
        adapter_1=config["adapt1"],
        m_length=config["mlength"],
        m_qual=config["mqual"]
    
    shell: "AdapterRemoval --file1 {input}  --discarded {output.Discarded}  --settings {output.Settings} --output1 {output.Truncated} --adapter1 {params.adapter_1} --threads {threads} --minlength {params.m_length} --minquality {params.m_qual} --trimns --trimqualities "

关键说明 #

• expand在Snakemake解析规则前就会执行，只能使用全局可用变量（config、全局列表/字典等），无法直接引用规则内部的params、wildcards。
• 目标文件路径必须和对应规则的输出路径完全一致，包括占位符拼写，否则Snakemake无法建立依赖关系，会判定文件缺失。

进阶文档参考 #

• 官方文档的Advanced Workflows章节：涵盖参数化工作流、动态规则生成、配置文件高级用法，比基础教程更深入。
• 官方文档的Configuration进阶部分：讲解如何用config实现复杂参数分支、多环境配置。
• Snakemake社区整理的Best Practices：包括规则模块化、参数管理、依赖优化的实用技巧，适合进阶学习。

内容的提问来源于stack exchange，提问作者angelina.Z.H