t/Tools/Sigcleave.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 => 18);
  10
  11         use_ok('Bio::PrimarySeq');
  12         use_ok('Bio::Tools::Sigcleave');
  13 }
  14
  15 #load n-terminus of MGR5_HUMAN as test seq
  16 my $protein = "MVLLLILSVLLLKEDVRGSAQSSERRVVAHMPGDIIIGALFSVHHQPTVDKVHERKCGAVREQYGI";
  17
  18 ok my $seq= Bio::PrimarySeq->new(-seq => $protein);
  19
  20 ok my $sig = Bio::Tools::Sigcleave->new();
  21 ok $sig->seq($seq);
  22 ok my $sout = $sig->seq;
  23 is $sout->seq, $protein;
  24 is $sig->threshold, 3.5;
  25 is $sig->threshold(5), 5;
  26 is $sig->matrix, 'eucaryotic';
  27 is $sig->matrix('procaryotic'), 'procaryotic';
  28 is $sig->matrix('eucaryotic'), 'eucaryotic';
  29
  30 like $sig->pretty_print, qr/Maximum score 7/;
  31 ok my %results = $sig->signals;
  32
  33 is $results{9}, 5.2, "unable to get raw sigcleave results";
  34
  35
  36 $sig = Bio::Tools::Sigcleave->new(-seq=>$protein,
  37                                  -threshold=>5);
  38 ok %results = $sig->signals;
  39 is $results{9}, 5.2, "unable to get raw sigcleave results";
  40 is $sig->result_count, 5;