This is a test version of Biostars. For the public version, visit https://www.biostars.org.
MACS2 peak calling

Hi all!

I have performed peak calling using MACS2. Is it correct that the bed file has not strand of the peak?

Thank you very much for your help

macs2 peak-calling
• 1,893 views • link
updated by seidel seidel • written by Maurice • Maurice

1 answer

Yes, that is correct. Peaks from experiments such as ChIP-seq or ATAC-seq are based on DNA, which is double-stranded and unlike RNA (where we have expression from a specific strand) not strand-specific.

• 0 views • link
written by ATpoint ATpoint

Thank you! And for what regards peaks annotation, how can the tool I'm going to use assign peaks to the correct location without knowing the strand of the peak?

• 0 views • link
written by Maurice • Maurice

You do not need strand information for this as, as I said, strands do not apply in DNA-seq.

• 0 views • link
written by ATpoint ATpoint

What ATpoint said. But given that, getting the distance between features depends on the tool you're using and how you ask questions. The "correct location" is simply a point on a numbered line of bases, regardless of strand. If you want to ask questions about distance between your peak and other features, make sure you know how your tool compares stranded and nonstranded features.

• 0 views • link
written by seidel seidel

Log in to answer this question.

More posts like this