
Bug: with perl > 5.36, the compiler is more restrictive with errors. This package has implicit function definitions. The fix is to ignore the errors as vep, the main package for that, runs without issue. Also this package is no longer maintained We also add test that were missing from the package. Not all test can be run as they depends on other Perl modules in Bioper.
85 lines
2.7 KiB
Diff
85 lines
2.7 KiB
Diff
diff --git a/Bio/Ext/Align/test.pl b/Bio/Ext/Align/test.pl
|
|
index 72411f3..1deb77b 100755
|
|
--- a/Bio/Ext/Align/test.pl
|
|
+++ b/Bio/Ext/Align/test.pl
|
|
@@ -8,13 +8,10 @@ my $DEBUG = $ENV{'BIOPERLDEBUG'} || 0;
|
|
BEGIN {
|
|
eval { require Test; };
|
|
use Test;
|
|
- plan tests => 9;
|
|
+ plan tests => 4;
|
|
}
|
|
|
|
use Bio::Ext::Align;
|
|
-use Bio::Tools::dpAlign;
|
|
-use Bio::Seq;
|
|
-use Bio::AlignIO;
|
|
|
|
$loaded = 1;
|
|
ok(1); # modules loaded
|
|
@@ -34,64 +31,3 @@ $alb = &Bio::Ext::Align::Align_Sequences_ProteinSmithWaterman($seq1,$seq2,
|
|
$seq2->seq,15,50,STDERR) if $DEBUG;
|
|
|
|
|
|
-warn( "Testing Local Alignment case...\n") if $DEBUg;
|
|
-
|
|
-$alnout = new Bio::AlignIO(-format => 'pfam', -fh => \*STDERR);
|
|
-$aln = &Bio::Ext::Align::Align_DNA_Sequences("AATGCCATTGACGG",
|
|
- "CAGCCTCGCTTAG",3,-1,3,1,
|
|
- Bio::Tools::dpAlign::DPALIGN_LOCAL_MILLER_MYERS);
|
|
-
|
|
-$out = Bio::SimpleAlign->new();
|
|
-
|
|
-$out->add_seq(Bio::LocatableSeq->new(-seq => $aln->aln1,
|
|
- -start => $aln->start1,
|
|
- -end => $aln->end1,
|
|
- -id => "one"));
|
|
-
|
|
-$out->add_seq(Bio::LocatableSeq->new(-seq => $aln->aln2,
|
|
- -start => $aln->start2,
|
|
- -end => $aln->end2,
|
|
- -id => "two"));
|
|
-$alnout->write_aln($out) if $DEBUG;
|
|
-
|
|
-$aln = &Bio::Ext::Align::Align_Protein_Sequences("WLGQRNLVSSTGGNLLNVWLKDW","WMGNRNVVNLLNVWFRDW",0,
|
|
- Bio::Tools::dpAlign::DPALIGN_LOCAL_MILLER_MYERS);
|
|
-$out = Bio::SimpleAlign->new();
|
|
-ok($aln);
|
|
-
|
|
-$out->add_seq(Bio::LocatableSeq->new(-seq => $aln->aln1,
|
|
- -start => $aln->start1,
|
|
- -end => $aln->end1,
|
|
- -id => "one"));
|
|
-
|
|
-$out->add_seq(Bio::LocatableSeq->new(-seq => $aln->aln2,
|
|
- -start => $aln->start2,
|
|
- -end => $aln->end2,
|
|
- -id => "two"));
|
|
-$alnout->write_aln($out) if $DEBUG;
|
|
-ok(1);
|
|
-
|
|
-warn( "Testing Global Alignment case...\n") if $DEBUG;
|
|
-
|
|
-$factory = new Bio::Tools::dpAlign('-alg' => Bio::Tools::dpAlign::DPALIGN_GLOBAL_MILLER_MYERS);
|
|
-$s1 = new Bio::Seq(-id => "one", -seq => "AATGCCATTGACGG", -alphabet => 'dna');
|
|
-$s2 = new Bio::Seq(-id => "two", -seq => "CAGCCTCGCTTAG", -alphabet => 'dna');
|
|
-$aln = $factory->pairwise_alignment($s1, $s2);
|
|
-$alnout->write_aln($aln) if $DEBUG;
|
|
-$factory->align_and_show($s1, $s2) if $DEBUG;
|
|
-
|
|
-ok(1);
|
|
-
|
|
-$s1 = new Bio::Seq(-id => "one", -seq => "WLGQRNLVSSTGGNLLNVWLKDW",
|
|
- -alphabet => 'protein');
|
|
-$s2 = new Bio::Seq(-id => "two", -seq => "WMGNRNVVNLLNVWFRDW",
|
|
- -alphabet => 'protein');
|
|
-$aln = $factory->pairwise_alignment($s1, $s2);
|
|
-$alnout->write_aln($aln) if $DEBUG;
|
|
-$factory->align_and_show($s1, $s2) if $DEBUG;
|
|
-ok(1);
|
|
-
|
|
-$prof = $factory->sequence_profile($s1);
|
|
-warn( "Optimal Alignment Score = %d\n", $factory->pairwise_alignment_score($prof, $s2)) if $DEBUG;
|
|
-
|
|
-ok($factory->pairwise_alignment_score($prof,$s2),77);
|