BWA mem vs. Sampe alignment
Entering edit mode
5.2 years ago
kirannbishwa01 ★ 1.3k

I have done BWA alignment before using - mem but need to rethink over it again. My main objective is to call polymorphisms from my genome resequenced data for two different populations for a plant model which has its reference genome sequenced.

I have illumina PE reads but after trimming the reads now range from 50 to 101 bp but majority above 80 bp.

1) why/which, -mem vs. -sampe would be the best choice for calling polymorphisms using the data I have?

2) I have done bwa mem with - M which flags the short split hits as secondary. why do we need this flag if we are aligning genomic reseq data? isn't it for data like RNAseq. If it can be used for genomic reads, how is it going to be useful.


sequencing alignment snp BWA • 2.2k views

Login before adding your answer.

Traffic: 1542 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6