From: Heng Li <lh3@live.co.uk>
Date: Tue, 9 Nov 2010 19:12:06 +0000 (+0000)
Subject:  * samtools-0.1.9-7 (r804)
X-Git-Url: https://git.donarmstrong.com/?a=commitdiff_plain;ds=sidebyside;h=ed191cc20851f34d1bf87db5eba6da6d6c9cbd92;p=samtools.git

 * samtools-0.1.9-7 (r804)
 * fixed a bug in the gap caller
---

diff --git a/bam2bcf_indel.c b/bam2bcf_indel.c
index 2254cf0..de9df3e 100644
--- a/bam2bcf_indel.c
+++ b/bam2bcf_indel.c
@@ -66,12 +66,12 @@ int bcf_call_gap_prep(int n, int *n_plp, bam_pileup1_t **plp, int pos, bcf_calla
 	if (ref == 0 || bca == 0) return -1;
 	// determine if there is a gap
 	for (s = N = 0; s < n; ++s) {
-		N += n_plp[s]; // N is the total number of reads
 		for (i = 0; i < n_plp[s]; ++i)
 			if (plp[s][i].indel != 0) break;
 		if (i < n_plp[s]) break;
 	}
 	if (s == n) return -1; // there is no indel at this position.
+	for (s = N = 0; s < n; ++s) N += n_plp[s]; // N is the total number of reads
 	{ // find out how many types of indels are present
 		int m;
 		uint32_t *aux;
diff --git a/bamtk.c b/bamtk.c
index 5395d96..e76a9d4 100644
--- a/bamtk.c
+++ b/bamtk.c
@@ -9,7 +9,7 @@
 #endif
 
 #ifndef PACKAGE_VERSION
-#define PACKAGE_VERSION "0.1.9-6 (r803)"
+#define PACKAGE_VERSION "0.1.9-7 (r804)"
 #endif
 
 int bam_taf2baf(int argc, char *argv[]);