blastn for selected database and/or organism
1
0
Entering edit mode
aush ▴ 40
@aush-9128
Last seen 2.9 years ago
United States

Is there any way to blast my nucleotide sequence only against specific NCBI database, for example, refseq_rna?

Similarly, I'd like to limit my search to a specific organism/taxid.

https://screenshotscdn.firefoxusercontent.com/images/0d67820a-787a-4d5b-9eb0-1143113f4d20.png

 

As I understand, blastSequences() doesn't have these options now. Is there any other solution?

blast blastsequences • 1.3k views
ADD COMMENT
1
Entering edit mode
@james-w-macdonald-5106
Last seen 1 hour ago
United States

Note that blastSequences is just doing a query to the NCBI blast API, as noted in the help page for the function. Because of that, you can only make the queries that are amenable to the API. You can for example query on a different database (again, as noted in the help page):

>  blastSequences(x = "GGCCTTCATTTACCCAAAATG", database = "reseq_rna")
estimated response time 45 seconds
list()

But there isn't any facility, so far as I can tell, to specify the organism.

ADD COMMENT
0
Entering edit mode

Thank you James! I have contacted NCBI support, and the solution turned out to be quite simple: besides QUERY, their API supports a field ENTREZ_QUERY. Unfortunately it is not included in annotate for now, but I simply added &ENTREZ_QUERY=human%5borgn%5d string to my QUERY, and this worked!

ADD REPLY

Login before adding your answer.

Traffic: 759 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6