#...no - an internal stop codon error in the database...
#...no - an internal stop codon error in the database...
if($num_tga<$num_stops){
if($num_tga<$num_stops){
if($sourceeq'havana'){
if($sourceeq'havana'){
$support->log_warning("INTERNAL STOPS HAVANA: Transcript $tsi ($tname) has non \'$scodon\' stop codons [$mod_date]:\nSequence = $orig_seq\nStops at $positions)\n\n");
$support->log_warning("INTERNAL STOPS HAVANA: Transcript $tsi ($tname) from gene $gnamehas non \'$scodon\' stop codons [$mod_date]:\nSequence = $orig_seq\nStops at $positions)\n\n");
}
}
else{
else{
$support->log_warning("INTERNAL STOPS EXTERNAL: Transcript $tsi ($tname) has non \'$scodon\' stop codons[$mod_date]:\nSequence = $orig_seq\nStops at $positions)\n\n");
$support->log_warning("INTERNAL STOPS EXTERNAL: Transcript $tsi ($tname) from gene $gnamehas non \'$scodon\' stop codons[$mod_date]:\nSequence = $orig_seq\nStops at $positions)\n\n");
}
}
}
}
...
@@ -315,7 +316,7 @@ sub check_for_stops {
...
@@ -315,7 +316,7 @@ sub check_for_stops {
$support->log_verbose("Transcript $tsi ($tname) has potential selenocysteines but has been discounted by annotators:\n\t".$seen_transcripts->{$tsi}.") [$mod_date]\n");
$support->log_verbose("Transcript $tsi ($tname) has potential selenocysteines but has been discounted by annotators:\n\t".$seen_transcripts->{$tsi}.") [$mod_date]\n");
}
}
else{
else{
$support->log("POTENTIAL SELENO ($seq) in $tsi ($tname) found at $pos) [$mod_date]\n");
$support->log("POTENTIAL SELENO ($seq) in $tsi ($tname, gene $gname) found at $pos [$mod_date]\n");