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
|
# -*-Perl-*- Test Harness script for Bioperl
# $Id: Gel.t 15112 2008-12-08 18:12:38Z sendu $
use strict;
BEGIN {
use lib '.';
use Bio::Root::Test;
test_begin(-tests => 9);
use_ok('Bio::PrimarySeq');
use_ok('Bio::Restriction::Analysis');
use_ok('Bio::Tools::Gel');
}
my $seq1 = Bio::PrimarySeq->new(-id=>'groundhog day',
-seq=>'AAAAAAAAAGAATTCTTTTTTTTTTTTTTGAATTCGGGGGGGGGGGGGGGGGGGG');
my $ra=Bio::Restriction::Analysis->new(-seq=>$seq1);
is my @cuts = $ra->fragments('EcoRI'), 3;
ok my $gel = Bio::Tools::Gel->new(-seq=>\@cuts,-dilate=>10);
ok my %bands = $gel->bands;
my @bands = (26, 27, 30);
my $c = 0;
foreach my $band (sort {$b <=> $a} keys %bands){
#print $band,"\t", sprintf("%.1f", $bands{$band}), "\n";
is $bands[$c], sprintf("%.0f", $bands{$band});
$c++;
}
|
