counsyl#57 - c. dup converted back to delins

SACGF · May 25, 2023 · 7591c13 · 7591c13
1 parent 116217d
commit 7591c13
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/pyhgvs/__init__.py b/pyhgvs/__init__.py
@@ -155,8 +155,9 @@ def hgvs_justify_indel(chrom, offset, ref, alt, strand, genome):
         return chrom, offset, ref, alt
 
     # Get genomic sequence around the lesion.
-    start = max(offset - 100, 0)
-    end = offset + 100
+    size = max(len(ref), len(alt)) + 1
+    start = max(offset - size, 0)
+    end = offset + size
     seq = str(genome[str(chrom)][start - 1:end]).upper()
     cds_offset = offset - start