Question: missingoutputexception, snakemake after completing 95% analysis
0
gravatar for arshil
6 months ago by
arshil0
arshil0 wrote:

Hi, can anyone look at the code and tell whats wrong with the code. it runs and does fastqc analysis till 95 % and than show missingf output files and stops. can anyone help me out with this .

from util.varsub import varsub
shell.prefix("set -euo pipefail;")
configfile: "config.yaml"
varsub(config)

shell.prefix("source ~/.bash_profile")
freeze = config['freeze']

SAMPLES, = glob_wildcards(config['datadirs']['fastq'] + "/" + "{file}_1.fq.gz")
READS = ["1", "2"]

rule all:
   input:
       expand(config['datadirs']['qc'] + "/" + "{file}_{read}.fastqc.html", file = SAMPLES, read = READS)


rule fastqc:
   input:
      f1 = config['datadirs']['fastq'] + "/" + "{file}_{read}.fq.gz"
   output: config['datadirs']['qc'] + "/" + "{file}_{read}.fastqc.html", config['datadirs']['qc'] + "/" + "{file}_{read}.fastqc.zip"
   shell: """
      mkdir -p qc_output
      fastqc  --thread 8 --outdir qc_output --nogroup {input.f1}
          """
fastqc rna-seq snakemake • 354 views
ADD COMMENTlink modified 5 months ago by Biostar ♦♦ 20 • written 6 months ago by arshil0

Please show us the complete error message you receive.

ADD REPLYlink written 6 months ago by finswimmer13k

It's likely that there's a typo in the the output file name. Just double check that.

ADD REPLYlink written 6 months ago by Devon Ryan93k

thanks, yeah there was issue with the output naming.

ADD REPLYlink written 6 months ago by arshil0

Happens to all of us :)

ADD REPLYlink written 6 months ago by Devon Ryan93k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1019 users visited in the last hour