t/Tools/tRNAscanSE.t

   1 # -*-Perl-*- Test Harness script for Bioperl
   2 # $Id$
   3
   4 use strict;
   5
   6 BEGIN {
   7     use Bio::Root::Test;
   8
   9     test_begin(-tests => 14);
  10
  11         use_ok('Bio::Tools::tRNAscanSE');
  12 }
  13
  14 my $verbose = test_debug();
  15
  16 my $parser = Bio::Tools::tRNAscanSE->new(-verbose => $verbose,
  17                                          -file => test_input_file('yeast.tRNAscanSE'));
  18
  19 isa_ok($parser, 'Bio::Tools::tRNAscanSE') ;
  20
  21 my @genes;
  22 while( my $gene = $parser->next_prediction ) {
  23     push @genes, $gene;
  24 }
  25
  26 is (scalar(@genes), 287);
  27 is($genes[2]->seq_id, 'I', 'seq_id');
  28 my ($codon) = $genes[2]->get_tag_values('Codon');
  29 is($codon, 'TTG', 'codon');
  30 is($genes[2]->start, 181135, 'start');
  31 is($genes[2]->end, 181248, 'end');
  32 is($genes[2]->strand, 1, 'strand');
  33
  34 my @exons = $genes[2]->get_SeqFeatures ;
  35 is ( scalar(@exons), 2, 'exons' );
  36 is($exons[0]->end,181172, 'end' );
  37 is($exons[0]->start,$genes[2]->start, 'start');
  38 is($exons[1]->start,181205, 'start');
  39 is($exons[1]->end,$genes[2]->end, 'end');
  40 is($exons[0]->seq_id, $genes[2]->seq_id, 'seq_id');