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

no effective changes

This commit is contained in:
Heng Li 2013-02-01 12:57:48 -05:00
parent 5d372cef65
commit 86f2e134ba
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
bwt.c
View File

@ -296,9 +296,10 @@ int bwt_smem1(const bwt_t *bwt, int len, const uint8_t *q, int x, int min_intv,
if (q[i] < 4) { // an A/C/G/T base if (q[i] < 4) { // an A/C/G/T base
c = 3 - q[i]; // complement of q[i] c = 3 - q[i]; // complement of q[i]
bwt_extend(bwt, &ik, ok, 0); bwt_extend(bwt, &ik, ok, 0);
if (ok[c].x[2] != ik.x[2]) // change of the interval size if (ok[c].x[2] != ik.x[2]) { // change of the interval size
kv_push(bwtintv_t, *curr, ik); kv_push(bwtintv_t, *curr, ik);
if (ok[c].x[2] < min_intv) break; // the interval size is too small to be extended further if (ok[c].x[2] < min_intv) break; // the interval size is too small to be extended further
}
ik = ok[c]; ik.info = i + 1; ik = ok[c]; ik.info = i + 1;
} else { // an ambiguous base } else { // an ambiguous base
kv_push(bwtintv_t, *curr, ik); kv_push(bwtintv_t, *curr, ik);