This is a test version of Biostars. For the public version, visit https://www.biostars.org.
call peak for all replicates in all conditions using macs2

Hi all,

I am doing ATAC-seq analysis, where I want to make a matrix of all consensus peaks from all replicates in all conditions. Can I just call all of them in one call using macs2?

macs2 callpeak --verbose 3 -t *.bed \
    -f BEDPE \
    -g ${macs_genomesize} -B \
    -q ${macs_q} \
    --extsize 200 --nomodel --shift -100 --nolambda \
    --keep-dup all --outdir ${peak_all_BEDPE_dir} -n "peak_all" --call-summits

or do I need to perform intersect for replicates and then merge for all conditions? what is a more properly to deal with multiple replicates and conditions

macs2 atacseq
• 1,851 views • link
updated by zau saa • zau saa • written by QX • QX

I would call peaks for individual replicates and conditions rather than pool them together. This will also help you check replicate agreement.

• 0 views • link
written by rfran010 • rfran010

after calling individually, how can you call the consensus peaks for all conditions? can you share the command for that?

• 0 views • link
written by QX • QX

Fully agree to rfran010's reply. And, compared to intersect peaks for all replicates, I think you may set consensus peaks for each condition to intervals which are more than 500bp and identified as peaks by at least 50% replicates.

• 0 views • link
written by zau saa • zau saa

do you know how to set this 500 and 50%? is that a macs2 option?

• 0 views • link
written by QX • QX

0 answers

No answers yet.

Log in to answer this question.

More posts like this