[Bioperl-l] BPlite to parse multi query sequence report

nkuipers nkuipers@uvic.ca
Tue Jan 14 18:39:59 EST 2003


This should do the trick.

#parse the BLAST report
	my $parser = Bio::Tools::BPlite->new( -file => $toBeParsed );
	{
	while ( my $subject = $parser->nextSbjct ) {
		while ( my $hsp = $subject->nextHSP ) { #your code here
		}
	}
	last if $parser->_parseHeader == -1;
	redo;
	}

>===== Original Message From "E Klee" <klee_17@hotmail.com> =====
>Ok, I am still struggling w/ the BioPerl Blast parsers.  Can I use BPlite to
>parse a blast output file that contains multiple query sequences compared to
>multiple subject sequences.  My experience with this so far is openning a
>new BPlite report will only reference the output for the first query
>sequence in my file.  Using a while() statement to step through all subjects
>via nextSbjct() is only finding the hits to the first query sequence.  Any
>help?
>
>-E
>
>_________________________________________________________________
>MSN 8: advanced junk mail protection and 2 months FREE*.
>http://join.msn.com/?page=features/junkmail
>
>_______________________________________________
>Bioperl-l mailing list
>Bioperl-l@bioperl.org
>http://bioperl.org/mailman/listinfo/bioperl-l




More information about the Bioperl-l mailing list