SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
bwa perfect match only guil Bioinformatics 4 08-30-2013 09:18 AM
samtools pileup output, base quality does not match original ones wwmm933 Bioinformatics 3 11-30-2011 09:52 PM
the output of bwa milkingroom Bioinformatics 1 09-21-2011 01:48 PM
Question on BWA output rboettcher Bioinformatics 3 04-18-2011 03:21 AM
bwa output scami Bioinformatics 1 12-22-2010 02:10 AM

Reply
 
Thread Tools
Old 09-04-2009, 05:28 AM   #1
mingkunli
Member
 
Location: Germany

Join Date: Jan 2009
Posts: 41
Default output all possible match using bwa

I am moving from bowtie to bwa, as bwa considers gap while doing mapping.
But i don't know how to get all the possible hits for a single read on the reference genome(besides the perfect one), like using -a with bowtie. I tried -N, but it does't work. Although X1 in the sam file can give me the number of suboptimal hits, I need the details of the hit like the optimal one.
Anyone can give me some clue, if bwa can't do this, what else can?

ok, I got the result myself, after .sai, using bwa samse -n INT to display maximum n match for each read

Last edited by mingkunli; 09-04-2009 at 07:16 AM.
mingkunli is offline   Reply With Quote
Old 09-04-2009, 08:15 AM   #2
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Quote:
Originally Posted by mingkunli View Post
I am moving from bowtie to bwa, as bwa considers gap while doing mapping.
But i don't know how to get all the possible hits for a single read on the reference genome(besides the perfect one), like using -a with bowtie. I tried -N, but it does't work. Although X1 in the sam file can give me the number of suboptimal hits, I need the details of the hit like the optimal one.
Anyone can give me some clue, if bwa can't do this, what else can?

ok, I got the result myself, after .sai, using bwa samse -n INT to display maximum n match for each read
I think you are looking for the -e option. Set this to the maximum length indel you want detect. My experience shows BWA has difficulty with longer indels (>10bp).

A solution to the indel problem as well as getting "all hits" is to use BFAST, admittedly my own software. I feel like I am doing a lot of self-promotion of late but I guess it is because I am passionate about my solution.

Nils
nilshomer is offline   Reply With Quote
Old 09-05-2009, 08:06 AM   #3
mingkunli
Member
 
Location: Germany

Join Date: Jan 2009
Posts: 41
Default

:)
Thanks for reminding me of the -e option. I will try your BFAST if longer than 10bp indels are expected in my reads.
mingkunli 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 07:23 PM.


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