Species-Specific Output From Blast -Db Nt -Remote ...

Question: Species-Specific Output From Blast -Db Nt -Remote ...

8.6 years ago by

Mayo Clinic College of Medicine

kevin.l.neff • 280 wrote:

I'm using the C++ version of blastn, with -db nt and -remote. I need to limit my search to a specific organism. I've had success using -gilist <file> to do this on searches with a local database, but that option is excluded for remote searches. Should I use a script to filter the search results, like -gilist, or is there an easier way?

blast • 5.7k views

ADD COMMENT • link •

modified 6.4 years ago by chris.mit7 • 60 • written 8.6 years ago by kevin.l.neff • 280

I would also like to know this, but for local blast vs nr!

ADD REPLY • link written 8.6 years ago by Lee Katz • 3.1k

6.4 years ago by

chris.mit7 • 60

United States

chris.mit7 • 60 wrote:

Add -entrez_query "Homo Sapiens[Organism]" to your search or whatever organism you want.

Such as:

blastn -remote -db nr -query stuff_to_search.fasta -entrez_query "Homo Sapiens[Organism]" -evalue 1e-20 -num_alignments 10 > stuff_to_search.blastn

ADD COMMENT • link written 6.4 years ago by chris.mit7 • 60

how about if more than one species.

thanks

ADD REPLY • link written 5.6 years ago by Manoj • 80

8.6 years ago by

Paulo Nuin • 3.7k

Canada

Paulo Nuin • 3.7k wrote:

Use the -gilist option in the command line. Create a file with the GIs you want, for instance all bacterias, all mammals, one frog, one lizard, etc and call that file using the option

... -gilist prots.txt ...

ADD COMMENT • link written 8.6 years ago by Paulo Nuin • 3.7k

Sadly, this doesn't work with -remote option.

ADD REPLY • link written 8.6 years ago by kevin.l.neff • 280

Yes, you are right. IMHO, you would be best to download the nr DB and run locally.

ADD REPLY • link written 8.6 years ago by Paulo Nuin • 3.7k

How would you generate this gi list for an entire genus or other taxonomy?

ADD REPLY • link written 8.6 years ago by Lee Katz • 3.1k

Best option would be to go to NCBI website and browse via Taxonomy. This way it is possible to go down/up the tree and check the GIs, as mostly there's one set be genus, species, etc.

ADD REPLY • link written 8.6 years ago by Paulo Nuin • 3.7k

Yes, I agree. However, my web server only allows 10 GB data, so I need to come up with an alternative. Dreamhost or some other service may be the easiest solution, but I wondered if there was a way to do this without bringing in the big guns. Seems like no...

ADD REPLY • link written 8.6 years ago by kevin.l.neff • 280

http://www.biostars.org/post/show/15002/make-a-custom-blast-library-using-the-output-of-another-blast-result/#45602

ADD REPLY • link written 8.6 years ago by kevin.l.neff • 280

Please log in to add an answer.

Similar posts • Search »