I am creating a coverage plot +/- 250bp around a transcription start site with various ATAC-seq replicates. I have created these plots very easily for S2 (drosophila embryonic) cells using microarray data. However, I am having a hard time locating precise peak centers in the ATAC data which is very noisy. Ideally, I would like a smooth curve centered at 0bp. Right now, my plots seem to be really skewed right or left depending on the replicate (I am actually finding the peak centers by hand, using a BW visualization file which is clearly not working nor is that efficient) Is there any better, consistent, way of doing this? I am calling peaks using HOMER and creating the coverage plot using perl and python scripts.
Question: Peak center for ATAC-seq data
2.3 years ago by
christacaggiano • 20
christacaggiano • 20 wrote:
ADD COMMENT • link •
2.2 years ago by
igor ♦ 7.3k
igor ♦ 7.3k wrote:
This sounds like two separate questions. Do you want to make TSS enrichment plots or find peak centers?
To visualize enrichment at TSS, you can try these methods:
- HOMER (which you are already using): http://homer.salk.edu/homer/ngs/quantification.html
- ChIPseeker: http://guangchuangyu.github.io/2014/04/visualization-methods-in-chipseeker/
- deepTools: http://deeptools.readthedocs.io/en/latest/content/tools/plotProfile.html
ADD COMMENT • link
Please log in to add an answer.
Powered by Biostar version 2.3.0
Traffic: 1267 users visited in the last hour