[Bioperl-l] extracting genbank content

[Chris Fields]

Depends on what you need to do.  If you are just interested in pulling out certain bits of data from each record, using SeqIO is a good option.  But if you want to access the records as a flat database (not iteration, but indexed for fast access), use Bio::Index::GenBank or Bio::DB::Flat to make a simple flat file database and access them by ID.

chris

On May 17, 2010, at 7:42 AM, Ross KK Leung wrote:

> Dear all,
> 
> 
> 
> When there are more than one genbank records in a file, except by splitting
> the file into separate records, what can I do to transverse the records?
> 
> 
> 
> $obj=Bio::SeqIO->new(-file=>$gbfile,-format=>"genbank");
> 
> 
> $seqobj=$obj->next_seq();        
> 
> 
> 
> Do I just use another $obj->next_seq() so it will point to another record?
> 
> 
> 
> Thanks for your advice.
> 
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l