[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