[Bioperl-l] GFF synopsis code question

Sean Davis sdavis2 at mail.nih.gov
Thu Mar 18 07:30:38 EST 2004


I have been learning to use Bio::DB::GFF and had a question about the
synopsis code from Bio::DB::GFF.  The section of code is:

-------From Synopsis-------
  # pull out all transcript features
  my @transcripts = $segment->features('transcript');

  # for each transcript, total the length of the introns
  my %totals;
  for my $t (@transcripts) {
    my @introns = $t->Intron;
    $totals{$t->name} += $_->length foreach @introns;
  }

  # Sort the exons of the first transcript by position
  my @exons = sort {$a->start <=> $b->start} $transcripts[0]->Exon;

  # Get a region 1000 bp upstream of first exon
  my $upstream = $exons[0]->segment(-1000,0);

  # get its DNA
  my $dna = $upstream->dna;
-------------------------------

It seems here that $exons[0] is a feature, so does not have a segment
method?  Instead, one can get the segment of interest by:

my $upstream = $db->segment($exons[0],-1000=>0)

Am I mistaken?  I only bring it up as a documentation issue, as the code I
give works for me.

Sean 



More information about the Bioperl-l mailing list