From af1d6afba9ae1537074ce01cd33e0a9d1b91d8d0 Mon Sep 17 00:00:00 2001 From: Heng Li Date: Thu, 30 Nov 2017 11:34:10 -0500 Subject: [PATCH] r600: bugfix - missing secondary alignments (#71) This should very rarely happen to typical data, but has a higher chance in artifactual data. --- hit.c | 2 +- main.c | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/hit.c b/hit.c index c46c30f..a7c2c8a 100644 --- a/hit.c +++ b/hit.c @@ -234,7 +234,7 @@ void mm_select_sub(void *km, float pri_ratio, int min_diff, int best_n, int *n_, if (p == i || r[i].inv) { // primary or inversion r[k++] = r[i]; } else if ((r[i].score >= r[p].score * pri_ratio || r[i].score + min_diff >= r[p].score) && n_2nd < best_n) { - if (!(r[i].qs == r[p].qs && r[i].qe == r[p].qe && r[i].rs == r[p].rs && r[i].re == r[p].re)) // not identical hits + if (!(r[i].qs == r[p].qs && r[i].qe == r[p].qe && r[i].rid == r[p].rid && r[i].rs == r[p].rs && r[i].re == r[p].re)) // not identical hits r[k++] = r[i], ++n_2nd; else if (r[i].p) free(r[i].p); } else if (r[i].p) free(r[i].p); diff --git a/main.c b/main.c index 74b0ef9..442cba0 100644 --- a/main.c +++ b/main.c @@ -6,7 +6,7 @@ #include "mmpriv.h" #include "getopt.h" -#define MM_VERSION "2.5-r574-dirty" +#define MM_VERSION "2.5-r600-dirty" #ifdef __linux__ #include