t/SeqIO/kegg.t

   1 # -*-Perl-*- Test Harness script for Bioperl
   2 # $Id$
   3
   4 use strict;
   5
   6 BEGIN {
   7         use lib '.';
   8     use Bio::Root::Test;
   9
  10     test_begin(-tests => 16);
  11
  12         use_ok('Bio::SeqIO::kegg');
  13 }
  14
  15 my $verbose = test_debug();
  16
  17 my $io = Bio::SeqIO->new(-format => 'kegg',
  18                                                                  -verbose => $verbose,
  19                                                                  -file => test_input_file('AHCYL1.kegg'));
  20 isa_ok($io, 'Bio::SeqIO');
  21 my $kegg = $io->next_seq();
  22 ok($kegg);
  23 isa_ok($kegg, 'Bio::Seq::RichSeq');
  24 is($kegg->accession, '10768');
  25 is($kegg->display_id, 'AHCYL1');
  26 is($kegg->alphabet, 'dna');
  27 is($kegg->seq,'atgtcgatgcctgacgcgatgccgctgcccggggtcggggaggagctgaagcaggccaaggagatcgaggacgccgagaagtactccttcatggccaccgtcaccaaggcgcccaagaagcaaatccagtttgctgatgacatgcaggagttcaccaaattccccaccaaaactggccgaagatctttgtctcgctcgatctcacagtcctccactgacagctacagttcagctgcatcctacacagatagctctgatgatgaggtttctccccgagagaagcagcaaaccaactccaagggcagcagcaatttctgtgtgaagaacatcaagcaggcagaatttggacgccgggagattgagattgcagagcaagacatgtctgctctgatttcactcaggaaacgtgctcagggggagaagcccttggctggtgctaaaatagtgggctgtacacacatcacagcccagacagcggtgttgattgagacactctgtgccctgggggctcagtgccgctggtctgcttgtaacatctactcaactcagaatgaagtagctgcagcactggctgaggctggagttgcagtgttcgcttggaagggcgagtcagaagatgacttctggtggtgtattgaccgctgtgtgaacatggatgggtggcaggccaacatgatcctggatgatgggggagacttaacccactgggtttataagaagtatccaaacgtgtttaagaagatccgaggcattgtggaagagagcgtgactggtgttcacaggctgtatcagctctccaaagctgggaagctctgtgttccggccatgaacgtcaatgattctgttaccaaacagaagtttgataacttgtactgctgccgagaatccattttggatggcctgaagaggaccacagatgtgatgtttggtgggaaacaagtggtggtgtgtggctatggtgaggtaggcaagggctgctgtgctgctctcaaagctcttggagcaattgtctacattaccgaaatcgaccccatctgtgctctgcaggcctgcatggatgggttcagggtggtaaagctaaatgaagtcatccggcaagtcgatgtcgtaataacttgcacaggaaataagaatgtagtgacacgggagcacttggatcgcatgaaaaacagttgtatcgtatgcaatatgggccactccaacacagaaatcgatgtgaccagcctccgcactccggagctgacgtgggagcgagtacgttctcaggtggaccatgtcatctggccagatggcaaacgagttgtcctcctggcagagggtcgtctactcaatttgagctgctccacagttcccacctttgttctgtccatcacagccacaacacaggctttggcactgatagaactctataatgcacccgaggggcgatacaagcaggatgtgtacttgcttcctaagaaaatggatgaatacgttgccagcttgcatctgccatcatttgatgcccaccttacagagctgacagatgaccaagcaaaatatctgggactcaacaaaaatgggccattcaaacctaattattacagatactaa');
  28 ok($kegg->translate->seq);
  29
  30 is(($kegg->annotation->get_Annotations('description'))[0]->text,
  31    'S-adenosylhomocysteine hydrolase-like 1 [EC:3.3.1.1]');
  32
  33 is(($kegg->annotation->get_Annotations('pathway'))[0]->text,
  34    'Metabolism; Amino Acid Metabolism; Methionine metabolism');
  35
  36 is( (grep {$_->database eq 'KO'}
  37      $kegg->annotation->get_Annotations('dblink'))[0]->comment,
  38     'adenosylhomocysteinase' );
  39
  40 is( (grep {$_->database eq 'PATH'}
  41      $kegg->annotation->get_Annotations('dblink'))[0]->primary_id,
  42     'hsa00271' );
  43
  44 is( ($kegg->annotation->get_Annotations('aa_seq'))[0]->text,
  45 'MSMPDAMPLPGVGEELKQAKEIEDAEKYSFMATVTKAPKKQIQFADDMQEFTKFPTKTGRRSLSRSISQSSTDSYSSAASYTDSSDDEVSPREKQQTNSKGSSNFCVKNIKQAEFGRREIEIAEQDMSALISLRKRAQGEKPLAGAKIVGCTHITAQTAVLIETLCALGAQCRWSACNIYSTQNEVAAALAEAGVAVFAWKGESEDDFWWCIDRCVNMDGWQANMILDDGGDLTHWVYKKYPNVFKKIRGIVEESVTGVHRLYQLSKAGKLCVPAMNVNDSVTKQKFDNLYCCRESILDGLKRTTDVMFGGKQVVVCGYGEVGKGCCAALKALGAIVYITEIDPICALQACMDGFRVVKLNEVIRQVDVVITCTGNKNVVTREHLDRMKNSCIVCNMGHSNTEIDVTSLRTPELTWERVRSQVDHVIWPDGKRVVLLAEGRLLNLSCSTVPTFVLSITATTQALALIELYNAPEGRYKQDVYLLPKKMDEYVASLHLPSFDAHLTELTDDQAKYLGLNKNGPFKPNYYRY');
  46
  47 is( ($kegg->annotation->get_Annotations('position'))[0]->text, '1:join(26813..26932,50794..50905,52974..53117,54122..54222,54657..
  48             54759,56523..56617,57185..57291,58104..58220,58427..58490,59255..
  49             59343,59706..59776,60133..60227,60312..60410,60811..60879,61308..
  50             61386,62491..62611,63434..63440)');
  51
  52 is($kegg->species->common_name , 'H.sapiens');