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

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

2.7 years ago by

Mayo Clinic College of Medicine

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

ADD COMMENT • link •

modified 5 months ago by chris.mit7 • 40 • written 2.7 years ago by kevin.l.neff • 60

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

ADD REPLY • link written 2.7 years ago by Lee Katz • 2.5k

2.7 years ago by

Paulo Nuin • 3.5k

Canada

Paulo Nuin • 3.5k 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 2.7 years ago by Paulo Nuin • 3.5k

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

ADD REPLY • link written 2.7 years ago by kevin.l.neff • 60

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

ADD REPLY • link written 2.7 years ago by Paulo Nuin • 3.5k

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

ADD REPLY • link written 2.7 years ago by Lee Katz • 2.5k

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 2.7 years ago by Paulo Nuin • 3.5k

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 2.7 years ago by kevin.l.neff • 60

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

ADD REPLY • link written 2.7 years ago by kevin.l.neff • 60

5 months ago by

chris.mit7 • 40

United States

chris.mit7 • 40 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 5 months ago by chris.mit7 • 40

Please log in to add an answer.

Similar posts • Search »