ignore ref bases other than A/T/C/G

OpenGene · Aug 20, 2018 · 592c8e0 · 592c8e0
1 parent 6428fd4
commit 592c8e0
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 1 deletion.
diff --git a/src/cluster.cpp b/src/cluster.cpp
@@ -465,6 +465,9 @@ int Cluster::makeConsensus(vector<bam1_t* >& reads, bam1_t* out, vector<char*>&
                 refbase = refdata[refpos];
         }
 
+        if(refbase!='A' && refbase!='T' && refbase!='C' && refbase!='G')
+            refbase = 0;
+
         // the secondary base is a single base
         if(secNum ==1){
             // low quality secondary

diff --git a/src/common.h b/src/common.h
@@ -1,7 +1,7 @@
 #ifndef COMMON_H
 #define COMMON_H
 
-#define VERSION_NUMBER "0.5.0"
+#define VERSION_NUMBER "0.6.0"
 
 #define _DEBUG false