SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics
Similar Threads
Thread Thread Starter Forum Replies Last Post
Can paired-end mapping produce more reads than single-end ? warrenemmett Bioinformatics 13 03-20-2012 11:10 PM
Question about BWA mapping quality oiiio Bioinformatics 6 07-25-2011 04:33 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 03-27-2010, 08:51 AM   #1
pparg
Member
 
Location: NY

Join Date: Aug 2008
Posts: 19
Default BWA paired end mapping quality

I used BWA to map my PE sequencing data to reference genome. I try to use paired mapping quality to filter bad read pairs out for downstream analysis.
How BWA calculate paired mapping quality? I understand it calculates single-end mapping quality like MAQ does. But I am not sure how it proceeds after having the single mapping quality for both ends? Simply add up or something more complicated? I’ve checked the source code, but the program does not make much sense without a good understanding of the variable names/notations. FYI, the relevant source code is located in the ‘static int pairing’ function of the bwape.c file.
I would really appreciate your input.
pparg
pparg is offline   Reply With Quote
Old 03-29-2010, 10:34 AM   #2
pparg
Member
 
Location: NY

Join Date: Aug 2008
Posts: 19
Default

Hello, does anybody have any ideas on this? Thank you!
pparg is offline   Reply With Quote
Old 05-16-2011, 12:24 AM   #3
Seq84
Member
 
Location: Italy

Join Date: Feb 2011
Posts: 19
Default

Hi all, I'm interested too! Could someone post a link or a brief description of BWA quality mapping scoring ?

Thanks in advance.
Seq84 is offline   Reply With Quote
Old 06-09-2011, 10:56 PM   #4
Azazel
Member
 
Location: Japan

Join Date: Oct 2010
Posts: 52
Default +1, I have the exact same question, too

+1, I have the exact same question, too

I'd also like to know how the mapping quality for paired end reads is computed, is it just the sum of the quality of the two separate reads?
Azazel is offline   Reply With Quote
Old 06-10-2011, 09:16 AM   #5
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Unfortunately, the best documentation is from the original paper (single end) as well as the code (paired end). Try modifying the code to print out the relevant variables to understand the calculation etc.
nilshomer is offline   Reply With Quote
Old 06-20-2011, 11:51 PM   #6
trickytank
Member
 
Location: Melbourne

Join Date: Dec 2010
Posts: 19
Default

Hey I'm interested in this too. In particular, what if one read maps to one location on the reference, but the the other read maps to somewhere differently (such that it does not have the correct orientation and/or distance)? What I really want to know if such pairs are down weighted by low mapping quality in some way?
trickytank is offline   Reply With Quote
Old 06-21-2011, 04:24 AM   #7
dp05yk
Member
 
Location: Brock University

Join Date: Dec 2010
Posts: 66
Default

It says in the paper that BWA will find all single-end alignments for each mate and sort them in ascending order of chromosomal coordinates. Then it uses an estimated insert size to determine which of the chromosomal coordinates are best for both mates.

The insert size is determined in the function infer_isize, and I believe the pairing is determined in the function pairing :-) both are contained in bwape.c.
dp05yk is offline   Reply With Quote
Old 11-14-2011, 04:04 PM   #8
angerusso
Member
 
Location: US

Join Date: Oct 2011
Posts: 47
Default

Hello All,

I have a WholeExome paired end sample and I reached the step where I am performing the alignment to human genome (hg19.fa) on a 10 node cluster.

I am running the command:
bwa aln hg19.fa sample1_1.fastq > sample1_1.sai
bwa aln hg19.fa sample1_2.fastq > sample1_2.sai

But it's taking forever. I understand this could due to couple of reasons, main reason being that I am not doing any pre-filtering. I saw that packages like GenomeQuest do lot of pre-filtering which can make the alignment faster.

I am total new-bie and i am wondering if I can get help here regarding how and what kind of pre-filtering can I run with this sample before using bwa for alignment. I am kind of in a hurry to get some results so any result will be extremely appreciated.

Thanks,
angel
angerusso is offline   Reply With Quote
Old 11-14-2011, 05:28 PM   #9
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

I've run bwa on exome capture DNA with no filtering at all. And takes a while, but it doesn't take forever, and every minute or so it updates the screen telling me how many more reads its finished processing.

Using multiple processors with the -t option would certainly speed things along, if your computer has that capacity.
swbarnes2 is offline   Reply With Quote
Old 11-14-2011, 06:51 PM   #10
angerusso
Member
 
Location: US

Join Date: Oct 2011
Posts: 47
Default

Thanks swbarnes2 very much for your reply.

I hope my files will finish by tomorrow. The size of one paired-end fastq file in my case is 63GB.

I will try the multi-threading mode you mentioned tomorrow.

Angel
angerusso is offline   Reply With Quote
Reply

Tags
bwa, mapping quality, paired end

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 06:30 PM.


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