MACS2 peak caller omitting large chunks of peaks
0
0
Entering edit mode
11 weeks ago

Hi all,

We performed ATAC-seq on a cell line for which we got about 100,000 peaks using MACS2 peak caller. However, in some of the chromosomes there were very few or no peaks called at all. Is this normal? Would we reasonable expect that there would be similar number of peaks called across the different chromosomes. What could potentially be wrong?

Many thanks in advance!

peak-calling ATAC-seq MACS2 • 167 views
ADD COMMENT
2
Entering edit mode

Given the different size of chromosomes (in human chr22 is about 1/5th of the size of chr1), I do not see how the assumption of equal peak numbers per chromosome would hold true, or why this would be assumed in the first place.

So what is the issue here? Do you see peaks on the IGV that macs2 does not call or are you more like wondering why some regions or chromosomes have few peaks? Can you add numbers to back that up?

ADD REPLY

Login before adding your answer.

Traffic: 1472 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6