From ccab2888c70c04a45138e7b853d4fb0ca4062fe7 Mon Sep 17 00:00:00 2001
From: Martijn Vermaat <martijn@vermaat.name>
Date: Thu, 29 Mar 2012 13:23:58 +0000
Subject: [PATCH] Fix crossmapping on reverse strand if CDS start is exon start
This bug was noticed and reported by S Venkata Suresh Kumar.
git-svn-id: https://humgenprojects.lumc.nl/svn/mutalyzer/trunk@503 eb6bd6ab-9ccd-42b9-aceb-e2899b4a52f1
---
mutalyzer/Crossmap.py | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/mutalyzer/Crossmap.py b/mutalyzer/Crossmap.py
index 9354927b..d806e854 100644
--- a/mutalyzer/Crossmap.py
+++ b/mutalyzer/Crossmap.py
@@ -226,7 +226,7 @@ class Crossmap() :
# CDS start.
i = y - c
# Find CDS start.
- while d * (self.RNA[i] - ((i + 1) % 2)) < d * self.CDS[c] :
+ while d * (self.RNA[i] - ((i + 1) % 2)) < d * self.CDS[c] + c :
i += d
cPos = d * (self.RNA[i] - self.CDS[c] + (d * 2)) # Get the right
# boundary.
--
GitLab