* samtools-0.1.9-4 (r801)

 * fixed a minor issue in printing indel VCF

diff --git a/bam2bcf.c b/bam2bcf.c
index c2d994479a9094ea330291db34290e15faca8850..384e481cbc5e6d9d4cae4f2e50a66008f29f1c84 100644 (file)
--- a/bam2bcf.c
+++ b/bam2bcf.c
@@ -181,7 +181,9 @@ int bcf_call2bcf(int tid, int pos, bcf_call_t *bc, bcf1_t *b, bcf_callret1_t *bc
                kputc(ref[pos], &s);
                for (i = 1; i < 4; ++i) {
                        if (bc->a[i] < 0) break;
-                       if (i > 1) kputc(',', &s);
+                       if (i > 1) {
+                               kputc(',', &s); kputc(ref[pos], &s);
+                       }
                        if (bca->indel_types[bc->a[i]] < 0) { // deletion
                                for (j = -bca->indel_types[bc->a[i]]; j < bca->indelreg; ++j)
                                        kputc(ref[pos+1+j], &s);

diff --git a/bamtk.c b/bamtk.c
index 743f33772f9de77edf9f7227a7d25a2fe1ecaf40..35eb2a58aab826c3b59228cc89447be39f0120d0 100644 (file)
--- a/bamtk.c
+++ b/bamtk.c
@@ -9,7 +9,7 @@
 #endif
 
 #ifndef PACKAGE_VERSION
-#define PACKAGE_VERSION "0.1.9-3 (r797)"
+#define PACKAGE_VERSION "0.1.9-4 (r801)"
 #endif
 
 int bam_taf2baf(int argc, char *argv[]);