What Is The Recommended Way To Use Bwa On Files Containing Haplotypes Of An Exome
2
0
Entering edit mode
11.2 years ago
agatorano ▴ 50

Hello

I am unsure how best to use BWA on two files containing haplotypes of an exome. There are three BWA types. I assume I am to be using the MEM algorithm but I am unsure what many of my options should be. For example what should my gap costs be? match? Are these things well established? what are rules to use when trying to decide on the parameters for this algorithm?

Thank you for your help.

bwa exome • 4.5k views
ADD COMMENT
3
Entering edit mode
11.2 years ago
Leandro Lima ▴ 970

Hello!

Maybe this post can help you.

What is the best pipeline for human whole exome sequencing?

ADD COMMENT
0
Entering edit mode

This is really helpful in the long term but even they admit to not knowing the parameters to use and why.

ADD REPLY
0
Entering edit mode

Also, this post is from time before BWA mem algoritm so it does not answer to the question of OP in that way.

ADD REPLY
2
Entering edit mode
11.2 years ago
onter ▴ 170

Typically the default parameters have been tested and optimized on the data they are prescribed to. For example, BWA-MEM is excellent for 75+bp reads. There are also tools online for determining which kind pipeline is suitable for your sample - one of these is GCAT on bioplanet.com which allows you to view benchmarks on pipelines as well as run your own. Below is a GCAT summary of a BWA pipeline comparison with GATK, Freebayes, and ATLAS on a 30x illumina exome:

View Full GCAT Variant Calls Report: Bwa+Atlas (illumina-100bp-pe-exome-30x)

If you are uncertain about what pipeline or parameters to use, I really suggest you give GCAT a shot

ADD COMMENT
0
Entering edit mode

This is extremely interesting

ADD REPLY

Login before adding your answer.

Traffic: 3205 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6