[Bioperl-l] BPlite::HSP Error messages

Peter Kos kos@rite.or.jp" <kos@rite.or.jp
Mon, 1 Apr 2002 16:41:13 +0900


Hi all,

I have BlastX outputs from local NCBI BLAST. Now I just try to count 
them like

while(my $hsp = $sbjct->nextHSP) {
	if (($hsp->bits >= $min_score) && ($hsp->P <= $e_value)){
	$hits{$seqid}++;
	}

every time I touch (perhaps) the nextHSP, I get some error messages 
like

-------------------- WARNING ---------------------
MSG: Frame (-1) did not match strand of query match (1)
---------------------------------------------------

but the whole thing works well of course. I do not think it is a big 
problem, if my query does not match in all six frames. As I figured 
out, it comes from ..../Bio/Tools/BPlite/HSP.pm

As I have thousands of these files, my intended messages are just 
wiped away by the tide of these error messages. I could not SOLVE it, 
but since it is warn, rather than print, I circumvent it by

>( ~/scripts/blstparser.pl > /dev/tty ) | & cat > /dev/null

as I do not need these warnings at all. However, it is usually not 
very clever to neglect all possible error messages, so I am not 
really delighted by this solution. I can use grep before it is 
directed to /dev/null, I know. Is there any intelligent bioperl-like 
way to solve it?

Oh, yes, this is the good old bioperl 0.7 or something. I admit that 
I have not yet been brave enough to attempt the installation of the 
1.0 till so many questions are open about Solaris and firewall and 
other things that cause so much difficulties to me.

Thanks
Peter