* samtools-0.1.2-7

 * fixed a bug in fixmate: SE reads are flagged as BAM_FMUNMAP

diff --git a/bam_mate.c b/bam_mate.c
index bb5360521adaa1cc8b785b3d7d0f1bbd58bae8fd..f37f8f836b04b499cfeae8dc616e992bfe2c6310 100644 (file)
--- a/bam_mate.c
+++ b/bam_mate.c
@@ -40,8 +40,10 @@ void bam_mating_core(bamFile in, bamFile out)
                                has_prev = 0;
                        } else { // unpaired or singleton
                                pre->core.mtid = -1; pre->core.mpos = -1; pre->core.isize = 0;
-                               pre->core.flag |= BAM_FMUNMAP;
-                               pre->core.flag &= ~BAM_FMREVERSE & ~BAM_FPROPER_PAIR;
+                               if (pre->core.flag & BAM_FPAIRED) {
+                                       pre->core.flag |= BAM_FMUNMAP;
+                                       pre->core.flag &= ~BAM_FMREVERSE & ~BAM_FPROPER_PAIR;
+                               }
                                bam_write1(out, pre);
                        }
                } else has_prev = 1;

diff --git a/bamtk.c b/bamtk.c
index 65b7d69aaa939fe0b8d100f0af771f1f2c21f833..9521fdd2445e01e325fc6aa494e19a340380f377 100644 (file)
--- a/bamtk.c
+++ b/bamtk.c
@@ -3,7 +3,7 @@
 #include "bam.h"
 
 #ifndef PACKAGE_VERSION
-#define PACKAGE_VERSION "0.1.2-7"
+#define PACKAGE_VERSION "0.1.2-8"
 #endif
 
 int bam_taf2baf(int argc, char *argv[]);