@@ -69,13 +69,6 @@ include { FASTQ_QC_TRIM_FILTER_SETSTRANDEDNESS } from '../../subwor
6969~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
7070*/
7171
72- // Header files for MultiQC
73- ch_pca_header_multiqc = file(" $projectDir /workflows/rnaseq/assets/multiqc/deseq2_pca_header.txt" , checkIfExists : true )
74- sample_status_header_multiqc = file(" $projectDir /workflows/rnaseq/assets/multiqc/sample_status_header.txt" , checkIfExists : true )
75- ch_clustering_header_multiqc = file(" $projectDir /workflows/rnaseq/assets/multiqc/deseq2_clustering_header.txt" , checkIfExists : true )
76- ch_biotypes_header_multiqc = file(" $projectDir /workflows/rnaseq/assets/multiqc/biotypes_header.txt" , checkIfExists : true )
77- ch_dummy_file = ch_pca_header_multiqc
78-
7972workflow RNASEQ {
8073
8174 take :
@@ -99,6 +92,13 @@ workflow RNASEQ {
9992
10093 main :
10194
95+ // Header files for MultiQC
96+ ch_pca_header_multiqc = file(" $projectDir /workflows/rnaseq/assets/multiqc/deseq2_pca_header.txt" , checkIfExists : true )
97+ sample_status_header_multiqc = file(" $projectDir /workflows/rnaseq/assets/multiqc/sample_status_header.txt" , checkIfExists : true )
98+ ch_clustering_header_multiqc = file(" $projectDir /workflows/rnaseq/assets/multiqc/deseq2_clustering_header.txt" , checkIfExists : true )
99+ ch_biotypes_header_multiqc = file(" $projectDir /workflows/rnaseq/assets/multiqc/biotypes_header.txt" , checkIfExists : true )
100+ ch_dummy_file = ch_pca_header_multiqc
101+
102102 ch_multiqc_files = Channel . empty()
103103 ch_trim_status = Channel . empty()
104104 ch_map_status = Channel . empty()
@@ -544,7 +544,7 @@ workflow RNASEQ {
544544 // Get RSeqC modules to run
545545 def rseqc_modules = params. rseqc_modules ? params. rseqc_modules. split(' ,' ). collect{ it. trim(). toLowerCase() } : []
546546 if (params. bam_csi_index) {
547- for (rseqc_module in [' read_distribution' , ' inner_distance' , ' tin' ]) {
547+ [' read_distribution' , ' inner_distance' , ' tin' ]. each { rseqc_module ->
548548 if (rseqc_modules. contains(rseqc_module)) {
549549 rseqc_modules. remove(rseqc_module)
550550 }
@@ -571,12 +571,12 @@ workflow RNASEQ {
571571 .map {
572572 meta, strand_log ->
573573 def rseqc_inferred_strand = getInferexperimentStrandedness(strand_log, params. stranded_threshold, params. unstranded_threshold)
574- rseqc_strandedness = rseqc_inferred_strand. inferred_strandedness
574+ def rseqc_strandedness = rseqc_inferred_strand. inferred_strandedness
575575
576576 def status = ' fail'
577577 def multiqc_lines = []
578578 if (meta. salmon_strand_analysis) {
579- salmon_strandedness = meta. salmon_strand_analysis. inferred_strandedness
579+ def salmon_strandedness = meta. salmon_strand_analysis. inferred_strandedness
580580
581581 if (salmon_strandedness == rseqc_strandedness && rseqc_strandedness != ' undetermined' ) {
582582 status = ' pass'
0 commit comments