fix ambiguous basepair bug.Thanks to Qiaoping Yuan

mehrdadbakhtiari · Mar 24, 2019 · 962232b · 962232b
1 parent 667c2d0
commit 962232b
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/advntr/vntr_finder.py b/advntr/vntr_finder.py
@@ -747,9 +747,9 @@ def get_hash(string):
                 continue
             window_hash = None
             for i in range(0, len(sequence) - keyword_size):
-                if sequence[i].upper() == 'N' or sequence[i - 1 + keyword_size].upper() == 'N':
+                if sequence[i].upper() not in 'ACTG' or sequence[i - 1 + keyword_size].upper() not in 'ACTG':
                     continue
-                if window_hash is None or sequence[i - 1].upper() == 'N':
+                if window_hash is None or sequence[i - 1].upper() not in 'ACTG':
                     if 'N' in sequence[i:i + keyword_size].upper():
                         window_hash = None
                         continue