SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
BWA MEM mate pair incorrect insert sizes? hartmaier Bioinformatics 7 07-24-2014 11:38 AM
BWA-MEM: output mapped reads larger than input reads Kennels Bioinformatics 9 09-05-2013 05:27 PM
BWA rescue of multi-mapping or unmapped reads kjlee Bioinformatics 6 07-13-2012 10:29 PM
Paired-end reads mapped 50% while each pair reads mapped 80% zack80.liu Bioinformatics 3 03-03-2011 01:06 AM
bwa for mate pair reads talk24 Bioinformatics 1 03-29-2010 07:37 PM

Reply
 
Thread Tools
Old 09-05-2014, 06:06 AM   #1
Flo1000
Member
 
Location: Germany

Join Date: Sep 2014
Posts: 12
Default BWA MEM - which reads were mapped by mate-pair rescue?

Is there a way to tell which reads were mapped by mate rescue?
I learned that in BWA results there should be a XT:A:M tag in that case, but in my sam files I can't find any XT tags at all...
Flo1000 is offline   Reply With Quote
Old 09-08-2014, 12:00 AM   #2
Flo1000
Member
 
Location: Germany

Join Date: Sep 2014
Posts: 12
Default

For those who are interested: Apparently, there is no way to tell if a read was mapped by mate pair rescue. From the BWA mailing list:

Quote:
For bwa-backtrack, rescue is performed after mapping, but for bwa-mem, rescue is tightly integrated with mapping. The hits resulted from rescue and from single-end mapping are mostly indistinguishable.
Flo1000 is offline   Reply With Quote
Old 09-08-2014, 08:05 AM   #3
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

What do you plan to use that information for, anyway?
Brian Bushnell is offline   Reply With Quote
Old 09-08-2014, 10:57 PM   #4
Flo1000
Member
 
Location: Germany

Join Date: Sep 2014
Posts: 12
Default

The authors of this paper suggest that "Simply removing SNVs for which 7% or more of the mapped reads were mate-rescued increased the validation rate [..] without eliminating any TPs". I wanted to give it a try within our pipeline.
Flo1000 is offline   Reply With Quote
Old 09-09-2014, 08:01 AM   #5
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

Assuming that paper was based on bwa-aln rather than bwa-mem, its results should have little bearing on bwa-mem, since they do rescue differently. But you can avoid rescue altogether simply by mapping all the reads as single-ended rather than paired to see if you get better results.
Brian Bushnell is offline   Reply With Quote
Reply

Tags
bwa, bwa mem, mate pair rescue, mate rescue

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 04:12 AM.


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