1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
|
# -*-Perl-*- Test Harness script for Bioperl
# $Id: GenBank.t 15649 2009-04-17 15:08:36Z cjfields $
use strict;
BEGIN {
use lib '.';
use Bio::Root::Test;
test_begin(-tests => 21,
-requires_modules => [qw(IO::String
LWP::UserAgent
HTTP::Request::Common)],
-requires_networking => 1);
use_ok('Bio::DB::Query::GenBank');
use_ok('Bio::DB::GenBank');
}
my %expected_lengths = (
'MUSIGHBA1' => 408,
'AF303112' => 1611,
'AF041456' => 1156,
'AY080910' => 798,
'AY080909' => 1042,
'AF155220' => 1172,
'AF442768' => 2547,
);
my ($gb, $seq, $seqio, $seqin, $query);
# test query facility
ok $query = Bio::DB::Query::GenBank->new('-db' => 'nucleotide',
'-query' => 'Onchocerca volvulus[Organism]',
'-mindate' => '2002/1/1',
'-maxdate' => '2002/12/31'), 'Bio::DB::Query::GenBank';
SKIP: {
cmp_ok $query->count, '>', 0;
my @ids = $query->ids;
cmp_ok @ids, '>', 0;
is @ids, $query->count;
ok $gb = Bio::DB::GenBank->new('-delay' => 0);
eval {$seqio = $gb->get_Stream_by_query($query);};
skip "Couldn't connect to complete GenBank query tests. Skipping those tests", 5 if $@;
my $done = 0;
while (my $s = $seqio->next_seq) {
is $s->length, $expected_lengths{$s->display_id}, $s->display_id;
undef $gb; # test the case where the db is gone,
# but the pipeline should remain until seqio goes away
$done++;
}
skip('No seqs returned', 5) if !$done;
is $done, 4;
}
$seq = $seqio = undef;
# test query facility (again)
ok $query = Bio::DB::Query::GenBank->new('-db' => 'nucleotide',
'-ids' => [qw(J00522 AF303112 2981014)]);
SKIP: {
cmp_ok $query->count, '>', 0;
my @ids = $query->ids;
cmp_ok @ids, '>', 0;
is @ids, $query->count;
$gb = Bio::DB::GenBank->new('-delay' => 0);
eval {$seqio = $gb->get_Stream_by_query($query);};
skip "Couldn't connect to complete GenBank query tests. Skipping those tests: $@", 4 if $@;
my $done = 0;
while (my $s = $seqio->next_seq) {
is $s->length, $expected_lengths{$s->display_id}, $s->display_id;
$done++;
}
skip('No seqs returned', 4) if !$done;
is $done, 3;
$seqio->close(); # the key to preventing errors during make test, no idea why
}
$seq = $seqio = undef;
# and yet again, for bug 2133
$query = Bio::DB::Query::GenBank->new('-query' => 'AF303112',
'-ids' => [qw(J00522 AF303112 2981014)]);
is $query->query, 'J00522[PACC]|AF303112[PACC]|2981014[UID]';
|
