<![CDATA[
#!perl -w
# Read in a series of GenBank records
use strict;
use warnings;
use Boulder::Stream;
use Boulder::Genbank;
my $debug = grep(/-d/, @ARGV); # check for -d option
# parse a file of Genbank records
my $stream = new Boulder::Genbank(-accessor=>'File', -fetch => 'library.gb');
my $stone;
while ($stone = $stream->get) {
if ( $debug ) {
print "Accession is: ",$stone->Accession,"\n";
print "Definition is: ",$stone->Definition,"\n";
my $dna = $stone->Sequence;
print "The dna is ", length($dna), " bp long.\n";
} else {
$stream->write_record($stone);
}
}
]]>