[Bioperl-l] no results with standalone tblastn

Torsten Seemann torsten.seemann at infotech.monash.edu.au
Mon Feb 6 21:17:40 EST 2006


> I have downloaded all the databases, as you recommended me. And it is 
> working, but I don't get any results, if I try it online it works fine.
> my result file looks like that:
> 
> TBLASTN 2.2.13 [Nov-27-2005]
> Reference: Altschul, Stephen F., Thomas L. Madden, Alejandro A. Schaffer,
> Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997),
> "Gapped BLAST and PSI-BLAST: a new generation of protein database search
> programs",  Nucleic Acids Res. 25:3389-3402.
> Query=
>          (8 letters)
> Database: All GenBank+EMBL+DDBJ+PDB sequences (but no EST, STS,
> GSS,environmental samples or phase 0, 1 or 2 HTGS sequences)
>            3,749,503 sequences; 16,556,997,203 total letters
> Searching..................................................done
> Sequences producing significant alignments:                Score    
> E      (bits) Value

Is your query only 8 amino acids long?

This report looks like it did have alignments that were not displayed, 
otherwise it would print "**** No hits ****".

This mailing list is not here to solve your BLAST problems unless it is 
a problem with the Perl module running BLAST.

You first need to try and get your problem working on the command line 
*without* Perl. eg.

/home/Hubert/blast/blast-2.2.13/bin/blastall -p tblastn -d nt -i 
YOUR_FASTA_FILE_WITH_SEQUENCE_IN_IT -o OUTPUT_FILE.txt -e 0.001
...

where "..." is the rest of the options you are setting in your Perl 
script. If it doesn't work that way, it will never work in Perl.

-- 
Torsten Seemann
Victorian Bioinformatics Consortium, Monash University, Australia
http://www.vicbioinformatics.com/


More information about the Bioperl-l mailing list