Skip to content

HTTPS clone URL

Subversion checkout URL

You can clone with
or
.
Download ZIP
Browse files

r364: bug in mem pairing (no effect with -A=1)

Forgot to adjust for matching score. This bug has no effect when -A takes the
default value.
  • Loading branch information...
commit e5355fe3a055307aaf8f6868770573a6c6aef0b4 1 parent bdf34f6
@lh3 authored
Showing with 2 additions and 2 deletions.
  1. +1 −1  bwamem_pair.c
  2. +1 −1  main.c
View
2  bwamem_pair.c
@@ -201,7 +201,7 @@ int mem_pair(const mem_opt_t *opt, int64_t l_pac, const uint8_t *pac, const mem_
if (dist > pes[dir].high) break;
if (dist < pes[dir].low) continue;
ns = (dist - pes[dir].avg) / pes[dir].std;
- q = (int)((v.a[i].y>>32) + (v.a[k].y>>32) + .721 * log(2. * erfc(fabs(ns) * M_SQRT1_2)) + .499); // .721 = 1/log(4)
+ q = (int)((v.a[i].y>>32) + (v.a[k].y>>32) + .721 * log(2. * erfc(fabs(ns) * M_SQRT1_2)) * opt->a + .499); // .721 = 1/log(4)
if (q < 0) q = 0;
p = kv_pushp(pair64_t, u);
p->y = (uint64_t)k<<32 | i;
View
2  main.c
@@ -3,7 +3,7 @@
#include "utils.h"
#ifndef PACKAGE_VERSION
-#define PACKAGE_VERSION "0.7.2-r363-beta"
+#define PACKAGE_VERSION "0.7.2-r364-beta"
#endif
int bwa_fa2pac(int argc, char *argv[]);
Please sign in to comment.
Something went wrong with that request. Please try again.