Snakemake does not give output files. Gives no errors
0
0
Entering edit mode
2.8 years ago
kamanovae ▴ 100

Hi!

My cod looks like this:

import os

(SAMPLES,) = glob_wildcards("../exom/NIST7035/{sample}_L001_R1_001.fastq")
SAMPLE_PATH=os.path.dirname(config["sample"])

rule all:
        input:
                expand(config["out"] + "/{sample}/{sample}_R1.fastq", sample=SAMPLES),
                expand(config["out"] + "/{sample}/{sample}_R2.fastq", sample=SAMPLES)

rule sort:
        output:
                output1=config["out"] + "/{sample}/{sample}_R1.fastq", output2=config["out"] + "/{sample}/{sample}_R2.fastq"
        shell:
                " zcat  {SAMPLE_PATH}/{wildcards.sample}*R1*.fastq.gz | paste - - - - |  sort -k1,1 -S 30G | tr '\t' '\n' > {output.output1} ;"
                " zcat  {SAMPLE_PATH}/{wildcards.sample}*R2*.fastq.gz | paste - - - - |  sort -k1,1 -S 30G | tr '\t' '\n' > {output.output2}"

The snakemake terminates without output files without throwing an error

Building DAG of jobs...
Using shell: /bin/bash
Provided cores: 2
Rules claiming more threads will be scaled down.
Job stats:
job      count    min threads    max threads
-----  -------  -------------  -------------
all          1              1              1
total        1              1              1

Select jobs to execute...

[Thu Feb  3 12:43:46 2022]
localrule all:
    jobid: 0
    resources: tmpdir=/tmp

[Thu Feb  3 12:43:46 2022]
Finished job 0.
1 of 1 steps (100%) done
Complete log: /storage1/GatkBwaTest/SnakemakeDir/.snakemake/log/2022-02-03T124346.250486.snakemake.log

Maybe someone knows what I'm doing wrong? Thanks!

Snakemake • 555 views
ADD COMMENT
0
Entering edit mode

check if SAMPLES are listing files, check if ouptut directories and expected files already exist in output directories. There are few other issues with the code. As far as I know, this part of the code {wildcards.sample}*R1*.fastq.gz should fail with warning.

ADD REPLY

Login before adding your answer.

Traffic: 2647 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6