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

bugfix: min_intv not working in SMEM

This commit is contained in:
Heng Li 2013-01-31 12:12:58 -05:00
parent 543c719a54
commit 6de74888fd
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
bwt.c
View File

@ -323,7 +323,7 @@ int bwt_smem1(const bwt_t *bwt, int len, const uint8_t *q, int x, int min_intv,
}
} // otherwise the match is contained in another longer match
}
if (c >= 0 && (ok[c].x[2] && (curr->n == 0 || ok[c].x[2] != curr->a[curr->n-1].x[2]))) {
if (c >= 0 && ok[c].x[2] >= min_intv && (curr->n == 0 || ok[c].x[2] != curr->a[curr->n-1].x[2])) {
ok[c].info = p->info;
kv_push(bwtintv_t, *curr, ok[c]);
}