From 0735c4ca30725f8f41070890b0e8eb6439846c54 Mon Sep 17 00:00:00 2001
From: Heng Li <lh3@live.co.uk>
Date: Thu, 18 Nov 2010 02:58:58 +0000
Subject: [PATCH] fixed a bug in printing MNPs

---
 bam2bcf.c | 2 +-
 bamtk.c   | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/bam2bcf.c b/bam2bcf.c
index fe447e6..de414d7 100644
--- a/bam2bcf.c
+++ b/bam2bcf.c
@@ -209,7 +209,7 @@ int bcf_call2bcf(int tid, int pos, bcf_call_t *bc, bcf1_t *b, bcf_callret1_t *bc
 		}
 		kputc('\0', &s);
 	} else if (bc->ori_ref == B2B_REF_MNP) {
-		for (j = 0; j < bca->indelreg; ++j) kputc(ref[pos+1+j], &s);
+		for (j = 0; j < bca->indelreg; ++j) kputc(ref[pos+j], &s);
 		kputc('\0', &s);
 		for (i = 1; i < 4; ++i) {
 			if (bc->a[i] < 0) break;
diff --git a/bamtk.c b/bamtk.c
index 591257b..89f7e80 100644
--- a/bamtk.c
+++ b/bamtk.c
@@ -9,7 +9,7 @@
 #endif
 
 #ifndef PACKAGE_VERSION
-#define PACKAGE_VERSION "0.1.10-2 (r832)"
+#define PACKAGE_VERSION "0.1.10-3 (r833)"
 #endif
 
 int bam_taf2baf(int argc, char *argv[]);
-- 
2.39.2