MACS peak calling of paired-end reads in SAM format
1
0
Entering edit mode
3.2 years ago
biostart ▴ 360

HI,

Just a quick question: if SAM format is used for paired-end reads, will these be automatically recognised as paired-end by MACS when doing peak-calling?

Thanks!

ChIP-Seq MACS SAM • 1.4k views
ADD COMMENT
1
Entering edit mode
3.2 years ago
ATpoint 55k

MACS was developed in the days of short single-end reads. As far as I know, it can autodetect the file type (SAM, BAM, BED, Eland...), but you have to indicate if paired-end reads are present of not. I would write your SAM files as BAM (samtools view -b in.sam -o out.bam), and then add to your command -f BAMPE. You can also check from the stuff MACS writes to STDERR if MACS says "read treatment tags" (indicating single end) or "read treatment fragments" (indicating paired-reads, because two reads are read that are sequenced from one fragment).

ADD COMMENT

Login before adding your answer.

Traffic: 1781 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