[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


> 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

More information about the Bioperl-l mailing list