[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