This is a test version of Biostars. For the public version, visit https://www.biostars.org.
MACS for a paired-end sequencing

Hello, I need some help with peak calling with MACS for a paired-end sequencing generated fastq files. I have separate R1/R2 .bam files and would like to know the best way to proceed further.

chip-seq
• 2,065 views • link
written by KT • KT

1 answer

If it is paired-end fastq file then why you have separate bam files?. The best would be to realign the paired-end data and you will have one bam file and that can be used for peak calling.

• 0 views • link
written by Prakash Prakash

Agreed. PE data have to be mapped together otherwise the advantage of the method is gone and you would be double-counting every read. Once you have a paired-end BAM file just use -f BAMPE as parameter in macs. It will take care of the rest.

• 0 views • link
written by ATpoint ATpoint

Thank you very much for suggestion. I will try that and update accordingly.

• 0 views • link
written by KT • KT

Log in to answer this question.

More posts like this