From 483ecf4c5e1952b6ec4dda43ced9598eb4ed75c7 Mon Sep 17 00:00:00 2001 From: Peter Vegh Date: Tue, 10 Dec 2024 18:41:46 +0000 Subject: [PATCH] Apply cutoff to asm length --- ediacara/Comparator.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/ediacara/Comparator.py b/ediacara/Comparator.py index 40e13a0..abdbf10 100644 --- a/ediacara/Comparator.py +++ b/ediacara/Comparator.py @@ -368,8 +368,8 @@ def interpret_asm(self): self.len_asm = len(self.asm) self.asm_ref_len_diff = self.len_asm - self.len_ref - length_band = 50 # bp, smaller errors should be picked up by coverage, VCF etc - if abs(self.asm_ref_len_diff) > length_band: + length_band = self.diff_cutoff # bp + if abs(self.asm_ref_len_diff) >= length_band: self.has_errors = True self.has_large_diff = True self.wrong_length = True