Protocol Online logo
Top : Forum Archives: : Protein and Proteomics

Problem with blastpgp - (May/22/2008 )

Hi,

When I run blastpgp.exe on 1a1v.fasta (protein), I get the following errors:

[command prompte>> blastpgp -b 0 -j 3 -h 0.001 -d %dbname% -i 1a1v.fasta -C 1a1v.blast ]

############################################################################
Sequences not found previously or not previously below threshold:
Database: db\nr
Posted date: May 22, 2008 9:59 AM
Number of letters in database: 999,999,866
Number of sequences in database: 2,979,045

Lambda K H
0.313 0.173 0.512
Lambda K H
0.267 0.0527 0.140

Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 3,498,838,404
Number of Sequences: 2979045
Number of extensions: 186639301
Number of successful extensions: 619258
Number of sequences better than 10.0: 500
Number of HSP's better than 10.0 without gapping: 5090
Number of HSP's successfully gapped in prelim test: 6339
Number of HSP's that attempted gapping in prelim test: 586519
Number of HSP's gapped (non-prelim): 16709
length of query: 476
length of database: 999,999,866
effective HSP length: 132
effective length of query: 344
effective length of database: 606,765,926
effective search space: 208727478544
effective search space used: 208727478544
T: 11
A: 40
X1: 16 ( 7.2 bits)
X2: 38 (14.6 bits)
X3: 64 (24.7 bits)
S1: 41 (21.0 bits)
S2: 78 (34.3 bits)
BLASTP 2.2.9 [May-01-2004]

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= 1A1V:B|PDBID|CHAIN|SEQUENCE
(8 letters)
Database: db\nr
2,979,045 sequences; 999,999,866 total letters

***** No hits found ******

***** No hits found ******
Database: db\nr
Posted date: May 22, 2008 9:59 AM
Number of letters in database: 999,999,866
Number of sequences in database: 2,979,045

Lambda K H
-1.00 -1.00 -1.00
Lambda K H
0.267 0.0410 0.140

############################################################################


Whould you please help me with solving it?

Thanks in advanced.

-nhabibi-

Are you sure, that your fasta file hast the correct format? -> Its content should start with >protein name
Did you try to paste the sequence directly and not to upload the sequence file?

-jazim-