Question: MACS peak calling of paired-end reads in SAM format
gravatar for biostart
22 months ago by
biostart340 wrote:


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?


sam chip-seq macs • 895 views
ADD COMMENTlink written 22 months ago by biostart340
gravatar for ATpoint
22 months ago by
ATpoint34k wrote:

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 COMMENTlink modified 22 months ago • written 22 months ago by ATpoint34k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1001 users visited in the last hour