Check if start and end bounds an exon (#66)

ave-dcd · Oct 16, 2024 · c04c3a0 · c04c3a0
1 parent 88e3078
commit c04c3a0
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/dcd_mapping/lookup.py b/src/dcd_mapping/lookup.py
@@ -186,7 +186,8 @@ async def get_transcripts(
     SELECT tx_ac
     FROM {uta.schema}.tx_exon_aln_v
     WHERE hgnc = '{gene_symbol}'
-      AND ({start} BETWEEN alt_start_i AND alt_end_i OR {end} BETWEEN alt_start_i AND alt_end_i)
+      AND ({start} BETWEEN alt_start_i AND alt_end_i OR {end} BETWEEN alt_start_i AND
+      alt_end_i OR (alt_start_i >= {start} AND alt_end_i <= {end}))
       AND alt_ac = '{chromosome_ac}'
       AND tx_ac NOT LIKE 'NR_%';
     """  # noqa: S608