summaryrefslogtreecommitdiffstats
path: root/varscan2codon.pl
diff options
context:
space:
mode:
authorBrian Cully <bjc@spork.org>2024-01-19 20:33:00 -0500
committerBrian Cully <bjc@spork.org>2024-01-19 20:33:00 -0500
commit337f04c364a2c291a4df23c31f5aff4070d08b47 (patch)
tree7fcf048e7743e8aaf6b823cce52a67a69d505037 /varscan2codon.pl
parent882b94f56de360933d712d601ae01cda520b80cd (diff)
downloadfav-337f04c364a2c291a4df23c31f5aff4070d08b47.tar.gz
fav-337f04c364a2c291a4df23c31f5aff4070d08b47.zip
don't bail on nucleotide check errors
print an error message and keep going
Diffstat (limited to 'varscan2codon.pl')
-rwxr-xr-xvarscan2codon.pl2
1 files changed, 1 insertions, 1 deletions
diff --git a/varscan2codon.pl b/varscan2codon.pl
index e45da77..6352a24 100755
--- a/varscan2codon.pl
+++ b/varscan2codon.pl
@@ -57,7 +57,7 @@ while (<$variants>) {
# verify the nucleotide at $start is what's expected from the csv.
my $check_nt = substr $genome, $start+$protein_start+$offset, 1;
if ($check_nt ne $orig_nt) {
- die "$name: nucleotide at position $pos is $check_nt, but expected $orig_nt.\n";
+ print STDERR "warning: $name (line $.): nucleotide at position $pos is $check_nt, but expected $orig_nt.\n";
}
my $codon = substr $genome, $start+$protein_start, 3;