Snakemake 'MissingOutputException'

Question

输入： 使用 SSP 软件计算测序数据的各种质量指标的 Snakefile。 SSP 的输入是一个 BAM 文件。

sample1.sorted.bam

输出：各种文件，但我唯一关心的是一个名为{prefix}.stats.txt的文件。

sample1.stats.txt

Snakefile： ($SCIF_DATA = /scif/data)

configfile: "config.yaml"
workdir: "/scif/data"

# define samples
SAMPLES, = glob_wildcards("raw_data/{sample}.fastq.gz")

rule all:
    input:
        expand("processed_data/qc/{sample}/{sample}.stats.txt", sample=SAMPLES),

rule quality_metrics:
    input:
        "processed_data/{sample}.sorted.bam"
    params:
        prefix="{sample}",
        gt="raw_data/hg38.chrom.sizes"
    output:
        "processed_data/qc/{sample}/{sample}.stats.txt"
    shell:
        "scif run ssp '-i $SCIF_DATA/{input} -o {params.prefix} --gt {params.gt} -p 50 --odir $SCIF_DATA/{params.prefix}'"

当我在终端上运行ssp '-i sample1.sorted.bam -o sample1 --gt {params.gt} -p 50 --odir sample1 时，我得到了正确的输出：

{path}/sample1/sample1.stats.txt

但是，当我运行我的 snakemake 工作流程时，我收到以下错误：

Waiting at most 5 seconds for missing files.
MissingOutputException in line 58 of /scif/data/Snakefile:
Missing files after 5 seconds:
processed_data/qc/THP-1_PU1-cMyc_PU1_sc_S40_R1_001/THP-1_PU1-cMyc_PU1_sc_S40_R1_001.stats.txt
This might be due to filesystem latency. If that is the case, consider to increase the wait time with --latency-wait.
Will exit after finishing currently running jobs.
Shutting down, this might take some time.

增加延迟等待时间无济于事。

有什么想法吗？

Answer 1

我认为您在前缀中缺少部分输出路径

params:
    prefix="processed_data/qc/{sample}"