[Bioperl-l] Genbank LOCUS entry

Matthias Adams matze at bioinf.uni-leipzig.de
Wed Feb 2 17:30:08 EST 2005


On Wed, Feb 02, 2005 at 01:52:01PM -0800, Hilmar Lapp wrote:
> You do get the LOCUS name - what you are asking for is the entire LOCUS 
> line. There no method for the entire line; but if that's all you need
"The id is equivalent to the LOCUS field of the GenBank/EMBL databanks"
(from http://doc.bioperl.org/releases/bioperl-1.4/Bio/Seq.html,
description of the display_id method)

> 	$ grep "^LOCUS" < my-genbank-file
> 
> on the command line.
As shown in the little example skript the Perl::Seq object is obtained
from the ncbi database, there is no my-genbank-file. Your solution is a
workaround. It's looks a little ugly to me to save it into a genbank file
an grep the LOCUS line, back to the same perl script...

Matze.





More information about the Bioperl-l mailing list