This is a test version of Biostars. For the public version, visit https://www.biostars.org.
bigwig to peak

Hi,

I've a bigwig containing a chip-seq signal for a histone mark (H3K4me1) and I want to produce a bed file containing the peaks (the same way MACS produced a bigwig and a peak bed file)? How can I do that? If it's impossible I can realign the raw data and do the whole peak calling step but I prefer to avoid that.

Thanks

peak bigwig chip-seq
• 26,790 views • link
updated by Ram Ram • written by Nicolas Rosewick Nicolas Rosewick

I ran into the same goal. However, I believe this is not working, as the bigwig might not be properly normalized the way that MACS does it

• 1 views • link
written by liruiradiant • liruiradiant

1 answer

MACS2 can call peak from bedgraph, you can convert bigwig to bedgraph first by bigwigtobedgraph from here and then call peak with MACS2.

bdgpeakcall: Call peaks from bedGraph output.
• 2 views • link
updated by Ram Ram • written by Ming Tommy Tang Ming Tommy Tang

Log in to answer this question.

More posts like this