t/Spidey.t

   1 # -*-Perl-*- Test Harness script for Bioperl
   2 # $Id$
   3
   4 use strict;
   5
   6 BEGIN {
   7     use lib 't/lib';
   8     use BioperlTest;
   9
  10     test_begin(-tests => 26);
  11
  12         use_ok('Bio::Tools::Spidey::Results');
  13 }
  14
  15 ok my $spidey = Bio::Tools::Spidey::Results->new(-file=> test_input_file('spidey.noalignment'));
  16
  17 $spidey->close();
  18
  19 my $exonset = $spidey->next_exonset;
  20 ok(!defined($exonset));
  21
  22 $spidey = Bio::Tools::Spidey::Results->new(-file=> test_input_file('spidey.test1'));
  23 $exonset = $spidey->next_exonset;
  24 my @exons = $exonset->sub_SeqFeature();
  25 is @exons, 6;
  26
  27 is $spidey->genomic_dna_length(), 145732769;
  28 is $spidey->splicesites(), 4;
  29 is $spidey->est_coverage(), 100;
  30 is $spidey->overall_percentage_id(), 99.7;
  31 is $spidey->missing_mrna_ends(), 'neither';
  32
  33 is $exonset->seq_id(), 'lcl|chr2';
  34 is $exonset->start(), 36356457;
  35 is $exonset->end(), 36375798;
  36
  37 my $exon = 0;
  38 is $exons[$exon]->est_hit()->seq_id(), 'lcl|tmpseq_0';
  39 is $exons[$exon]->start(), 36375691;
  40 is $exons[$exon]->end(), 36375798;
  41 is $exons[$exon]->strand(), -1;
  42 is $exons[$exon]->est_hit()->start(), 1;
  43 is $exons[$exon]->est_hit()->end(), 108;
  44 is $exons[$exon]->donor(), 1;
  45 is $exons[$exon]->acceptor(), 0;
  46
  47 $exon = 1;
  48 is $exons[$exon]->start(), 36369345;
  49 is $exons[$exon]->end(), 36369492;
  50 is $exons[$exon]->est_hit()->start(), 109;
  51 is $exons[$exon]->est_hit()->end(), 256;
  52 is $exons[$exon]->donor(), 0;
  53 is $exons[$exon]->acceptor(), 1;
  54
  55 $spidey->close();