1 #!/usr/local/bin/perl -w
8 die("must define a valid seqfile to read") unless ( defined $seqfile && -r
$seqfile);
10 my $seqio = new Bio
::SeqIO
(-format
=> 'genbank',
13 while( my $seq = $seqio->next_seq ) {
15 # defined a default name
16 my $fname = sprintf("%s.gff", $seq->display_id || "seq-$count");
17 my $gffout = new Bio
::Tools
::GFF
(-file
=> ">$fname" ,
20 foreach my $feature ( $seq->top_SeqFeatures() ) {
21 $gffout->write_feature($feature);