chore[call]: include ref snv bases in snv dict too

davidlougheed · Apr 18, 2023 · 0121ef0 · 0121ef0
1 parent cf09e80
commit 0121ef0
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/strkit/call/caller/snvs.py b/strkit/call/caller/snvs.py
@@ -173,7 +173,7 @@ def get_read_snvs_dbsnp(
             continue
         if ref is not None and len(ref) == 1 and alts and all(len(a) == 1 for a in alts):
             read_base = _find_base_at_pos(query_sequence, pairs, pos)
-            if read_base in alts:
+            if read_base == ref or read_base in alts:
                 snvs[pos] = read_base
 
     return snvs