Skip to content

Commit 188d9ff

Browse files
committed
fixed a sorting bug in make_bigwig_for_cleavage_sites
1 parent 38c55a3 commit 188d9ff

1 file changed

Lines changed: 2 additions & 2 deletions

File tree

main.nf

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -794,8 +794,8 @@ process make_bigwig_for_cleavage_sites {
794794
--tag ${params.tag_quantification_cs} \\
795795
--include_multimappers ${params.include_multimappers}
796796
797-
# 3. Sort the extracted BED
798-
sort -k1,1 -k2,2n cs.bed > cs.sorted.bed
797+
# 3. Sort the extracted BED by chr, start, end, AND strand (-k6,6)
798+
sort -k1,1 -k2,2n -k3,3n -k6,6 cs.bed > cs.sorted.bed
799799
800800
# 4. Generate Weighted BedGraphs per strand
801801
# Since intervals are exactly 1bp long, we just group by coordinate and sum the weights (col 5).

0 commit comments

Comments
 (0)