Question: How To Adjust Bwa Mem To Allow For Longer Gaps In Contig Alignment
gravatar for William
7.3 years ago by
William4.7k wrote:

My first guess would be to adjust the gap extension parameter. But that one is already set to 1, the lowest integer. I now get multiple split alignment, that are relatively very close to each other and I would rather have more of a end to end alignment.

Algorithm options:

   -t INT     number of threads [1]
   -k INT     minimum seed length [19]
   -w INT     band width for banded alignment [100]
   -d INT     off-diagonal X-dropoff [100]
   -r FLOAT   look for internal seeds inside a seed longer than {-k} * FLOAT [1.5]
   -c INT     skip seeds with more than INT occurrences [10000]
   -S         skip mate rescue
   -P         skip pairing; mate rescue performed unless -S also in use
   -A INT     score for a sequence match [1]
   -B INT     penalty for a mismatch [4]
   -O INT     gap open penalty [6]
   -E INT     gap extension penalty; a gap of size k cost {-O} + {-E}*k [1]
   -L INT     penalty for clipping [5]
   -U INT     penalty for an unpaired read pair [17]
ADD COMMENTlink modified 7 weeks ago by wljlinksly50 • written 7.3 years ago by William4.7k
gravatar for thackl
5.1 years ago by
thackl2.8k wrote:

 If you want to follow through on your idea about decreasing gap costs, you can do something like -A 2 -E 1. At least in recent versions of bwa -A .. which scales options -TdBOELU unless overridden scales all costs, while -E will overwrite the scaled value, essentially decreasing extension costs.

ADD COMMENTlink modified 5.1 years ago • written 5.1 years ago by thackl2.8k
gravatar for wljlinksly
7 weeks ago by
wljlinksly50 wrote:

Try to set -O to 20 and -E to 0

ADD COMMENTlink written 7 weeks ago by wljlinksly50
Please log in to add an answer.


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