Skip to content

Commit 9196038

Browse files
committed
Don't mix transcriptome bam stats with genome ones for multiqc
1 parent 1bc1b73 commit 9196038

File tree

1 file changed

+9
-16
lines changed
  • subworkflows/local/bam_dedup_umi

1 file changed

+9
-16
lines changed

subworkflows/local/bam_dedup_umi/main.nf

Lines changed: 9 additions & 16 deletions
Original file line numberDiff line numberDiff line change
@@ -101,21 +101,14 @@ workflow BAM_DEDUP_UMI {
101101
ch_dedup_transcriptome_bam = ended_transcriptome_dedup_bam.single_end
102102
.mix(UMITOOLS_PREPAREFORRSEM.out.bam)
103103

104-
// Collect files useful for MultiQC into one helpful emission
105-
106-
ch_stats = UMI_DEDUP_GENOME.out.stats
107-
.mix(UMI_DEDUP_TRANSCRIPTOME.out.stats)
108-
109-
ch_flagstat = UMI_DEDUP_GENOME.out.flagstat
110-
.mix(UMI_DEDUP_TRANSCRIPTOME.out.flagstat)
111-
112-
ch_idxstats = UMI_DEDUP_GENOME.out.idxstats
113-
.mix(UMI_DEDUP_TRANSCRIPTOME.out.idxstats)
104+
// Collect files useful for MultiQC into one helpful emission. Don't
105+
// automatically add transcriptome stats- difficult to separate in multiqc
106+
// without a bit more work
114107

115108
ch_multiqc_files = ch_dedup_log
116-
.mix(ch_stats)
117-
.mix(ch_flagstat)
118-
.mix(ch_idxstats)
109+
.mix(UMI_DEDUP_GENOME.out.stats)
110+
.mix(UMI_DEDUP_GENOME.out.flagstat)
111+
.mix(UMI_DEDUP_GENOME.out.idxstats)
119112
.transpose()
120113
.map{it[1]}
121114

@@ -129,9 +122,9 @@ workflow BAM_DEDUP_UMI {
129122
bam = UMI_DEDUP_GENOME.out.bam // channel: [ val(meta), path(bam) ]
130123
bai = bam_csi_index ? UMI_DEDUP_GENOME.out.csi : UMI_DEDUP_GENOME.out.bai // channel: [ val(meta), path(bai) ]
131124
dedup_log = ch_dedup_log // channel: [ val(meta), path(log) ]
132-
stats = ch_stats
133-
flagstat = ch_flagstat
134-
idxstats = ch_idxstats
125+
stats = UMI_DEDUP_GENOME.out.stats.mix(UMI_DEDUP_TRANSCRIPTOME.out.stats)
126+
flagstat = UMI_DEDUP_GENOME.out.flagstat.mix(UMI_DEDUP_TRANSCRIPTOME.out.flagstat)
127+
idxstats = UMI_DEDUP_GENOME.out.idxstats.mix(UMI_DEDUP_TRANSCRIPTOME.out.idxstats)
135128
multiqc_files = ch_multiqc_files
136129
transcriptome_bam = ch_dedup_transcriptome_bam // channel: [ val(meta), path(bam) ]
137130
versions = ch_versions // channel: [ path(versions.yml) ]

0 commit comments

Comments
 (0)