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
|
# -*-Perl-*- Test Harness script for Bioperl
# $Id: Seg.t 15112 2008-12-08 18:12:38Z sendu $
use strict;
BEGIN {
use lib '.';
use Bio::Root::Test;
test_begin(-tests => 15);
use_ok('Bio::Tools::Seg');
}
my ($infile, $parser) ;
$infile = test_input_file('seg.out');
ok ($parser = Bio::Tools::Seg->new(-file=>$infile), 'parser defined') ;
my @feat;
while ( my $feat = $parser->next_result ) {
push @feat, $feat;
}
is scalar(@feat), 3;
# seq 0
#>LBL_0012(32-46) complexity=2.47 (12/2.20/2.50)
#gdggwtfegwggppe
# seq 1
#>LBL_0012(66-80) complexity=2.31 (12/2.20/2.50)
#kfssrasakavakks
# seq 2
#>LBL_0012(123-138) complexity=2.31 (12/2.20/2.50)
#svivsqsqgvvkgvgv
my $raa_testdata = [
[ 'LBL_0012', 32, 46, 2.47 ],
[ 'LBL_0012', 66, 80, 2.31 ],
[ 'LBL_0012', 123, 138, 2.31 ],
] ;
for (0..( scalar(@feat)-1 )) {
is ( $feat[$_]->seq_id, $raa_testdata->[$_]->[0], "seq id for seq $_ identified" ) ;
is ( $feat[$_]->start, $raa_testdata->[$_]->[1], "start for seq $_ identified" ) ;
is ( $feat[$_]->end, $raa_testdata->[$_]->[2], "end for seq $_ identified" ) ;
is ( $feat[$_]->score, $raa_testdata->[$_]->[3], "score for seq $_ identified" ) ;
}
|
