From bee369f86cde4159ba9aca296aafa47f2c1fdd97 Mon Sep 17 00:00:00 2001 From: Andrey Prjibelski Date: Thu, 11 Apr 2024 00:29:52 +0300 Subject: [PATCH] revert to normal profile ranges --- src/long_read_profiles.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/long_read_profiles.py b/src/long_read_profiles.py index 5dc4168d..1df28ef7 100644 --- a/src/long_read_profiles.py +++ b/src/long_read_profiles.py @@ -233,10 +233,10 @@ def construct_profile(self, sorted_blocks, polya_position=-1, polyt_position=-1) exon_profile[i] = -2 start_pos = 0 - while start_pos < len(exon_profile) and exon_profile[start_pos] < 1: + while start_pos < len(exon_profile) and exon_profile[start_pos] == 0: start_pos += 1 end_pos = len(exon_profile) - 1 - while end_pos >= 0 and exon_profile[end_pos] < 1: + while end_pos >= 0 and exon_profile[end_pos] == 0: end_pos -= 1 return MappedReadProfile(exon_profile, read_profile, read_exons, (start_pos, end_pos + 1))