[Bioperl-l] retrieving top_SeqFeatures for RefSeq proteins fails

Torsten Seemann torsten.seemann at infotech.monash.edu.au
Tue Apr 11 22:37:20 EDT 2006


Eleni,

> my $factory = Bio::DB::GenBank->new(-no_redirect => 1);
> my $seq = $factory->get_Seq_by_acc(' NM_001008292');

You have a space character before the NM_ ?
Try removing it.

> returns:
> MSG: [gb| NM_001008292] is not a normal sequence entry but a RefSeq entry.
> Redirecting the request.

There's the space, right in the error message.

Is that the problem?

Brian's counter example didn't have the space.

-- 
Torsten Seemann
Victorian Bioinformatics Consortium, Monash University, Australia
http://www.vicbioinformatics.com/


More information about the Bioperl-l mailing list