This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Fix SAM flags on quality filtered paired-end BAM file

Hello,

I have a BAM file of paired-end reads that have been quality filtered. I would like to salvage those reads whose mates have been filtered out and continue using them in my analysis (I would also appreciate any input on whether this is acceptable practice). Since some of my downstream tools give me errors if they see only one mate of a paired end read, is there a way to fix the SAM flags (to 0 or 16, for example) so that these reads are treated as single-end reads?

Thank you!

sam-flag bam paired-end-read
• 6,398 views • link
updated by Ram Ram • written by le2336 le2336

1 answer

I believe you can use the samtools fixmate command to update the flags, but you need to sort your BAM file by read name first:

samtools -n sort sample.bam sampled.sorted
samtools fixmate sample.sorted.bam sample.fixed.bam
• 0 views • link
updated by Ram Ram • written by James Ashmore James Ashmore

This works perfectly. Thank you so much!

• 0 views • link
written by le2336 le2336

Just for information, you should be able to use pipes to avoid intermediate files and have a final bam file coordinate sorted and mate fixed:

samtools sort -o -n aln.bam - \
| samtools fixmate -O bam - - \
| samtools sort - aln.sorted.fixed

I would like to salvage those reads whose mates have been filtered out and continue using them in my analysis (I would also appreciate any input on whether this is acceptable practice)

That's what I also do. When I filter for, say, mapping quality I don't mind if one mate fails and the other passes.

• 0 views • link
written by dariober dariober

This is useful to know. Thank you!

• 0 views • link
written by le2336 le2336

Log in to answer this question.

More posts like this