r295: gap closure from the middle of non-HPC k
This WILL slightly affect the result of genomic mapping, but hopefully in the good direction.
This commit is contained in:
parent
43506edbc5
commit
28f86688ab
4
align.c
4
align.c
|
|
@ -164,8 +164,8 @@ static inline void mm_adjust_minier(const mm_idx_t *mi, uint8_t *const qseq0[2],
|
||||||
c = mm_get_hplen_back(mi, a->x<<1>>33, (int32_t)a->x);
|
c = mm_get_hplen_back(mi, a->x<<1>>33, (int32_t)a->x);
|
||||||
*r = (int32_t)a->x + 1 - c;
|
*r = (int32_t)a->x + 1 - c;
|
||||||
} else {
|
} else {
|
||||||
*r = (int32_t)a->x + 1;
|
*r = (int32_t)a->x - (mi->k>>1);
|
||||||
*q = (int32_t)a->y + 1;
|
*q = (int32_t)a->y - (mi->k>>1);
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|
|
||||||
Loading…
Reference in New Issue