t/Tools/Lucy.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 => 22);
  10
  11     use_ok('Bio::Tools::Lucy');
  12 }
  13
  14 my @params = (adv_stderr => 1, seqfile => test_input_file('lucy.seq'), rev_desig => 'R');
  15 # Bio::Tools::Lucy will find .qual, .info, and .stderr files in this folder
  16
  17 my $lucyObj = Bio::Tools::Lucy->new(@params);
  18 isa_ok $lucyObj,'Bio::Tools::Lucy';
  19 ok $lucyObj->seqfile();
  20 $lucyObj->adv_stderr(1);
  21 my $stderr = $lucyObj->adv_stderr();
  22 ok $stderr;
  23 my $names =$lucyObj->get_sequence_names();
  24 ok $names;
  25 my $seq = shift @$names;
  26 is $seq, 'TaLr1010B10R';
  27 is $lucyObj->length_raw("$seq"), 1060;
  28 is $lucyObj->length_clear("$seq"), 420;
  29 is $lucyObj->start_clear("$seq"), 86;
  30 is $lucyObj->end_clear("$seq"), 505;
  31 ok $lucyObj->avg_quality("$seq");
  32 ok $lucyObj->full_length("$seq");
  33 ok $lucyObj->polyA("$seq");
  34 is $lucyObj->direction("$seq"), 'R';
  35 ok $lucyObj->per_GC("$seq");
  36 ok $lucyObj->sequence("$seq");
  37 ok $lucyObj->quality("$seq");
  38 my $seqObj = $lucyObj->get_Seq_Obj("$seq");
  39 ok $seqObj;
  40 my $seqObjs = $lucyObj->get_Seq_Objs();
  41 ok $seqObjs;
  42
  43 my $rejects = $lucyObj->get_rejects();
  44 ok $rejects;
  45 my ($key) = (sort keys %$rejects);
  46 is $key, 'TaLr1011A07R';
  47 is $rejects->{$key}, 'Q';