Fixed a bug in the assignation of CAGE peaks of trasncripts in the ne…

…gative strand
ConesaLab · Jun 13, 2024 · 52def33 · 52def33
1 parent a7c2680
commit 52def33
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/sqanti3_qc.py b/sqanti3_qc.py
@@ -2204,7 +2204,10 @@ def find(self, chrom, strand, query, search_window=10000):
             if strand=='-' and start0<int(query) and end1<int(query):
                 continue
 ##
-            within_out = (start0<=query<end1)
+            if strand == "+":
+                within_out = (start0<=query<end1)
+            if strand == "-":
+                within_out = (start0<query<=end1)
             if within_out:
                 w = 'TRUE'
             else: