From 205eadf3467170d79563db324dbaf9349a615012 Mon Sep 17 00:00:00 2001 From: Heng Li Date: Sat, 23 Jan 2021 16:25:52 -0500 Subject: [PATCH] r334: ignore a high-occ seed in a small window --- CommandLines.h | 2 +- sketch.cpp | 3 ++- 2 files changed, 3 insertions(+), 2 deletions(-) diff --git a/CommandLines.h b/CommandLines.h index 56f23d5..d7e753b 100644 --- a/CommandLines.h +++ b/CommandLines.h @@ -3,7 +3,7 @@ #include -#define HA_VERSION "0.13-r333-dirty" +#define HA_VERSION "0.13-r334-dirty" #define VERBOSE 0 diff --git a/sketch.cpp b/sketch.cpp index f053073..72f9d39 100644 --- a/sketch.cpp +++ b/sketch.cpp @@ -67,7 +67,8 @@ static void select_mz(ha_mz1_v *p, int len, int sample_dist) } //ks_heapsort_mz(k, b); // sorting is not needed for now for (j = 0; j < k; ++j) - p->a[b[j].pos].rid = 0; + if (b[j].rid < pe - ps) + p->a[b[j].pos].rid = 0; } last0 = i; }