各位大佬们好,我在使用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.