[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