mpirun -np 22 mpiblast -e 0.02 -m 7 -v 20 -b 20 -I -p blastp -d nr -i prova.fas -o prova.xml
where prova.fas is this file:
>gi|23613525|ref|NP_704546.1| rifin [Plasmodium falciparum 3D7]
MKVHYINILLFALPLNILIYNQRNHNSTTHHTLKIPITRLLCECELYELANYDNDPEMKEVMQQFEDRTTQRFHEYDERM
KTTRQKCKDKCDKEIQKIILKDKLEKELMDKFATLHTDIQSDAIPTCICEKSLEDKMEKECLKCAQNLGGIVAPSTGVLG
EIAALAVNAWKTEAIAAATKAAIAKGTAKGLAAGAAKGVAEVIAQVESQFRLSTIGVKELGSIFNASNYTNETFISGYIY
AQYQGSQCGSLSMLLGKSKPFCTFVEGRIFATSVRVGRSFSPEDFIKTTVQTIVKNAKTTAEATKAQVASAEKAAVLETS
KKAIEATTTPYYTPIIVSIVAIVVIILIMVIIYKILRYRRKKKMKKKLQYIKLLKE
If I try to do the same with the blastall executable, I obtain the requested 20 hits, and nothing more.
Thank you Mike for help me. Please try to execute the same file with the same parameters than I am trying and comment me how are your results.
Thank you very much.
2006/5/30, Mike Cariaso <[EMAIL PROTECTED]>:
I use -b and -v successfully with mpiblast. Are you perhaps counting the number of HSPs instead of the number of hits?
----- Original Message ----
From: Gabri Aparicio <[EMAIL PROTECTED] >
To: [email protected]
Sent: Tuesday, May 30, 2006 10:43:54 AM
Subject: [Mpiblast-users] NUMBER OF OUTPUT HITS DOES NOT MATCH WITH THE INPUT PARAMETERS
Hi,
I want to obtain only a certain number of hits in the output XML. In the sequential BLASTALL I can solve it by using the -v and -b parameters, but the same line in MPIBLAST does not work properly.
I just want to obtain 20 blast hits and MPIBLAST always give me 274 blast hits. It seems MPIBLAST does not analize -v and -b parameters.
Do you have any idea of how to solve my problem?
Thanks.
