各位大佬们好,我在使用snakemake搭流程中,在空运行时出现 reason: Missing output files: /public/home/guoting/lsn/chipship/snakemake_test/02Clean_data/Control_1.trim.fastq.gz这种错误,请问造成这个错误的原因是什么呢?
这是测试代码:
rule trim:
    input:
        s1="/public/home/guoting/lsn/chipship/snakemake_test/01Rowdata/Control_1.fastq.gz"
    output:
        trim1="/public/home/guoting/lsn/chipship/snakemake_test/02Clean_data/Control_1.trim.fastq.gz"
    threads:15
    params:
        min_len = 30,
        q=25,
        compression=4
    log:
        fastp_error="/public/home/guoting/lsn/chipship/snakemake_test/log/test1.fastp.log"
    shell:
        "source /public/software/anaconda3/env.sh" 
        "&"
         "source activate  /public/home/guoting/.conda/envs/rnaseq"
          "&"
         "fastp --thread={threads} -I {input.s1} -O {output.trim1} --compression {params.compression} -q {params.q} -l {params.min_len} 2>{log.fastp_error}"
这是空运行后的错误:
Building DAG of jobs...
Job stats:
job      count
-----  -------
trim         1
total        1
[Sat Sep  9 19:04:07 2023]
rule trim:
    input: /public/home/guoting/lsn/chipship/snakemake_test/01Rowdata/Control_1.fastq.gz
    output: /public/home/guoting/lsn/chipship/snakemake_test/02Clean_data/Control_1.trim.fastq.gz
    log: /public/home/guoting/lsn/chipship/snakemake_test/log/test1.fastp.log
    jobid: 0
    reason: Missing output files: /public/home/guoting/lsn/chipship/snakemake_test/02Clean_data/Control_1.trim.fastq.gz
    resources: tmpdir=/tmp
Job stats:
job      count
-----  -------
trim         1
total        1
Reasons:
    (check individual jobs above for details)
    missing output files:
        trim
This was a dry-run (flag -n). The order of jobs does not reflect the order of execution.