Question: Are there maximum Penalty Scores in BWA-MEM?
gravatar for Mesohexaploid
6 months ago by
Mesohexaploid 10 wrote:

Apologies if this is a simplistic question, I am new to Galaxy.

I am working on a variant calling workflow for paired DNAseq reads, with BWA-MEM as my aligner. Since I have custom parameters, I would like to manually enter them into BWA-MEM through the "Select full list of options" category.

I was just wondering if there are any maximum values for scoring I should adhere to. For example, if I want an ungapped alignment, what's the highest penalty score I can assign to -O? In the bwa manual, there are values in square brackets behind each argument, but I'm not sure what they represent. The default settings?

alignment bwa-mem bwa • 185 views
ADD COMMENTlink modified 5 months ago • written 6 months ago by Mesohexaploid 10
gravatar for Mesohexaploid
5 months ago by
Mesohexaploid 10 wrote:

Hello again, I'm just going to answer my own question here in case there is anyone else who would like to know the answer to this (doubtful).

I asked this question in the non-Galaxy Biostars, and basically, there are no maximum penalty scores and the scores in the square brackets represent the default values.

The best way to generate an ungapped alignment is to set the penalty value for a gap (-O) really high, like to 100, while keeping all other settings default. This should result in BWA-MEM dropping a seed extension as soon as a gap is encountered.

ADD COMMENTlink written 5 months ago by Mesohexaploid 10
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 82 users visited in the last hour