This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Sam to bowtie native format

Hi everyone,

Does anyone know how to convert .sam files to the native bowtie format? I need this for peak calling using QuEST software for ChIP seq data.

Thanks a bunch!

teabonng

chip-seq quest
• 1,977 views • link
updated by h.mon h.mon • written by druggable • druggable

What is native bowtie format? Do you have a reference? If you are referring to the alignments written by bowtie then they are in SAM format, which is a standard.

• 0 views • link
written by GenoMax GenoMax

Thanks for the replies, now I understand.

• 1 views • link
written by druggable • druggable

Just out of curiosity, why do you want to use an almost 10-year old peak caller?

• 0 views • link
written by ATpoint ATpoint

1 answer

QuEST accepts sam files as inputs, which is Bowtie's output if you use the -S flag.

How to use QuEST for ATAC-Seq data

• 0 views • link
written by h.mon h.mon

From this link:

For a 4.25 GB ATAC Sam file, it ran for 30 hours.

Well, this even more should motivate you to use a more recent peak caller. For ATAC-seq, I typically use MACS2:

macs2 callpeak -t file.bam -n samplename --nomodel -f BAMPE

It assumes that your file is paired-end.

• 0 views • link
written by ATpoint ATpoint

Log in to answer this question.

More posts like this