[critica] protein fix 5' to 3' reverse complement

Chokyotager · Sep 12, 2021 · 96aaf28 · 96aaf28
1 parent 4ccd90e
commit 96aaf28
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/orffinder/orffinder.py b/orffinder/orffinder.py
@@ -215,7 +215,7 @@ def getORFProteins (sequence, translation_table=1, return_loci=False, **kwargs):
             locus["protein"] = Seq(forward[locus["start"] - 1 : locus["end"] - 1 - difference]).translate()
 
         else:
-            locus["protein"] = Seq(reverse[sequence_length - locus["start"] - 1 : sequence_length - locus["end"] - 1  - difference]).translate()
+            locus["protein"] = Seq(reverse[sequence_length - locus["start"] + 1 : sequence_length - locus["end"] + 1  - difference]).translate()
 
         proteins.append(locus["protein"])