Thomas 6 ヶ月 前
コミット
ecb49f0dab
1 ファイル変更7 行追加6 行削除
  1. 7 6
      src/collection/bam.rs

+ 7 - 6
src/collection/bam.rs

@@ -990,12 +990,13 @@ pub enum PileBase {
 
 /// Decode one HTSlib 4-bit nucleotide (0=A,1=C,2=G,3=T,4=N, …) to `PileBase`
 #[inline]
-fn decode(nibble: u8) -> PileBase {
-    match nibble & 0x0f {
-        0 => PileBase::A,
-        1 => PileBase::C,
-        2 => PileBase::G,
-        3 => PileBase::T,
+fn decode(n: u8) -> PileBase {
+    match n & 0x0f {
+        1 => PileBase::A,
+        2 => PileBase::C,
+        4 => PileBase::G,
+        8 => PileBase::T,
+        15 => PileBase::N, // also use N for ambiguous (e.g. 3,5,6,...)
         _ => PileBase::N,
     }
 }