Merge pull request #73 from bigbio/issue72

fixes issue #72
bigbio · Nov 15, 2022 · bbbbb34 · bbbbb34
2 parents 913b1e7 + 36f8b58
commit bbbbb34
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/pypgatk/cgenomes/cgenomes_proteindb.py b/pypgatk/cgenomes/cgenomes_proteindb.py
@@ -121,6 +121,8 @@ def get_mut_pro_seq(snp, seq):
 
                 if "Missense" in snp.type:
                     mut_aa = snp.aa_mut[-1]
+                    if not mut_aa.isalpha():
+                        return ''
                     index = int(positions[0]) - 1
                     mut_pro_seq = protein_seq[:index] + mut_aa + protein_seq[index + 1:]
                 elif "Nonsense" in snp.type:
@@ -146,6 +148,8 @@ def get_mut_pro_seq(snp, seq):
                     except ValueError:
                         return ''
                     mut_aa = snp.aa_mut[index + 1:]
+                    if not mut_aa.replace('*','').isalpha():
+                        return ''
                     if "deletion" in snp.type:
                         del_index1 = int(positions[0]) - 1
                         del_index2 = int(positions[1])