* samtools-0.1.9-7 (r804)

 * fixed a bug in the gap caller

diff --git a/bam2bcf_indel.c b/bam2bcf_indel.c
index 2254cf038d34d4c7f0938549b85e0c43befe62af..de9df3e6c2fb9042b4287ac60e828224aaa651a2 100644 (file)
--- 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 5395d96bc2a2ed1c21d94257eca91eda62195dcb..e76a9d4499acb9404910732ec0ff624ed11c7b1d 100644 (file)
--- 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[]);