From Code to Community: Sponsoring The Perl and Raku Conference 2025 Learn more

SENDU / bioperl-1.5.2_004-RC / t / largepseq.t 

            

              
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
86
87
88
89
90
91

              
# -*-Perl-*-
## Bioperl Test Harness Script for Modules
## $Id: largepseq.t,v 1.9 2001/11/02 18:48:09 jason Exp $
# Before `make install' is performed this script should be runnable with
# `make test'. After `make install' it should work as `perl test.t'
use strict;
BEGIN {     
    # to handle systems with no installed Test module
    # we include the t dir (where a copy of Test.pm is located)
    # as a fallback
    eval { require Test; };
    if( $@ ) {
        use lib 't';
    }
    use Test;
    plan tests => 22;
}
use Bio::Seq::LargePrimarySeq;
use Bio::Seq::LargeSeq;
use Bio::Location::Simple;
use Bio::Location::Fuzzy;
use Bio::Location::Split;
my $pseq = Bio::Seq::LargePrimarySeq->new();
ok $pseq;
$pseq->add_sequence_as_string('ATGGGGTGGGGTGAAACCCTTTGGGGGTGGGGTAAAT');
$pseq->add_sequence_as_string('GTTTGGGGTTAAACCCCTTTGGGGGGT');
ok $pseq->display_id('hello'), 'hello';
ok $pseq->seq, 'ATGGGGTGGGGTGAAACCCTTTGGGGGTGGGGTAAATGTTTGGGGTTAAACCCCTTTGGGGGGT' , "Sequence is " . $pseq->seq;
ok $pseq->subseq(3,7), 'GGGGT', "Subseq is ".$pseq->subseq(3,7);
my $location = new Bio::Location::Simple(-start => 4, -end => 8,
                                         -strand => 1);
ok($pseq->subseq($location), 'GGGTG');
my $splitlocation = new Bio::Location::Split;
$splitlocation->add_sub_Location( new Bio::Location::Simple('-start' => 1,
                                                            '-end'   => 15,
                                                            '-strand' => 1));
$splitlocation->add_sub_Location( new Bio::Location::Simple('-start' => 21,
                                                            '-end'   => 27,
                                                            '-strand' => -1));
ok( $pseq->subseq($splitlocation), 'ATGGGGTGGGGTGAACCCCCAA');
my $fuzzy = new Bio::Location::Fuzzy(-start => '<10',
                                     -end   => '18',
                                     -strand => 1);
ok( $pseq->subseq($fuzzy), 'GGTGAAACC');
ok($pseq->trunc(8,15)->seq, 'GGGGTGAA', 
    'trunc seq was ' . $pseq->trunc(8,15)->seq);
ok $pseq->alphabet('dna'), 'dna'; # so translate will not complain
ok $pseq->translate()->seq, 'MGWGETLWGWGKCLGLNPFGG';
my $seq = new Bio::Seq::LargeSeq(-primaryseq => $pseq );
ok $seq->display_id('hello'), 'hello';
ok $seq->seq, 'ATGGGGTGGGGTGAAACCCTTTGGGGGTGGGGTAAATGTTTGGGGTTAAACCCCTTTGGGGGGT' , "Sequence is " . $seq->seq;
ok $seq->subseq(3,7), 'GGGGT', "Subseq is ".$seq->subseq(3,7);
ok ($seq->trunc(8,15)->seq, 'GGGGTGAA', 
    'trunc seq was ' . $seq->trunc(8,15)->seq);
ok $seq->alphabet('dna'), 'dna'; # so translate will not complain
ok $seq->translate()->seq, 'MGWGETLWGWGKCLGLNPFGG';
$seq = new Bio::Seq::LargeSeq( -display_id => 'hello');
$seq->seq('ATGGGGTGGGGT');
ok $seq->display_id, 'hello';
ok $seq->seq, 'ATGGGGTGGGGT' , "Sequence is " . $seq->seq;
ok $seq->subseq(3,7), 'GGGGT', "Subseq is ".$seq->subseq(3,7);
ok ($seq->trunc(8,12)->seq, 'GGGGT', 
    'trunc seq was ' . $seq->trunc(8,12)->seq);
ok $seq->alphabet('dna'), 'dna'; # so translate will not complain
ok $seq->translate()->seq, 'MGWG';