Skip to content

Commit

Permalink
releave penalty score for terminal exon elongation, also relates to #270
Browse files Browse the repository at this point in the history
  • Loading branch information
andrewprzh committed Jan 9, 2025
1 parent a9f29a2 commit 22bedd8
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions src/long_read_assigner.py
Original file line number Diff line number Diff line change
Expand Up @@ -184,8 +184,8 @@ def select_similar_isoforms(self, combined_read_profile):
read_region = (read_split_exon_profile.read_features[0][0], read_split_exon_profile.read_features[-1][1])
for isoform_id, diff_introns in intron_matching_isoforms:
transcript_start, transcript_end = self.gene_info.transcript_region(isoform_id)
extra_left = 1 if read_region[0] + self.params.delta < transcript_start else 0
extra_right = 1 if read_region[1] - self.params.delta > transcript_end else 0
extra_left = 0.5 if read_region[0] + self.params.delta < transcript_start else 0
extra_right = 0.5 if read_region[1] - self.params.delta > transcript_end else 0
candidates.append((isoform_id, diff_introns + extra_right + extra_left))
# select isoforms that have similar number of potential inconsistencies
best_diff = min(candidates, key=lambda x: x[1])[1]
Expand Down

0 comments on commit 22bedd8

Please sign in to comment.