From f4aedddee6d55dffad46060f5bfd8099f8d22024 Mon Sep 17 00:00:00 2001
From: Heng Li <lh3@me.com>
Date: Mon, 8 Sep 2014 11:32:48 -0400
Subject: [PATCH] r819: bugfix - added too many sub-SMEMs

---
 bwamem.c | 2 +-
 main.c   | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/bwamem.c b/bwamem.c
index 8377122..4fd6b28 100644
--- a/bwamem.c
+++ b/bwamem.c
@@ -137,7 +137,7 @@ static void mem_collect_intv(const mem_opt_t *opt, const bwt_t *bwt, int len, co
 		if (end - start < split_len || p->x[2] > opt->split_width) continue;
 		bwt_smem1(bwt, len, seq, (start + end)>>1, p->x[2]+1, &a->mem1, a->tmpv);
 		for (i = 0; i < a->mem1.n; ++i)
-			if ((a->mem1.a[i].info>>32) - (uint32_t)a->mem1.a[i].info >= opt->min_seed_len)
+			if ((uint32_t)a->mem1.a[i].info - (a->mem1.a[i].info>>32) >= opt->min_seed_len)
 				kv_push(bwtintv_t, a->mem, a->mem1.a[i]);
 	}
 	// sort
diff --git a/main.c b/main.c
index f5e65c1..93bd84e 100644
--- a/main.c
+++ b/main.c
@@ -4,7 +4,7 @@
 #include "utils.h"
 
 #ifndef PACKAGE_VERSION
-#define PACKAGE_VERSION "0.7.10-r815-dirty"
+#define PACKAGE_VERSION "0.7.10-r819-dirty"
 #endif
 
 int bwa_fa2pac(int argc, char *argv[]);