fork/minimap2
1
0
Fork 0
mirror of https://github.com/lh3/minimap2.git synced 2026-10-04 23:38:12 +08:00
Code Issues Packages Projects Releases Wiki Activity

r105: fixed a bug in repeated right ext when zdrop

Browse Source
This commit is contained in:
Heng Li
2017-06-29 15:45:07 -04:00
parent b9075d39a8
commit 337c2a21cd
2 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files
align.c
+2 -2
View File
@@ -240,7 +240,7 @@ static void mm_align1(void *km, const mm_mapopt_t *opt, const mm_idx_t *mi, int
mm_append_cigar(r, ez->n_cigar, ez->cigar); mm_append_cigar(r, ez->n_cigar, ez->cigar);
mm_update_extra(r->p, qseq, tseq, ez->n_cigar, ez->cigar, 0); mm_update_extra(r->p, qseq, tseq, ez->n_cigar, ez->cigar, 0);
} }
if (ez->zdropped) { // truncated by Z-drop if (ez->zdropped) { // truncated by Z-drop; TODO: sometimes Z-drop kicks in because the next seed placement is wrong. This can be fixed in principle.
int j; int j;
for (j = i - 1; j >= 0; --j) for (j = i - 1; j >= 0; --j)
if ((int32_t)a[r->as + j].x < re + ez->max_t) if ((int32_t)a[r->as + j].x < re + ez->max_t)
@@ -261,7 +261,7 @@ static void mm_align1(void *km, const mm_mapopt_t *opt, const mm_idx_t *mi, int
} }
} }
if (i == r->cnt && qe < qe0 && re < re0) { // right extension if (!r->split && qe < qe0 && re < re0) { // right extension
qseq = &qseq0[rev][qe]; qseq = &qseq0[rev][qe];
mm_idx_getseq(mi, rid, re, re0, tseq); mm_idx_getseq(mi, rid, re, re0, tseq);
ksw_extz2_sse(km, qe0 - qe, qseq, re0 - re, tseq, 5, mat, opt->q, opt->e, bw, opt->zdrop, KSW_EZ_EXTZ_ONLY, ez); ksw_extz2_sse(km, qe0 - qe, qseq, re0 - re, tseq, 5, mat, opt->q, opt->e, bw, opt->zdrop, KSW_EZ_EXTZ_ONLY, ez);
main.c
+1 -1
View File
@@ -10,7 +10,7 @@
#include "minimap.h" #include "minimap.h"
#include "mmpriv.h" #include "mmpriv.h"
#define MM_VERSION "2.0-r104-pre" #define MM_VERSION "2.0-r105-pre"
void liftrlimit() void liftrlimit()
{ {