Question: Blast Multifasta Against Ncbi Nr Database, And Download All The Results
0
gravatar for cristianrohr768
6.6 years ago by
Spain
cristianrohr76830 wrote:

Hello, i have a fasta file with 7600 nucleotide sequences, i have a local nr database formated for blast, with makeblastdb, i don't use formatdb because it give errors (zero length sequences).

I blast all my 7600 sequences, but only 20! have a hit against the nr database.

This looks very strange to me, so i pick some sequences, and blast manually from the NCBI website, and obviously this give to me results for every sequence..

so, 2 questions, what could be wrong in my local blast?

and two, i can blast all my sequences using the NCBI webpage, but i can't figure out how to download all the results at once, instead i have to download manually one by one, and i can't do this for 7600 sequences, any workaround?

thanks in advance

Cristian

ncbi blast download • 4.3k views
ADD COMMENTlink modified 6.6 years ago by Neilfws48k • written 6.6 years ago by cristianrohr76830

What's the blast command you used to blast against your local NR database?

ADD REPLYlink written 6.6 years ago by Damian Kao15k

blastall -p blastn -d nr -i ps.fas -o blast_out.txt

ADD REPLYlink modified 6.6 years ago by Alex Reynolds29k • written 6.6 years ago by cristianrohr76830
2
gravatar for Neilfws
6.6 years ago by
Neilfws48k
Sydney, Australia
Neilfws48k wrote:

I can see two potential problems.

First, you have formatted the BLAST database using makeblastdb, which is from the newer BLAST+ software package. However, to run BLAST you are using the command blastall, which is from the older version of BLAST. You should be running blastn if you used makeblastdb, or blastall if formatdb was used. Messages about zero-length sequences are not necessarily errors; they may just be warnings.

Second: this is probably not an issue but NCBI uses nr to denote non-redundant protein, nt for non-redundant nucleotide. I just wanted to check that your nr is indeed nucleotide sequences.

ADD COMMENTlink written 6.6 years ago by Neilfws48k

+1 In addition to these points there may also be some formatting issues with the input, given the error messages. I think that formatting the db correctly and using the correct db for the program will probably solve this problem though.

ADD REPLYlink written 6.6 years ago by SES8.2k
1

The problem is to me obviously the blastn vs nr. Try "blastx" instead.

ADD REPLYlink written 6.6 years ago by Manu Prestat3.9k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2013 users visited in the last hour