hi,
I'm working with ChIP-seq data and using MACS2 version- macs2 2.1.0.20140616. I run MACS2 callpeak twice for two different data sets to find the enrichment peaks for each one of them. I got 22,830 peaks at the first and 28,911 peaks at the second, when I use the bdg files for running bdgdiff it find only 6 and 183 differential peaks. why does MACS2 remove so many of the peaks? in which parameters should I focus on? thanks,
Nehora.