We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent 2934805 commit 9031f99Copy full SHA for 9031f99
1 file changed
malariagen_data/anoph/genome_features.py
@@ -430,8 +430,8 @@ def canonical_transcript(
430
df_exons = df_exons[df_exons["type"] == "exon"].sort_values("start")
431
432
if len(df_exons) > 0:
433
- # Calculate total transcribed length
434
- exon_lengths = (df_exons["end"] - df_exons["start"]).sum()
+ # Calculate total transcribed length (1-based inclusive coordinates)
+ exon_lengths = (df_exons["end"] - df_exons["start"] + 1).sum()
435
transcript_lengths[transcript_id] = exon_lengths
436
debug(f" {transcript_id}: {len(df_exons)} exons, {exon_lengths} bp")
437
0 commit comments