From d9530265899effae829ac332f9da9790f4654f44 Mon Sep 17 00:00:00 2001 From: kirtanav98 <123595850+kirtanav98@users.noreply.github.com> Date: Tue, 2 Jul 2024 13:02:51 -0400 Subject: [PATCH] fixed splitvariants.py to group MEI's with insertions (#687) --- .../04_variant_resolution/scripts/split_variants.py | 4 ++-- 1 file changed, 2 insertions(+), 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 index 1832ddc3a..ec0418459 100644 --- 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