Skip to content

Commit d5034fa

Browse files
committed
Follow up fix for #828 - name sort SE BAM too
1 parent 197059c commit d5034fa

File tree

1 file changed

+8
-8
lines changed

1 file changed

+8
-8
lines changed

workflows/rnaseq.nf

Lines changed: 8 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -330,8 +330,13 @@ workflow RNASEQ {
330330
ch_transcriptome_sorted_bam.join(ch_transcriptome_sorted_bai, by: [0])
331331
)
332332

333-
// Only name sort paired-end BAM files
334-
DEDUP_UMI_UMITOOLS_TRANSCRIPTOME
333+
// Name sort BAM before passing to Salmon
334+
SAMTOOLS_SORT (
335+
DEDUP_UMI_UMITOOLS_TRANSCRIPTOME.out.bam
336+
)
337+
338+
// Only run prepare_for_rsem.py on paired-end BAM files
339+
SAMTOOLS_SORT
335340
.out
336341
.bam
337342
.branch {
@@ -343,15 +348,10 @@ workflow RNASEQ {
343348
}
344349
.set { ch_umitools_dedup_bam }
345350

346-
// Name sort BAM before passing to Salmon
347-
SAMTOOLS_SORT (
348-
ch_umitools_dedup_bam.paired_end
349-
)
350-
351351
// Fix paired-end reads in name sorted BAM file
352352
// See: https://github.com/nf-core/rnaseq/issues/828
353353
UMITOOLS_PREPAREFORRSEM (
354-
SAMTOOLS_SORT.out.bam
354+
ch_umitools_dedup_bam.paired_end
355355
)
356356
ch_versions = ch_versions.mix(UMITOOLS_PREPAREFORRSEM.out.versions.first())
357357

0 commit comments

Comments
 (0)