Handle multiple regions at a single locus

166cb68d · Anne-Sophie Denommé-Pichon · 50476c90 · 166cb68d
Commit 166cb68d authored Apr 16, 2022 by Anne-Sophie Denommé-Pichon
Hide whitespace changes
Inline Side-by-side

Showing with 2 additions and 2 deletions

getResults.py getResults.py +2 -2

No files found.
--- a/getResults.py
+++ b/getResults.py
@@ -43,9 +43,9 @@ def enumerate_variants(catalog_path):

        for variant in variants:
            region = variant['ReferenceRegion']
-            # TODO handle other regions?
            if isinstance(region, list):
-                region = region[0]
+                locus_index = variant['VariantId'].index(variant['LocusId'])
+                region = region[locus_index]
            chrom, pos = region.split(':')
            pos = pos.split('-')[0]
            yield variant['LocusId'], (chrom, pos)