use strict;
BEGIN {
use lib '.';
use Bio::Root::Test;
test_begin(-tests => 21);
use_ok 'Bio::PrimarySeq';
use_ok 'Bio::SeqFeature::Primer';
use_ok 'Bio::SeqFeature::Amplicon';
}
my ($amplicon, $amplicon_seq, $fwd_primer, $rev_primer, $template);
# Basic amplicon object
$amplicon = Bio::SeqFeature::Amplicon->new();
isa_ok $amplicon, 'Bio::SeqFeature::Amplicon';
isa_ok $amplicon, 'Bio::SeqFeature::SubSeq';
# Amplicon with explicit sequence (sequence string)
ok $amplicon = Bio::SeqFeature::Amplicon->new(
-seq => 'CCCCCAAAAAGGGGGTTTTT',
);
ok $amplicon_seq = $amplicon->seq;
isa_ok $amplicon_seq, 'Bio::PrimarySeq';
is $amplicon_seq->seq, 'CCCCCAAAAAGGGGGTTTTT';
# Amplicon with explicit sequence (sequence object)
$fwd_primer = Bio::SeqFeature::Primer->new(
-start => 1,
-end => 4,
-strand => 1
);
$rev_primer = Bio::SeqFeature::Primer->new(
-seq => 'GATTA',
-start => 16,
-end => 20,
-strand => -1
);
ok $amplicon = Bio::SeqFeature::Amplicon->new(
-start => 1,
-end => 20,
-seq => Bio::PrimarySeq->new( -seq => 'CCCCCAAAAAGGGGGTTTTT' ),
-fwd_primer => $fwd_primer,
);
ok $amplicon->rev_primer($rev_primer);
is_deeply $amplicon->fwd_primer(), $fwd_primer;
is_deeply $amplicon->rev_primer(), $rev_primer;
ok $amplicon_seq = $amplicon->seq;
isa_ok $amplicon_seq, 'Bio::PrimarySeq';
is $amplicon_seq->seq, 'CCCCCAAAAAGGGGGTTTTT';
# Amplicon with implicit sequence
$template = Bio::Seq->new( -seq => 'ATCGATCGATCCCCCAAAAAGGGGGTTTTTAGCTAGCTAT');
ok $amplicon = Bio::SeqFeature::Amplicon->new(
-start => 11,
-end => 30,
-strand => -1
);
ok $template->add_SeqFeature($amplicon);
ok $amplicon_seq = $amplicon->seq;
isa_ok $amplicon_seq, 'Bio::PrimarySeq';
is $amplicon_seq->seq, 'AAAAACCCCCTTTTTGGGGG';