the previous change does not work... Fixed.

This commit is contained in:
Heng Li 2013-02-08 14:51:24 -05:00
parent fdb0a7405f
commit 220fc39e9d
1 changed files with 1 additions and 1 deletions

View File

@ -402,8 +402,8 @@ void mem_chain2aln(const mem_opt_t *opt, int64_t l_pac, const uint8_t *pac, int
if (t->rbeg + t->len > a->re || t->qbeg + t->len > a->qe) if (t->rbeg + t->len > a->re || t->qbeg + t->len > a->qe)
break; break;
} }
if (a->score >= best.score) csub = best.score, best = *a;
if (i >= c->n) break; // all seeds are included; no need to proceed if (i >= c->n) break; // all seeds are included; no need to proceed
if (a->score > best.score) csub = best.score, best = *a;
k = i; k = i;
} }
if (a->score < best.score) *a = best; if (a->score < best.score) *a = best;