From 6de74888fd7b2d922f05a4dcbd0d1dfd125218cc Mon Sep 17 00:00:00 2001
From: Heng Li <lh3@me.com>
Date: Thu, 31 Jan 2013 12:12:58 -0500
Subject: [PATCH] bugfix: min_intv not working in SMEM

---
 bwt.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/bwt.c b/bwt.c
index c4a008b..689d8f8 100644
--- a/bwt.c
+++ b/bwt.c
@@ -323,7 +323,7 @@ int bwt_smem1(const bwt_t *bwt, int len, const uint8_t *q, int x, int min_intv,
 					}
 				} // otherwise the match is contained in another longer match
 			}
-			if (c >= 0 && (ok[c].x[2] && (curr->n == 0 || ok[c].x[2] != curr->a[curr->n-1].x[2]))) {
+			if (c >= 0 && ok[c].x[2] >= min_intv && (curr->n == 0 || ok[c].x[2] != curr->a[curr->n-1].x[2])) {
 				ok[c].info = p->info;
 				kv_push(bwtintv_t, *curr, ok[c]);
 			}