Quantcast
Viewing all articles
Browse latest Browse all 41826

Blast Gives Cryptic Errors

I have a list of proteins in fasta format (say goodProteins.fasta). I first make a compatible database using NCBI's formatdb v 2.2.18:

formatdb -i goodProteins.fasta -p T -o T

This gives me a number of files (.psq, .pin, .phr, .psi, .psd).

Then, I run NCBI's blastall :

blastall -F 'm S' -v 100000 -b 100000 -z 67924 -e 1e-5  -p blastp -d goodProteins.fasta -m 8 -o blast.out -i goodProteins.fasta

This starts the blast but then proceeds to give me errors such as

[NULL_Caption] ERROR: SeqPortNew: lcl|5579_goodProteins.fasta start(384) >= len(363)
[NULL_Caption] ERROR: SeqPortNew: lcl|5579_goodProteins.fasta start(404) >= len(363)
[NULL_Caption] ERROR: SeqPortNew: lcl|5579_goodProteins.fasta start(433) >= len(363)
[NULL_Caption] ERROR: SeqPortNew: lcl|5579_goodProteins.fasta start(453) >= len(363)
[NULL_Caption] ERROR: SeqPortNew: lcl|5579_goodProteins.fasta start(472) >= len(363)
[NULL_Caption] ERROR: SeqPortNew: lcl|5579_goodProteins.fasta start(491) >= len(363)

I have uploaded all the data files at http://www.filedropper.com/datatar so anyone can reproduce the errors.

I found this answer to a related question, but I don't see any spaces in my id's.


Viewing all articles
Browse latest Browse all 41826

Trending Articles