Biostar Beta. Not for public use.
bigwig to peak
1
Entering edit mode
10 months ago
Belgium, Brussels

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

ADD COMMENTlink
0
Entering edit mode

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

ADD REPLYlink
5
Entering edit mode
13 months ago
Ming Tang ♦ 2.5k
Houston/MD Anderson Cancer Center

MACS2 can call peak from bedgraph, you can convert bigwig to bedgraph first by bigwigtobedgraph from here http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/ and then call peak with MACS2 https://github.com/taoliu/MACS/

bdgpeakcall: | Call peaks from bedGraph output.
---|---

ADD COMMENTlink
0
Entering edit mode

thanks that's perfect

ADD REPLYlink
0
Entering edit mode

thanks that's perfect

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1