BEGIN {
if
( $@ ) {
}
plan
tests
=> 29;
}
ok(1);
my
$verbose
=
$ENV
{
'BIOPERLDEBUG'
} || 0;
my
@formats
=
qw(gcg fasta raw pir tab ace )
;
foreach
my
$format
(
@formats
) {
print
"======== $format ========\n"
if
$verbose
;
read_write(
$format
);
}
sub
read_write {
my
$format
=
shift
;
my
$seq
;
my
$str
= Bio::SeqIO->new(
-file
=> Bio::Root::IO->catfile
(
"t"
,
"data"
,
"test.$format"
),
-format
=>
$format
);
ok
$seq
=
$str
->next_seq();
print
"Sequence 1 of 2 from $format stream:\n"
,
$seq
->seq,
"\n\n"
if
$verbose
;
unless
(
$format
eq
'raw'
) {
ok
$seq
->id,
'roa1_drome'
,
"ID for format $format"
;
ok
$seq
->
length
, 358;
}
unless
(
$format
eq
'gcg'
) {
ok
$seq
=
$str
->next_seq();
print
"Sequence 2 of 2 from $format stream:\n"
,
$seq
->seq,
$seq
->seq,
"\n"
if
$verbose
;
}
my
$out
= Bio::SeqIO->new(
-file
=>
">"
. Bio::Root::IO->catfile
(
"t"
,
"data"
,
"$format.out"
),
-format
=>
$format
);
ok
$out
->write_seq(
$seq
);
if
(
$format
eq
'fasta'
) {
my
$id_type
;
ok(
$id_type
=
$out
->preferred_id_type(
'accession.version'
),
'accession.version'
);
}
}
END {
map
{
unlink
Bio::Root::IO->catfile(
"t"
,
"data"
,
"$_.out"
) }
@formats
}