#355 fix chromosome comparison check

bcgsc · Aug 16, 2023 · cf2395f · cf2395f
1 parent 6446465
commit cf2395f
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/mavis/validate/align.py b/src/mavis/validate/align.py
@@ -563,7 +563,7 @@ def supports_primary_event(alignment):
         std_reads = set()
         alignments = []
         for raw_read in reads_by_query.get(contig.seq, []):
-            if raw_read.reference_name not in {evidence.break1.chr, evidence.break2.chr}:
+            if raw_read.reference_name != evidence.break1.chr and raw_read.reference_name != evidence.break2.chr:
                 continue
             read = evidence.standardize_read(raw_read)
             read.cigar = _cigar.merge_internal_events(