Integrate output of Rose into workflow and disable eHMM

daisybio · Jan 9, 2024 · 6aa60a6 · 6aa60a6
1 parent 7b7f492
commit 6aa60a6
Showing 1 changed file with 10 additions and 1 deletion.
diff --git a/subworkflows/local/peaks.nf b/subworkflows/local/peaks.nf
@@ -71,6 +71,7 @@ workflow PEAKS {
         REMOVE_CHR(ch_footprints, [])
         ch_footprints_nochr = REMOVE_CHR.out.output
 
+        /*
         if (params.bam_design) {
             ch_bams = Channel.value(file(params.bam_design))
                                 .splitCsv(header: true)
@@ -101,7 +102,8 @@ workflow PEAKS {
 
             ch_footprints = EHMM.out.all
         }
-
+        */
+
         if (params.bam_design2) {
             // Parse design file
             ch_bams = Channel.value(file(params.bam_design2))
@@ -139,6 +141,13 @@ workflow PEAKS {
             params.ucsc_file
         )
 
+        // Adapt process output to previous design of ch_footprints
+        ch_footprints = ROSE.out
+            .map{meta, bed -> 
+                meta = meta + [id: meta.state + "_enhancers", antibody: "enhancers"]
+                [meta, bed]
+            }
+
         }
 
         if (params.blacklist) {