[Bioperl-l] Hspwriter and fasta
Vilanova,David,LAUSANNE,NRC/BS
david.vilanova@rdls.nestle.com
Fri, 27 Sep 2002 14:49:00 +0200
Hello,
I'm just trying to parse and write a Fasta report using HSPtableWriter but
this is what I'm getting :
0 0 0 0 0.0e+00 0 0 0.00
Here below is an extract of the script. I think it has something to do with
HSPtableWriter
my $writer = new Bio::SearchIO::Writer::HSPTableWriter;
my $out = Bio::SearchIO->new( -writer => $writer,
-file => ">out"
);
my $searchio = new Bio::SearchIO(-file=>"report.FASTA",
-format=>'fasta');
my $writer = Bio::SearchIO::Writer::HSPTableWriter->new(
-columns =>
[qw(
query_name
query_length
hit_name
hit_length
rank
frac_identical_query
expect
)]
);
while ( my $result = $searchio->next_result() )
{
$out->write_result($result);
}
David