Skip to content
Snippets Groups Projects
Commit 58f103c6 authored by Magali Ruffier's avatar Magali Ruffier
Browse files

ENSCORESW-232: code does not skip proteins with a *, B, Z and X any more

parent 467a594e
No related branches found
No related tags found
No related merge requests found
......@@ -146,13 +146,10 @@ sub dump_translation {
for my $dbid (@translation_ids) {
my $translation = $ta->fetch_by_dbID($dbid);
my $peptide_seq = $translation->seq();
if ( !($peptide_seq =~ m/[BZX]/ig ) ) {
if ( $peptide_seq !~ /\n$/ ) {
$peptide_seq .= "\n";
}
$peptide_seq =~ s/\*$//;
print TMP ">$dbid\n$peptide_seq";
if ( $peptide_seq !~ /\n$/ ) {
$peptide_seq .= "\n";
}
print TMP ">$dbid\n$peptide_seq";
}
close(TMP);
}
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment