zzh
/
minimap2
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

backup

This commit is contained in:
Heng Li 2017-04-09 15:09:39 -04:00
parent 79c9478f46
commit 1bf3ae6752
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
map.c
View File

@ -49,8 +49,8 @@ const mm_reg1_t *mm_map(const mm_idx_t *mi, int l_seq, const char *seq, int *n_r
const uint64_t *dreg;
dreg = sdust_core((const uint8_t*)seq, l_seq, opt->sdust_thres, 64, &n_dreg, b->sdb);
for (j = k = 0; j < b->mini.n; ++j) {
int32_t qpos = (uint32_t)b->mini.a[j].y>>1;
int32_t s = qpos - (mi->k - 1), e = s + mi->k;
int32_t qpos = (uint32_t)b->mini.a[j].y>>1, span = b->mini.a[j].x&0xff;
int32_t s = qpos - (span - 1), e = s + span;
while (u < n_dreg && (uint32_t)dreg[u] <= s) ++u;
if (u < n_dreg && dreg[u]>>32 < e) {
int v, l = 0;