MACS for a paired-end sequencing
1
0
Entering edit mode
2.5 years ago
KT • 0

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 • 726 views
ADD COMMENT
3
Entering edit mode
2.5 years ago
Prakash ★ 2.2k

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.

ADD COMMENT
0
Entering edit mode

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.

ADD REPLY
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 1226 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6