[Bioperl-l] Bio::SeqIO::swiss.pm - gene names
Corradi, John
John.Corradi@astrazeneca.com
Wed, 7 Feb 2001 16:35:57 -0500
I checked CVS and this has not been fixed. SwissProt entries list gene name
synonyms on the same line, delimited by " OR ". The parser in swiss.pm is
just getting the first symbol listed:
#Gene name
elsif (/^GN\s+(\S+)/) {
$seq->annotation->gene_name($1);
}
This works just fine for me:
elsif(/^GN\s+([^\.]+)/) {
$seq->annotation->gene_name($1);
}
Of course, you then have to split the string to get all the names. Maybe
annotation->{'gene_name'} should be an array reference?
Cheers,
John