From e8a9ef1d4d0518671487cd681aa0e3a59ab60b0a Mon Sep 17 00:00:00 2001 From: Jonas Fuchs <78491186+jonas-fuchs@users.noreply.github.com> Date: Mon, 11 Mar 2024 11:27:48 +0100 Subject: [PATCH] fixed mutations at genome end (#16) --- bamdash/__init__.py | 2 +- bamdash/scripts/data.py | 4 ++-- 2 files changed, 3 insertions(+), 3 deletions(-) diff --git a/bamdash/__init__.py b/bamdash/__init__.py index fbc5807..9bf7ff1 100644 --- 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 index 25cd9aa..0d54fdc 100644 --- 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":