Question about BWA mapping quality

oiiio · 07-15-2011, 08:19 AM

I need to verify that BWA takes into account mate pair data, and the BWA paper that I found seems to be somewhat ambiguous in how it derives mapping quality for these paired reads. If someone could point me at some source code or describe how BWA deals with paired data that would be extremely helpful.

Also, I would feel much happier if I knew how the actual implementation of this statement from the paper worked:

Quote:

BWA assigns a Phred-like quality score to each read (MapQ) based on match uniqueness, sequence identity, end-pairing, and insert size that is intended to indicate confidence of read placement accuracy in the genome

Again if someone could give a little bit of a more detailed explanation or point me at the source code file (or both) I would be very happy.

nilshomer · 07-18-2011, 09:58 AM

Check out the "pairing" function in "bwape.c".

volks · 07-19-2011, 01:50 AM

Quote:

Originally Posted by nilshomer

Check out the "pairing" function in "bwape.c".

which is quite hard to understand. can anyone translate?

oiiio · 07-19-2011, 08:05 AM

Thanks for the pointer Nils, I'll give it a look. And of course, I'd definitely like to hear any other input people may have as well.

nilshomer · 07-21-2011, 08:23 AM

Quote:

Originally Posted by volks

which is quite hard to understand. can anyone translate?

Your local bioinformatician would be happy to help and receive the business.

dp05yk · 07-25-2011, 12:23 PM

Quote:

Originally Posted by nilshomer

Your local bioinformatician would be happy to help and receive the business.

To his credit, the code is kind of a mess, even if the theory behind it is solid.

nilshomer · 07-25-2011, 05:33 PM

Quote:

Originally Posted by dp05yk

To his credit, the code is kind of a mess, even if the theory behind it is solid.

That may be so, but then whos mess it to clean up?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
question about the quality when mapping with gmap	lixiangru	RNA Sequencing	1	01-06-2015 04:47 PM
negative bwa mapping quality	Anney	Bioinformatics	4	07-11-2011 02:42 PM
BWA mapping quality scores?	kweber2	Genomic Resequencing	2	09-27-2010 04:01 PM
Interpretation of BWA mapping quality	christophpale	Bioinformatics	0	07-21-2010 04:15 AM
bwa mapping quality	totalnew	Bioinformatics	6	05-21-2010 05:50 AM

07-18-2011, 09:58 AM	#2
nilshomer Nils Homer Location: Boston, MA, USA Join Date: Nov 2008 Posts: 1,285	Check out the "pairing" function in "bwape.c".

07-19-2011, 08:05 AM	#4
oiiio Senior Member Location: USA Join Date: Jan 2011 Posts: 105	Thanks for the pointer Nils, I'll give it a look. And of course, I'd definitely like to hear any other input people may have as well.