zzh
/
fast-bwa
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

bugfix: a typo leading to wrong pairing

This commit is contained in:
Heng Li 2013-02-22 10:57:07 -05:00
parent 58e4cc207f
commit 81fe6f8e38
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
bwamem_pair.c
View File

@ -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].high) break;
if (dist < pes[dir].low) continue; if (dist < pes[dir].low) continue;
ns = (dist - pes[dir].avg) / pes[dir].std; ns = (dist - pes[dir].avg) / pes[dir].std;
q = (int)((v.a[i].y>>32) + (v.a[i].y>>32) - 4.343 * log(erfc(fabs(ns) * M_SQRT1_2)) + .499); q = (int)((v.a[i].y>>32) + (v.a[k].y>>32) - 4.343 * log(erfc(fabs(ns) * M_SQRT1_2)) + .499);
pair = (uint64_t)k<<32 | i; pair = (uint64_t)k<<32 | i;
x = (uint64_t)q<<32 | (hash_64(pair ^ id<<8) & 0xffffffffU); x = (uint64_t)q<<32 | (hash_64(pair ^ id<<8) & 0xffffffffU);
if (x < o.x) subo = o, o.x = x, o.y = pair; if (x < o.x) subo = o, o.x = x, o.y = pair;