SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
bwa samse -s parameter rahul Bioinformatics 6 03-01-2012 02:44 PM
bwa insert size estimation athena.uci Bioinformatics 2 11-07-2011 08:49 AM
a confusing parameter of BWA chenjy Bioinformatics 0 12-08-2010 12:15 AM
BWA parameter for microRNA m_elena_bioinfo Bioinformatics 0 11-19-2010 07:11 AM
population genetics parameter estimation of solexa population sequencing baohua100 Bioinformatics 1 07-16-2008 09:53 AM

Reply
 
Thread Tools
Old 03-11-2013, 04:36 AM   #1
opulcy97
Member
 
Location: moon

Join Date: Nov 2012
Posts: 19
Default BWA automatic parameter estimation

Hi,

I am using Illumina paired-end data (read length 250 bp) for producing a bam file using BWA and Samtools. As BWA does not support read length greater than 200 bp. I had to use BWA-MEM in this case with some testing parameters (such as gap opening penalty 30, gap extension penalty 5, etc.). However, one of my colleagues used BWA from Miseq machine and got a totally different alignment and mapping.

As I am trying to replicate her work process, I thought of using the same parameters as she used. However, when I went through Miseq manual, it appears that Miseq-BWA automatically adjusts parameters based on read lengths and error rates, and then estimates insert size distribution (MiSeq manual page 22: BWA). From the bam file, from Miseq-BWA, we found that the genome has some deletions and this is actually authenticated by Sanger sequencing. So in this case, I am more or less sure that MiSeq-BWA is doing the right thing. However, the bam file created by my pipeline using BWA-MEM does not show the same result.

Now my questions are:

1. Is Miseq-BWA and the available BWA tools differently implemented? If not, then which one I should be using, BWA-SW or BWA-MEM?

2. Is there any way of adjusting parameters based on read lengths and error rates, and then estimating insert size distribution? If I do not specify any parameter, BWA should be taking the default one, I guess.

I know the questions are quite broad. But it would be great help if anyone has some advice for me.

Last edited by opulcy97; 03-11-2013 at 04:36 AM. Reason: Typo
opulcy97 is offline   Reply With Quote
Old 03-11-2013, 10:36 AM   #2
Bukowski
Senior Member
 
Location: UK

Join Date: Jan 2010
Posts: 390
Default

bwa-mem was introduced in 0.7.0 - I'm pretty sure the version with the MiSeq software is a 0.6 derivative. I'm not sure that helps but you're not doing a like for like comparison.
Bukowski is offline   Reply With Quote
Reply

Tags
alignment, bwa, bwa-mem, miseq

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 05:56 AM.


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