Hi,
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!
•
written
by
1888
1 answer
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.
Log in to answer this question.
-
Overestimate of ATAC-seq differential peaks
written by Darryl •
-
Define enhancer regions
written by Faiza •
-
macs2 bdgpeakcall only gives a few peaks
written by baijiangshan9726 •
-
Can I successfully go from a bed file (narrowpeak) back to a bam file containing the mapped ATAC-se…
written by Jillian •
-
How to get Read Counts from MACS2 output files?
written by cthangav •
-
Differential Peak Analysis from ENCODE Pipeline normalized outputs
written by sjquon •
-
WGBS - Methylation state at CpG to peak calls
written by Nico Gerstner •
-
Peak calling for broad histone-modification regions
written by biostart
-
How do I view the height / distribution of ChIP-seq peaks?
written by James Ashmore
-
bigwig : peak distance from specific genomic region
written by Nicolas Rosewick
