This is a test version of Biostars. For the public version, visit https://www.biostars.org.
missingoutputexception, snakemake after completing 95% analysis

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}
          """
rna-seq snakemake fastqc

Please show us the complete error message you receive.

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

thanks, yeah there was issue with the output naming.

Happens to all of us :)

0 answers

No answers yet.

Log in to answer this question.