[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