2
## Bioperl Test Harness Script for Modules
3
## $Id: Gel.t,v 1.1 2006/01/27 10:52:04 heikki Exp $
8
# to handle systems with no installed Test module
9
# we include the t dir (where a copy of Test.pm is located)
11
eval { require Test; };
19
use Bio::Restriction::Analysis;
25
my $seq1 = Bio::PrimarySeq->new(-id=>'groundhog day',
26
-seq=>'AAAAAAAAAGAATTCTTTTTTTTTTTTTTGAATTCGGGGGGGGGGGGGGGGGGGG');
29
my $ra=Bio::Restriction::Analysis->new(-seq=>$seq1);
30
ok my @cuts = $ra->fragments('EcoRI'), 3;
33
ok my $gel = Bio::Tools::Gel->new(-seq=>\@cuts,-dilate=>10);
34
ok my %bands = $gel->bands;
35
my @bands = (26, 27, 30);
37
foreach my $band (sort {$b <=> $a} keys %bands){
38
#print $band,"\t", sprintf("%.1f", $bands{$band}), "\n";
39
ok $bands[$c], sprintf("%.0f", $bands{$band});