Question: MACS peak calling of paired-end reads in SAM format
0
gravatar for biostart
11 months ago by
biostart300
Germany
biostart300 wrote:

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!

sam chip-seq macs • 563 views
ADD COMMENTlink written 11 months ago by biostart300
1
gravatar for ATpoint
11 months ago by
ATpoint19k
Germany
ATpoint19k 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 11 months ago • written 11 months ago by ATpoint19k
Please log in to add an answer.

Help
Access

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