* samtools-0.1.5-25 (r436)

 * in flagstats, do not report singletons if both ends are unmapped

diff --git a/bam_stat.c b/bam_stat.c
index c1c4a432ce8167448f31577d2ce0f650610ff626..ea9deee9d712c71bfd51f7574a5cae055a2c2ace 100644 (file)
--- a/bam_stat.c
+++ b/bam_stat.c
@@ -16,7 +16,7 @@ typedef struct {
                        if ((c)->flag & BAM_FPROPER_PAIR) ++(s)->n_pair_good;           \
                        if ((c)->flag & BAM_FREAD1) ++(s)->n_read1;                                     \
                        if ((c)->flag & BAM_FREAD2) ++(s)->n_read2;                                     \
-                       if ((c)->flag & BAM_FMUNMAP) ++(s)->n_sgltn;                            \
+                       if (((c)->flag & BAM_FMUNMAP) && !((c)->flag & BAM_FUNMAP)) ++(s)->n_sgltn;     \
                        if (!((c)->flag & BAM_FUNMAP) && !((c)->flag & BAM_FMUNMAP)) { \
                                ++(s)->n_pair_map;                                                                              \
                                if ((c)->mtid != (c)->tid) {                                                    \

diff --git a/bamtk.c b/bamtk.c
index 0b41561eb37100de3244535d41069d9b8f4d3bfc..08d0fdc2e5055bd4dc2f8813d6f93fbe934f2dc7 100644 (file)
--- a/bamtk.c
+++ b/bamtk.c
@@ -9,7 +9,7 @@
 #endif
 
 #ifndef PACKAGE_VERSION
-#define PACKAGE_VERSION "0.1.5-24 (r435)"
+#define PACKAGE_VERSION "0.1.5-25 (r436)"
 #endif
 
 int bam_taf2baf(int argc, char *argv[]);