This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Alignment with BWA: how to discard multi mappers?

Hello,

This topic has been discussed on many occasions however I couldn't find the right answer afterall.

I am aligning Illumina generated paired-end reads on reference genome with BWA. My understanding is that when BWA finds multiple equally perfect alignments for one read, it will report/map that read randomly. However, I couldn't find an option for discarding these reads, i.e. similar to -m1 option from bowtie (there is a specific reason why I want to use bwa over bowtie).

Did I miss something or...?

Thanks

bwa paired-end
• 4,710 views • link
updated by igor igor • written by ThePresident ThePresident

You can screen the bam file afterwards for reads with no other hits.

• 1 views • link
written by Asaf Asaf

1 answer

There is no parameter you can use, but you can filter the BAM with samtools. See this previous discussion: BWA mem how to know if a read is mapped uniquely? and Bwa-Mem: Discriminate Between Reads Mapping Uniquely And Those Mapping In Multiple Positions

• 1 views • link
written by igor igor

Log in to answer this question.

More posts like this