[Bioperl-l] Blastxml
Vilanova,David,LAUSANNE,NRC/BS
david.vilanova@rdls.nestle.com
Thu, 30 May 2002 10:14:31 +0200
Hello,
I'm writting a script to parse an xml blast output.
I use standaloneblast to launch the blast but I'm failing to pass the output
of the blast to searchIO
Does anyone know to correct it ?
use Bio::SeqIO;
use Bio::Tools::Run::StandAloneBlast;
use Bio::SearchIO
...
...
my $blast = Bio::Tools::Run::StandAloneBlast->new(
database => 'swissprot',
program => 'blastx',
'_READMETHOD' => 'Blast',
);
#More blast options
$blast->m(7); #create xml output
my $searchin = new Bio::SearchIO(-format => 'blastxml);
while ($seq = $Seq_in->next_seq() )
{
my $searchin = $blast->blastall($seq);
while (my $result = $searchin->next_result())
{....}
...
}
How can I parse the xml blast output ???
Thanks ,
David Vilanova