Snakemake合并Shell命令行报错：求合法合并方式以避免中间文件

你的语法错误根源在于把shell命令拆成了两个独立的引号字符串，Snakemake会把它们当作两个分开的指令处理，导致管道符|无法正确连接前后的命令流水线。下面是正确的合并方式和相关注意事项：

问题原因解析 #

你当前的shell字段将bedtools intersect和bcftools annotate命令拆成了两个用引号包裹的片段，中间的管道符没有被包含在一个完整的命令序列里。这种写法会让Snakemake解析时误认为是两个独立的shell命令，从而抛出语法错误。

正确的合并方法 #

1. 写成单个完整字符串 #

把所有命令（包括管道符）放在一个连续的字符串中，不需要拆分引号：

rule bedtools:
    input:
        invcf="/path/to/my.vcf.gz",
        bedgz="/path/to/my.bed.gz"
    output:
        outvcf="my.filtered.vcf.gz"
    shell: "/Tools/bedtools2/bin/bedtools intersect -a {input.invcf} -b {input.bedgz} -header -wa | /Tools/bcftools/bcftools annotate -c CHROM,FROM,TO,GENE -h <(echo '##INFO=<ID=GENE,Number=1,Type=String,Description=\"Gene name\">') > {output.outvcf}"

2. 使用多行字符串提升可读性 #

如果命令过长，用三引号包裹多行字符串，这样结构更清晰，便于维护：

rule bedtools:
    input:
        invcf="/path/to/my.vcf.gz",
        bedgz="/path/to/my.bed.gz"
    output:
        outvcf="my.filtered.vcf.gz"
    shell: """
        /Tools/bedtools2/bin/bedtools intersect -a {input.invcf} -b {input.bedgz} -header -wa |
        /Tools/bcftools/bcftools annotate -c CHROM,FROM,TO,GENE -h <(echo '##INFO=<ID=GENE,Number=1,Type=String,Description=\"Gene name\">') > {output.outvcf}
    """

关键注意点：进程替换的兼容性 #

你用到了bash特有的进程替换语法<(...)，但Snakemake默认使用的是sh（不支持该语法），所以必须确保用bash执行命令。有两种解决方式：

方式A：在rule内指定bash执行

shell: """
    bash -c '
        /Tools/bedtools2/bin/bedtools intersect -a {input.invcf} -b {input.bedgz} -header -wa |
        /Tools/bcftools/bcftools annotate -c CHROM,FROM,TO,GENE -h <(echo \"##INFO=<ID=GENE,Number=1,Type=String,Description=\\\"Gene name\\\">\") > {output.outvcf}
    '
"""

注意这里的引号转义：外层用双引号包裹，内层bash命令用单引号，而命令里的双引号需要转义。

方式B：全局设置默认shell为bash

在Snakefile最开头添加一行，让所有规则的shell命令都用bash执行：

shell: "/bin/bash"

workdir: "/path/to/workdir/"
rule all:
    input: "my.filtered.vcf.gz"
# ... 后续规则 ...

这种方式更简洁，适合多个规则都需要bash特性的场景。

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