Use read counts to compute scaling factor for p-value bedGraph files - fix #7

emi80 · emi80 · commit cbc150394e8c · 2021-12-24T10:09:45.000+01:00
diff --git a/chipseq-pipeline.nf b/chipseq-pipeline.nf
@@ -444,7 +444,7 @@ process narrowPeakCallNoInput {
 }
 
 crossedBams.map{ c, t ->
-  [t[0], t[4], c[4]]
+  [t[0], t[6], c[6]]
 }.cross(pileupBedGraphFiles)
 .map { r, s ->
   def (treat, control) = r[1..-1] as long[]

Original file line number	Diff line number	Diff line change
`@@ -444,7 +444,7 @@ process narrowPeakCallNoInput {`
`444`	`444`	`}`
`445`	`445`
`446`	`446`	`crossedBams.map{ c, t ->`
`447`		`- [t[0], t[4], c[4]]`
	`447`	`+ [t[0], t[6], c[6]]`
`448`	`448`	`}.cross(pileupBedGraphFiles)`
`449`	`449`	`.map { r, s ->`
`450`	`450`	`def (treat, control) = r[1..-1] as long[]`