t/SeqIO/gbxml.t

   1 use strict;
   2
   3 BEGIN {
   4    use Bio::Root::Test;
   5
   6    test_begin(-tests => 14,
   7               -requires_modules => [ qw(XML::SAX) ]
   8               );
   9
  10    use_ok('Bio::SeqIO::gbxml');
  11 }
  12
  13 my $verbose = test_debug();
  14
  15 my $in = Bio::SeqIO->new(-format  => 'gbxml',
  16                           -verbose => $verbose,
  17                           # This class can't parse dbEST data yet...
  18                           # -file    => test_input_file('roa1.gbxml'));
  19                           # So let's try a <GBSeq> file:
  20                           -file    => test_input_file('EG352462.gbxml'));
  21 isa_ok($in, 'Bio::SeqIO');
  22 my $seq = $in->next_seq();
  23 is($seq->molecule,   'mRNA',                                            'molecule');
  24 is($seq->alphabet,   'dna',                                             'alphabet');
  25 is($seq->primary_id,  116038450,                                        'primary_id');
  26 is($seq->display_id, 'EG352462',                                        'display_id');
  27 is($seq->version,     1,                                                'version');
  28 is($seq->is_circular, 0,                                                'is_circular');
  29
  30 is(substr($seq->description, 0, 10), 'SAAH-aad23',                      'description');
  31 is(substr($seq->seq,         0, 10), 'aataaaatta',                      'sequence');
  32
  33 my @class = $seq->species->classification;
  34 is($class[$#class],'Eukaryota',                                         'classification');
  35
  36 my ($feat) = $seq->get_SeqFeatures;
  37 is_deeply([ $feat->get_tag_values('clone_lib') ], [ 'Agen 0058' ],      'feat - clone_lib');
  38 is_deeply([ $feat->get_tag_values('db_xref') ],   [ 'taxon:79327' ],    'feat - db_xref');
  39 is_deeply([ $feat->get_tag_values('lab_host') ],  [ 'DH10B cells' ],    'feat - lab_host');
  40
  41
  42