missingoutputexception, snakemake after completing 95% analysis
0
0
Entering edit mode
4.8 years ago
arshil • 0

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 • 1.4k views
ADD COMMENT
0
Entering edit mode

Please show us the complete error message you receive.

ADD REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

thanks, yeah there was issue with the output naming.

ADD REPLY
0
Entering edit mode

Happens to all of us :)

ADD REPLY

Login before adding your answer.

Traffic: 2534 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