# -*-Perl-*- Test Harness script for Bioperl
# $Id$
use strict;
BEGIN {
use lib '.';
use Bio::Root::Test;
test_begin(-tests => 47);
use_ok('Bio::SeqFeature::Primer');
use_ok('Bio::PrimarySeq');
}
my ($primer, $primer_seq, $location, $start, $end, $strand, $id, $tm, $tme, $template, $seq);
# Implied primer sequence
$template = Bio::Seq->new( -seq => 'AAAAACCCCCGGGGGTTTTT' );
ok $primer = Bio::SeqFeature::Primer->new( -start => 6, -end => 10 ), 'Implied primer sequence';
isa_ok $primer, 'Bio::SeqFeature::Primer';
isa_ok $primer, 'Bio::SeqFeature::SubSeq';
ok $template->add_SeqFeature($primer); # $primer->attach_seq($template);
ok $primer_seq = $primer->seq;
isa_ok $primer_seq, 'Bio::PrimarySeqI';
is $primer_seq->seq, 'CCCCC';
# Bio::PrimarySeq primer
$template = Bio::Seq->new( -seq => 'AAAAACCCCCGGGGGTTTTT' );
$seq = Bio::PrimarySeq->new(-seq => 'CTTTTCATTCTGACTGCAACG');
ok $primer = Bio::SeqFeature::Primer->new(-seq => $seq), 'PrimarySeq primer';
ok $template->add_SeqFeature($primer); # $primer->attach_seq($template);
ok $primer_seq = $primer->seq;
isa_ok $primer_seq, 'Bio::PrimarySeqI';
is $primer_seq->seq, 'CTTTTCATTCTGACTGCAACG';
# Initialize with a sequence string
ok $primer = Bio::SeqFeature::Primer->new(
-seq => 'CTTTTCATTCTGACTGCAACG',
-start => 3,
-id => 'myid',
);
is $primer->start, 3;
ok $primer_seq = $primer->seq;
is $primer_seq->isa('Bio::PrimarySeqI'), 1;
is $primer_seq->seq, 'CTTTTCATTCTGACTGCAACG';
is $primer_seq->id, 'myid';
is $primer->primary_tag, 'Primer';
ok $primer->display_name('test');
is $primer->display_name, 'test';
# Coordinates
ok $primer->start(2);
is $primer->start, 2;
ok $primer->end(19);
is $primer->end, 19;
ok $primer->strand(-1);
is $primer->strand, -1;
ok $location = $primer->location;
isa_ok $location, 'Bio::LocationI';
# Melting temperatures
ok $tm = $primer->Tm;
is int($tm), 52;
ok $tm = $primer->Tm(-salt => 0.05, -oligo => 0.0000001);
ok $tme = $primer->Tm_estimate;
is int($tme), 58;
ok $tm = $primer->Tm_estimate(-salt => 0.05);
# Legacy
# * initializing with -sequence
# * passing a string to location()
{
local $SIG{'__WARN__'} = sub { }; # Silence deprecation warnings
ok $primer = Bio::SeqFeature::Primer->new(
-sequence => 'CTTTTCATTCTGACTGCAACG',
);
ok $location = $primer->location('3,25');
is $location, '3,25';
}
# Chad's tests
$seq = Bio::Seq->new(
-seq => 'gcatcgatctagctagcta' ,
-id => 'chads_nifty_sequence',
);
$primer = Bio::SeqFeature::Primer->new(
-seq => $seq,
-TARGET => '5,3'
);
isa_ok $primer, 'Bio::SeqFeature::Primer';
is $primer->seq->display_id, 'chads_nifty_sequence';
is $primer->seq->seq, 'gcatcgatctagctagcta';
$primer = Bio::SeqFeature::Primer->new(
-seq => 'aaaaaacgatcgatcgtagctagct',
-id => 'chads_nifty_primer',
-TARGET => '5,3',
);
isa_ok $primer, 'Bio::SeqFeature::Primer';
isa_ok $primer->seq(), 'Bio::PrimarySeq';
is $primer->seq->id, 'chads_nifty_primer';
is $primer->seq->seq, 'aaaaaacgatcgatcgtagctagct';