From 836d46423962b150372f0a95162f7796ba243238 Mon Sep 17 00:00:00 2001
From: Heng Li <lh3@me.com>
Date: Mon, 14 Apr 2014 09:55:55 -0400
Subject: [PATCH] r713: a bug in retrieving ref seq on rev

---
 bntseq.c | 4 ++--
 main.c   | 2 +-
 2 files changed, 3 insertions(+), 3 deletions(-)

diff --git a/bntseq.c b/bntseq.c
index 403d088..eddae84 100644
--- a/bntseq.c
+++ b/bntseq.c
@@ -397,8 +397,8 @@ uint8_t *bns_fetch_seq(const bntseq_t *bns, const uint8_t *pac, int64_t *beg, in
 	far_end = far_beg + bns->anns[*rid].len;
 	if (is_rev) { // flip to the reverse strand
 		int64_t tmp = far_beg;
-		far_beg = (bns->l_pac<<1) - 1 - far_end;
-		far_end = (bns->l_pac<<1) - 1 - tmp;
+		far_beg = (bns->l_pac<<1) - far_end;
+		far_end = (bns->l_pac<<1) - tmp;
 	}
 	*beg = *beg > far_beg? *beg : far_beg;
 	*end = *end < far_end? *end : far_end;
diff --git a/main.c b/main.c
index 07584b9..5f31591 100644
--- a/main.c
+++ b/main.c
@@ -4,7 +4,7 @@
 #include "utils.h"
 
 #ifndef PACKAGE_VERSION
-#define PACKAGE_VERSION "0.7.8-r705-dirty"
+#define PACKAGE_VERSION "0.7.8-r713-dirty"
 #endif
 
 int bwa_fa2pac(int argc, char *argv[]);