SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



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

Reply
 
Thread Tools
Old 07-15-2011, 07:19 AM   #1
oiiio
Senior Member
 
Location: USA

Join Date: Jan 2011
Posts: 105
Default Question about BWA mapping quality

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.
oiiio is offline   Reply With Quote
Old 07-18-2011, 08:58 AM   #2
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Check out the "pairing" function in "bwape.c".
nilshomer is offline   Reply With Quote
Old 07-19-2011, 12:50 AM   #3
volks
Member
 
Location: hd.de

Join Date: Jun 2010
Posts: 81
Red face

Quote:
Originally Posted by nilshomer View Post
Check out the "pairing" function in "bwape.c".
which is quite hard to understand. can anyone translate?
volks is offline   Reply With Quote
Old 07-19-2011, 07:05 AM   #4
oiiio
Senior Member
 
Location: USA

Join Date: Jan 2011
Posts: 105
Default

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.
oiiio is offline   Reply With Quote
Old 07-21-2011, 07:23 AM   #5
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Quote:
Originally Posted by volks View Post
which is quite hard to understand. can anyone translate?
Your local bioinformatician would be happy to help and receive the business.
nilshomer is offline   Reply With Quote
Old 07-25-2011, 11:23 AM   #6
dp05yk
Member
 
Location: Brock University

Join Date: Dec 2010
Posts: 66
Default

Quote:
Originally Posted by nilshomer View Post
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.
dp05yk is offline   Reply With Quote
Old 07-25-2011, 04:33 PM   #7
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Quote:
Originally Posted by dp05yk View Post
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?
nilshomer is offline   Reply With Quote
Reply

Tags
bwa, mapping quality, mapping quality scores

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:57 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO