my
$error
;
my
$DEBUG
=
$ENV
{
'BIOPERLDEBUG'
} || 0;
BEGIN {
$error
= 0;
if
( $@ ) {
}
plan
tests
=> 12;
}
ok(1);
my
$parser
= new Bio::Tools::EMBOSS::Palindrome
(
-file
=> Bio::Root::IO->catfile(
qw( t data humts1.pal)
));
my
$seq
=
$parser
->next_seq;
ok(
$seq
);
ok(
$seq
->display_id,
'HUMTS1'
);
ok(
$seq
->
length
, 18596);
my
@features
=
$seq
->get_SeqFeatures();
ok(
scalar
@features
, 23);
ok(
$features
[0]->feature1->start, 126);
ok(
$features
[0]->feature1->end, 142);
ok(
$features
[0]->feature1->strand, 1);
ok(
$features
[0]->feature1->seq_id,
'HUMTS1'
);
ok(
$features
[0]->feature2->start, 201);
ok(
$features
[0]->feature2->end, 217);
ok(
$features
[0]->feature2->strand, -1);
if
(
$DEBUG
) {
my
$out
= new Bio::Tools::GFF(
-gff_version
=> 2);
$out
->write_feature(
$features
[0]);
}