fix: return unknown when last codon is not stop codon

chrovis · Oct 14, 2024 · 2c0c034 · 2c0c034
1 parent bef787d
commit 2c0c034
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/src/varity/vcf_to_hgvs/protein.clj b/src/varity/vcf_to_hgvs/protein.clj
@@ -333,6 +333,11 @@
                    (count ref-exon-seq) (:name rg) (:name2 rg))
         {:type :unknown, :pos nil, :ref nil, :alt nil})
 
+    (not= (last ref-prot-seq) \*)
+    (do (log/warnf "Last codon is not stop codon: %s (%s, %s)"
+                   (str (last ref-prot-seq)) (:name rg) (:name2 rg))
+        {:type :unknown, :pos nil, :ref nil, :alt nil})
+
     :else
     (let [alt-prot-seq* (format-alt-prot-seq seq-info)
           ppos (protein-position pos rg)