bugfix: long-existing out-of-boundary bug
This commit is contained in:
parent
91a4a0c8ea
commit
a471f1918b
2
bwase.c
2
bwase.c
|
|
@ -245,7 +245,7 @@ char *bwa_cal_md1(int n_cigar, bwa_cigar_t *cigar, int len, bwtint_t pos, ubyte_
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
} else { // no gaps
|
} else { // no gaps
|
||||||
for (z = u = 0; z < (bwtint_t)len; ++z) {
|
for (z = u = 0; z < (bwtint_t)len && x+z < l_pac; ++z) {
|
||||||
c = pacseq[(x+z)>>2] >> ((~(x+z)&3)<<1) & 3;
|
c = pacseq[(x+z)>>2] >> ((~(x+z)&3)<<1) & 3;
|
||||||
if (c > 3 || seq[y+z] > 3 || c != seq[y+z]) {
|
if (c > 3 || seq[y+z] > 3 || c != seq[y+z]) {
|
||||||
ksprintf(str, "%d", u);
|
ksprintf(str, "%d", u);
|
||||||
|
|
|
||||||
Loading…
Reference in New Issue