fixed mutations at genome end (#16)

jonas-fuchs · Mar 11, 2024 · e8a9ef1 · e8a9ef1
1 parent 09bf627
commit e8a9ef1
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/bamdash/__init__.py b/bamdash/__init__.py
@@ -1,3 +1,3 @@
 """interactively visualize coverage and tracks"""
 _program = "bamdash"
-__version__ = "0.2.3"
+__version__ = "0.2.4"
diff --git a/bamdash/scripts/data.py b/bamdash/scripts/data.py
@@ -418,8 +418,8 @@ def analyse_a3_signature(variant, seq):
                 return "YES", f"{site[0]}>{site[1]}<"
             else:
                 return "NO", f"{site[0]}>{site[1]}<"
-        # for G ref check the following nt
-        if variant["reference"] == "G":
+        # for G ref check the following nt if the position is not at the end of the sequence
+        if variant["reference"] == "G" and variant["position"] != len(seq):
             site = str(seq[pos-1:pos+1])
             # check if the site and mutation type is in line with A3A activity
             if site == "GA" and variant["mutation"] == "A":