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
|
# -*-Perl-*- Test Harness script for Bioperl
# $Id: Utilities.t 16090 2009-09-15 21:57:56Z cjfields $
use strict;
BEGIN {
use lib '.';
use Bio::Root::Test;
test_begin(-tests => 13);
use_ok('Bio::Align::Utilities', qw(:all));
use_ok('Bio::SimpleAlign');
use_ok('Bio::PrimarySeq');
use_ok('Bio::LocatableSeq');
use_ok('Bio::AlignIO');
}
my $DEBUG = test_debug();
my $aa_align = Bio::SimpleAlign->new();
$aa_align->add_seq(Bio::LocatableSeq->new(-id => "n1", -seq => "MLIDVG-MLVLR"));
$aa_align->add_seq(Bio::LocatableSeq->new(-id => "n2", -seq => "MLIDVRTPLALR"));
$aa_align->add_seq(Bio::LocatableSeq->new(-id => "n3", -seq => "MLI-VR-SLALR"));
my %dnaseqs = ();
$dnaseqs{'n1'} = Bio::PrimarySeq->new(-id => "n1", -seq => 'atgctgatagacgtaggcatgctagtactgaga');
$dnaseqs{'n2'} = Bio::PrimarySeq->new(-id => "n2", -seq => 'atgctgatcgacgtacgcaccccgctagcactcaga');
$dnaseqs{'n3'} = Bio::PrimarySeq->new(-id => "n3", -seq => 'atgttgattgtacgctcgcttgcacttaga');
my $dna_aln;
ok( $dna_aln = &aa_to_dna_aln($aa_align, \%dnaseqs));
if( $DEBUG ) {
Bio::AlignIO->new(-format=>'clustalw')->write_aln($dna_aln);
}
is $dna_aln->length, 36;
is $dna_aln->num_residues, 99;
is $dna_aln->num_sequences, 3;
is $dna_aln->consensus_string(50), "atgctgat?gacgtacgc????cgctagcact?aga";
$dna_aln->verbose(-1);
my $replicates;
ok $replicates = &bootstrap_replicates($dna_aln,3);
is scalar @$replicates, 3;
my $repl_aln = pop @$replicates;
is $repl_aln->num_sequences, 3;
##use IO::String;
##use Bio::AlignIO;
##my $string;
##my $out = IO::String->new($string);
##
##my $strout = Bio::AlignIO->new(-fh => $out,'-format' => 'pfam');
##$strout->write_aln($repl_aln);
##is $string, "";
|
