SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Annotate contigs with BLAST hit names; remove contigs with no hit Bueller_007 Bioinformatics 10 02-27-2013 10:22 AM
how to set up max mismatches be 15 via bwa shuang Bioinformatics 1 02-02-2012 03:17 PM
how to use bwa? and how to set paremeters for 75bp reads Tlexander Bioinformatics 3 04-17-2011 02:35 AM
What is reciprocal best hit? Dinesh Bioinformatics 1 02-08-2011 03:51 AM
How can i set the "mismatch" parameter in bwa genelab Bioinformatics 14 01-12-2011 02:08 PM

Reply
 
Thread Tools
Old 08-07-2010, 08:59 PM   #1
genelab
Member
 
Location: honston

Join Date: Nov 2009
Posts: 27
Default How to set bwa to resport all the best repeat hit?

Hello all,
How to set bwa to resport all the best repeat hits? The default is only to report the random choosed one from all the best repeat hits.
All reply will be great appreciated. thx
genelab is offline   Reply With Quote
Old 08-07-2010, 09:49 PM   #2
john_mu
Member
 
Location: Stanford, CA

Join Date: May 2010
Posts: 88
Default

I thought it did this by default?

from the manual:

By default, if the best hit is no so repetitive (controlled by -R), BWA also finds all hits contains one more mismatch; otherwise, BWA finds all equally best hits only. Base quality is NOT considered in evaluating hits. In paired-end alignment, BWA pairs all hits it found. It further performs Smith-Waterman alignment for unmapped reads with mates mapped to rescue mapped mates, and for high-quality anomalous pairs to fix potential alignment errors.
__________________
SpliceMap: De novo detection of splice junctions from RNA-seq
Download SpliceMap Comment here
john_mu is offline   Reply With Quote
Old 08-08-2010, 12:14 AM   #3
genelab
Member
 
Location: honston

Join Date: Nov 2009
Posts: 27
Default

Quote:
Originally Posted by john_mu View Post
I thought it did this by default?

from the manual:

By default, if the best hit is no so repetitive (controlled by -R), BWA also finds all hits contains one more mismatch; otherwise, BWA finds all equally best hits only. Base quality is NOT considered in evaluating hits. In paired-end alignment, BWA pairs all hits it found. It further performs Smith-Waterman alignment for unmapped reads with mates mapped to rescue mapped mates, and for high-quality anomalous pairs to fix potential alignment errors.
Thanks for your kindly reply.
This may be the default setting of "bwa aln", however, when i use 'bwa samse' to convert the result.sai ,generated by 'bwa aln', to sam format, the read in the sam result hava at most one hit.

form the samse part of manual:
Samse generate alignments in the SAM format given single-end reads. Repetitive hits will be randomly chosen
genelab is offline   Reply With Quote
Old 08-08-2010, 09:51 AM   #4
john_mu
Member
 
Location: Stanford, CA

Join Date: May 2010
Posts: 88
Default

If I remember correctly, it outputs the other ones in the "XA" tag at the end of the line.

if you set the -n option to a high number, you should get most of the hits.

I think the hits are chosen randomly in this case because all the repetitive hits are equally good.
__________________
SpliceMap: De novo detection of splice junctions from RNA-seq
Download SpliceMap Comment here
john_mu is offline   Reply With Quote
Reply

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 11:51 AM.


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