@@ -30,12 +30,12 @@ workflow DEDUPLICATE {
3030 }
3131
3232 if ( params. deduplication_mergejobs ) {
33-
33+
3434 // No splitting of .bam files by contig, deduplicate all in one
3535 input_for_deduplication = ch_bam_bai
3636
3737 } else {
38-
38+
3939 // Create genomic regions file for splitting the bam before deduplication
4040 BUILD_INTERVALS ( fasta_fai )
4141 ch_versions = ch_versions. mix( BUILD_INTERVALS . out. versions. first() )
@@ -93,9 +93,9 @@ workflow DEDUPLICATE {
9393 ch_markduplicates_input. fasta,
9494 ch_markduplicates_input. fasta_fai
9595 )
96- ch_versions = ch_versions. mix( PICARD_MARKDUPLICATES . out. versions. first() )
96+ ch_versions = ch_versions. mix( PICARD_MARKDUPLICATES . out. versions. first() )
9797
98- ch_dedupped_bam = PICARD_MARKDUPLICATES . out. bam
98+ ch_dedupped_bam = PICARD_MARKDUPLICATES . out. bam
9999
100100 } else if ( params. deduplication_tool == " dedup" ) {
101101 ch_dedup_input = input_for_deduplication
@@ -105,7 +105,7 @@ workflow DEDUPLICATE {
105105 }
106106
107107 DEDUP ( ch_dedup_input )
108- ch_versions = ch_versions. mix( DEDUP . out. versions. first() )
108+ ch_versions = ch_versions. mix( DEDUP . out. versions. first() )
109109
110110 ch_dedupped_bam = DEDUP . out. bam
111111 }
@@ -147,29 +147,29 @@ workflow DEDUPLICATE {
147147 ch_input_for_samtools_merge. fasta,
148148 ch_input_for_samtools_merge. fasta_fai
149149 )
150- ch_versions = ch_versions. mix( SAMTOOLS_MERGE_DEDUPPED . out. versions )
150+ ch_versions = ch_versions. mix( SAMTOOLS_MERGE_DEDUPPED . out. versions )
151151
152152 ch_input_for_samtools_sort_dedupped = SAMTOOLS_MERGE_DEDUPPED . out. bam
153153
154154 }
155155
156156
157-
158157 // Sort the merged bam and index
159158 SAMTOOLS_SORT_DEDUPPED ( ch_input_for_samtools_sort_dedupped )
160159 ch_versions = ch_versions. mix( SAMTOOLS_SORT_DEDUPPED . out. versions )
161160 ch_dedup_bam = SAMTOOLS_SORT_DEDUPPED . out. bam
162161
163162 SAMTOOLS_INDEX_DEDUPPED ( ch_dedup_bam )
164163 ch_versions = ch_versions. mix( SAMTOOLS_INDEX_DEDUPPED . out. versions )
165- ch_dedup_bai = params. fasta_largeref ? SAMTOOLS_INDEX_DEDUPPED . out. csi : SAMTOOLS_INDEX_DEDUPPED . out. bai
164+ ch_dedup_bai = params. fasta_largeref ? SAMTOOLS_INDEX_DEDUPPED . out. csi : SAMTOOLS_INDEX_DEDUPPED . out. bai
166165
167166 // Finally run flagstat on the dedupped bam
168167 ch_input_for_samtools_flagstat = ch_dedup_bam. join( ch_dedup_bai )
169168
170169 SAMTOOLS_FLAGSTAT_DEDUPPED (
171170 ch_input_for_samtools_flagstat
172171 )
172+
173173 ch_versions = ch_versions. mix( SAMTOOLS_FLAGSTAT_DEDUPPED . out. versions )
174174 ch_multiqc_files = ch_multiqc_files. mix( SAMTOOLS_FLAGSTAT_DEDUPPED . out. flagstat )
175175 ch_dedup_flagstat = SAMTOOLS_FLAGSTAT_DEDUPPED . out. flagstat
0 commit comments