Biostar Beta. Not for public use.
Peak center for ATAC-seq data
1
Entering edit mode
14 months ago
UCSF

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
1
Entering edit mode

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

ADD REPLYlink
1
Entering edit mode
22 months ago
igor 7.7k
United States

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
0
Entering edit mode

can i use chIp seeker for ATAC seq data?

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1