diff --git a/src/RNAseq.config b/src/RNAseq.config index c8f39d48b40d19ea1fc5df10aeea62956ddc4abc..d15b85fde860b15260633dbcbf30de0d19a7fca8 100644 --- a/src/RNAseq.config +++ b/src/RNAseq.config @@ -104,6 +104,10 @@ profiles { container = "lbmc/htseq:0.11.2" cpus = 1 } + withName: multiqc { + container = "ewels/multiqc:1.9" + cpus = 1 + } } } } diff --git a/src/RNAseq.nf b/src/RNAseq.nf index 8022c44a2579f2c0b36b7e2f8466d401dcbdf910..e96e2fe590e5173ef02d6e199efb6612e2506c34 100644 --- a/src/RNAseq.nf +++ b/src/RNAseq.nf @@ -45,7 +45,7 @@ process fastqc_raw { set file_id, file(reads) from INPUT_FASTQC_RAW output: - set file_id, "*.{html,zip}" into OUTPUT_FASTQC_RAW + file "*_fastqc.{html,zip}" into OUTPUT_FASTQC_RAW when: params.do_fastqc @@ -104,7 +104,7 @@ process fastqc_cut { set file_id, file(reads) from CUTADAPT_OUTPUT_FASTQC output: - set file_id, "*.{html,zip}" into OUTPUT_FASTQC_CUT + file "*.{html,zip}" into OUTPUT_FASTQC_CUT when: params.do_fastqc @@ -265,6 +265,9 @@ if (! params.do_dedup) { HISAT_ALIGNED_DEDUP.into{DEDUP_GENOME} } + + + /* HTseq process sort_bam { @@ -349,3 +352,30 @@ bash ${script} -b ${bam} \ """ } */ + +/* MultiQC */ + + +process multiqc { + tag "multiqc" + publishDir "${params.output}/multiqc", mode: 'copy' + + input: + file ('fastqc/*') from OUTPUT_FASTQC_RAW.collect().ifEmpty([]) + file ('*') from CUTADAPT_LOG.collect().ifEmpty([]) + file ('fastqc/*') from OUTPUT_FASTQC_CUT.collect().ifEmpty([]) + file ('*') from FILTER_LOG.collect().ifEmpty([]) + file ('fastqc/*') from OUTPUT_FASTQC_FILTER.collect().ifEmpty([]) + file ('*') from HISAT_LOG.collect().ifEmpty([]) + file ('fastqc/*') from OUTPUT_FASTQC_GENOME.collect().ifEmpty([]) + + + output: + file "multiqc_report.html" into multiqc_report + file "multiqc_data" + + script: +""" +multiqc ./ +""" +}