Question: What is the difference between bed files from bigWig or bam, and bed files of narrow/broad peaks?
gravatar for 1888
7 months ago by
United States
188860 wrote:


Can someone point to procedures on how to get a bed file from a bigWig file, that somewhat reflects broad/narrow peak regions called in the same cell e.g. from Encode (if this is possible)?

I know it is possible to get a bed file from a bigWig file in Encode for example using the utility 'bigWigToBedGraph', but the resulting bed regions are very different from the regions in the bed files in broad/narrow peaks - is there a way to converge to similar regions?

Thank you!

alignment bed • 267 views
ADD COMMENTlink modified 7 months ago by Devon Ryan97k • written 7 months ago by 188860
gravatar for Devon Ryan
7 months ago by
Devon Ryan97k
Freiburg, Germany
Devon Ryan97k wrote:

bigWig files are generally whole-genome coverage tracks, so you would need to call peaks on them. You can convert them to bedGraph format and then input that into MACS2 bdgpeakcall. Note however that you may have to manually play around with the parameters quite a bit to get reasonable peak calls out of that. In general it's preferable to try and find the original BAM files from which the bigWig files were generated and just give those to a peak caller, since that's what the peak callers are written to work with.

ADD COMMENTlink written 7 months ago by Devon Ryan97k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1496 users visited in the last hour