2
## Bioperl Test Harness Script for Modules
3
## $Id: largepseq.t,v 1.9 2001/11/02 18:48:09 jason Exp $
5
# Before `make install' is performed this script should be runnable with
6
# `make test'. After `make install' it should work as `perl test.t'
10
# to handle systems with no installed Test module
11
# we include the t dir (where a copy of Test.pm is located)
13
eval { require Test; };
20
use Bio::Seq::LargePrimarySeq;
21
use Bio::Seq::LargeSeq;
22
use Bio::Location::Simple;
23
use Bio::Location::Fuzzy;
24
use Bio::Location::Split;
26
my $pseq = Bio::Seq::LargePrimarySeq->new();
28
$pseq->add_sequence_as_string('ATGGGGTGGGGTGAAACCCTTTGGGGGTGGGGTAAAT');
29
$pseq->add_sequence_as_string('GTTTGGGGTTAAACCCCTTTGGGGGGT');
31
ok $pseq->display_id('hello'), 'hello';
33
ok $pseq->seq, 'ATGGGGTGGGGTGAAACCCTTTGGGGGTGGGGTAAATGTTTGGGGTTAAACCCCTTTGGGGGGT' , "Sequence is " . $pseq->seq;
35
ok $pseq->subseq(3,7), 'GGGGT', "Subseq is ".$pseq->subseq(3,7);
36
my $location = new Bio::Location::Simple(-start => 4, -end => 8,
38
ok($pseq->subseq($location), 'GGGTG');
40
my $splitlocation = new Bio::Location::Split;
42
$splitlocation->add_sub_Location( new Bio::Location::Simple('-start' => 1,
46
$splitlocation->add_sub_Location( new Bio::Location::Simple('-start' => 21,
50
ok( $pseq->subseq($splitlocation), 'ATGGGGTGGGGTGAACCCCCAA');
52
my $fuzzy = new Bio::Location::Fuzzy(-start => '<10',
56
ok( $pseq->subseq($fuzzy), 'GGTGAAACC');
59
ok($pseq->trunc(8,15)->seq, 'GGGGTGAA',
60
'trunc seq was ' . $pseq->trunc(8,15)->seq);
63
ok $pseq->alphabet('dna'), 'dna'; # so translate will not complain
64
ok $pseq->translate()->seq, 'MGWGETLWGWGKCLGLNPFGG';
67
my $seq = new Bio::Seq::LargeSeq(-primaryseq => $pseq );
69
ok $seq->display_id('hello'), 'hello';
71
ok $seq->seq, 'ATGGGGTGGGGTGAAACCCTTTGGGGGTGGGGTAAATGTTTGGGGTTAAACCCCTTTGGGGGGT' , "Sequence is " . $seq->seq;
73
ok $seq->subseq(3,7), 'GGGGT', "Subseq is ".$seq->subseq(3,7);
74
ok ($seq->trunc(8,15)->seq, 'GGGGTGAA',
75
'trunc seq was ' . $seq->trunc(8,15)->seq);
77
ok $seq->alphabet('dna'), 'dna'; # so translate will not complain
78
ok $seq->translate()->seq, 'MGWGETLWGWGKCLGLNPFGG';
80
$seq = new Bio::Seq::LargeSeq( -display_id => 'hello');
81
$seq->seq('ATGGGGTGGGGT');
82
ok $seq->display_id, 'hello';
84
ok $seq->seq, 'ATGGGGTGGGGT' , "Sequence is " . $seq->seq;
86
ok $seq->subseq(3,7), 'GGGGT', "Subseq is ".$seq->subseq(3,7);
87
ok ($seq->trunc(8,12)->seq, 'GGGGT',
88
'trunc seq was ' . $seq->trunc(8,12)->seq);
90
ok $seq->alphabet('dna'), 'dna'; # so translate will not complain
91
ok $seq->translate()->seq, 'MGWG';