#!/usr/bin/perl use strict; use Bio::DB::GFF; my $db = Bio::DB::GFF->new(-user => 'bio', -pass => 'bio', -dsn => 'elegans', -aggregators => 'gene_model{coding_exon,5_UTR,3_UTR/CDS}'); my $gene_stream = $db->get_seq_stream('gene_model:curated'); while (my $gene = $gene_stream->next_seq) { print $gene->name,"\n"; for my $part ($gene->get_SeqFeatures) { print "\t",join("\t",$part->method,$part->start,$part->end),"\n"; } print "\n"; }