entry.within edge case

Variants on chromosomes not in a tree return false
ACEnglish · Jan 11, 2025 · 624d72a · 624d72a
1 parent f90f194
commit 624d72a
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/truvari/variant_record.py b/truvari/variant_record.py
@@ -960,6 +960,8 @@ def within_tree(self, tree):
         Extract entry and tree boundaries to call `truvari.coords_within`
         """
         qstart, qend = self.boundaries()
+        if self.chrom not in tree:
+            return False
         m_ovl = tree[self.chrom].overlap(qstart, qend)
         if len(m_ovl) != 1:
             return False