Skip to content

Commit

Permalink
remove bam to bed
Browse files Browse the repository at this point in the history
  • Loading branch information
lpantano committed Jul 15, 2024
1 parent cc5ca5d commit e6d9190
Showing 1 changed file with 1 addition and 20 deletions.
21 changes: 1 addition & 20 deletions subworkflows/local/bam_peaks_call_qc_annotate_macs2_homer.nf
Original file line number Diff line number Diff line change
Expand Up @@ -2,7 +2,6 @@
// Call peaks with MACS2, annotate with HOMER and perform downstream QC
//

include { BEDTOOLS_BAMTOBED } from '../../modules/nf-core/bedtools/bamtobed/main'
include { MACS2_CALLPEAK } from '../../modules/nf-core/macs2/callpeak/main'
include { HOMER_ANNOTATEPEAKS } from '../../modules/nf-core/homer/annotatepeaks/main'

Expand All @@ -29,29 +28,11 @@ workflow BAM_PEAKS_CALL_QC_ANNOTATE_MACS2_HOMER {

ch_versions = Channel.empty()

//
// Convert bam to bed
//
BEDTOOLS_BAMTOBED (
ch_bam.map { meta, ip_bam, control_bam -> [ meta, ip_bam ] }
)
ch_versions = ch_versions.mix(BEDTOOLS_BAMTOBED.out.versions.first())

// Create channels: [meta, ip_bed, []]
BEDTOOLS_BAMTOBED
.out
.bed
.map {
meta, ip_bed ->
[ meta, ip_bed, [] ]
}
.set { ch_bed }

//
// Call peaks with MACS2
//
MACS2_CALLPEAK (
ch_bed,
ch_bam,
macs_gsize
)
ch_versions = ch_versions.mix(MACS2_CALLPEAK.out.versions.first())
Expand Down

0 comments on commit e6d9190

Please sign in to comment.