fixed splitvariants.py to group MEI's with insertions (#687)

broadinstitute · Jul 2, 2024 · d953026 · d953026
1 parent b80073d
commit d953026
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/sv-pipeline/04_variant_resolution/scripts/split_variants.py b/src/sv-pipeline/04_variant_resolution/scripts/split_variants.py
@@ -12,8 +12,8 @@ def process_bed_file(input_bed, n_per_split, bca=True):
     condition_prefixes = {
         'gt5kb': {'condition': lambda line: (line[SVTYPE_FIELD] == 'DEL' or line[SVTYPE_FIELD] == 'DUP') and (int(line[END_FIELD]) - int(line[START_FIELD]) >= 5000)},
         'lt5kb': {'condition': lambda line: (line[SVTYPE_FIELD] == 'DEL' or line[SVTYPE_FIELD] == 'DUP') and (int(line[END_FIELD]) - int(line[START_FIELD]) < 5000)},
-        'bca': {'condition': lambda line: bca and line[SVTYPE_FIELD] not in ['DEL', 'DUP', 'INS']},
-        'ins': {'condition': lambda line: bca and line[SVTYPE_FIELD] == 'INS'}
+        'bca': {'condition': lambda line: bca and line[SVTYPE_FIELD] not in ['DEL', 'DUP'] and not line[SVTYPE_FIELD].startswith('INS')},
+        'ins': {'condition': lambda line: bca and line[SVTYPE_FIELD].startswith('INS')}
     }
 
     # Create trackers for the current file information