Question: Difference in number of peaks called
gravatar for anu014
22 months ago by
anu014160 wrote:

Hi Biostars,

I have used macs2 for calling peaks in chip-seq data. Firstly, I called peaks for antibiotic-treated sample & wild-type sample individually. After I ran macs2 with antibiotic-treated sample as treatment group (-t) & wild-type sample as control group (-c). Lastly, I ran the tool using antibiotic-treated sample as control group (-c) & wild-type sample as treatment group (-t). I got following number of peaks:

Antibiotic-treated S :  1007
WT S: 100
Antibiotic-treated vs WT : 46
WT vs Antibiotic-treated : 90

I am not getting the stats. If 100 were picked up individually & 90 peaks are called in WT vs Antibiotic treated run, then in Antibiotic-treated vs WT , the remaining 10 peaks should be checked against peaks present in Antibiotic-treated group. Why am I getting 46 peaks?

command parameters: # Command line: callpeak -f BAMPE -g 6748384 -B --mfold 1 60 --broad --broad-cutoff 0.1 --bw 300 --buffer-size 100000 --qvalue 0.05 --slocal 1000 --llocal 10000 --cutoff-analysis --fe-cutoff 1.0

Thank you in advance! :)

peakcalling chip-seq macs2 • 803 views
ADD COMMENTlink modified 12 months ago by colindaven1.2k • written 22 months ago by anu014160
gravatar for colindaven
12 months ago by
Hannover Medical School
colindaven1.2k wrote:

Shouldn't you be calling each against the Input (same chip-seq process, but not applying the antibody) for these data, not against WT or control ? Have you looked at the plotfingerprints method from Deeptools - in Galaxy ?

ADD COMMENTlink written 12 months ago by colindaven1.2k
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: 696 users visited in the last hour