diff options
| -rwxr-xr-x | varscan2codon.pl | 2 |
1 files changed, 2 insertions, 0 deletions
diff --git a/varscan2codon.pl b/varscan2codon.pl index 8cd36c6..71d880f 100755 --- a/varscan2codon.pl +++ b/varscan2codon.pl @@ -61,6 +61,8 @@ while (<$variants>) { } my $codon = substr $genome, $start+$protein_start, 3; + die "$. - couldn't find amino acid for codon '$codon'\n" unless ($codon2aa{$codon}); + my $orig_aa = $aa2letter{$codon2aa{$codon}}; substr $codon, $offset, 1, $new_nt; my $new_aa = $aa2letter{$codon2aa{$codon}}; |