Question: Peak center for ATAC-seq data
gravatar for christacaggiano
2.3 years ago by
christacaggiano20 wrote:

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.

ADD COMMENTlink modified 2.2 years ago by igor7.3k • written 2.3 years ago by christacaggiano20

If you call peaks using macs2, it also give the summit (--call-summit) , which is the peak center.

ADD REPLYlink modified 24 months ago • written 2.2 years ago by geek_y9.1k
gravatar for igor
2.2 years ago by
United States
igor7.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:

ADD COMMENTlink modified 2.2 years ago • written 2.2 years ago by igor7.3k

can i use chIp seeker for ATAC seq data?

ADD REPLYlink written 19 months ago by krushnach80460
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: 1267 users visited in the last hour