[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