This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Making Specialised A Local Blast Db

Hi,

I need to make a local Blast DNA (nt) database that only includes sequences from plants. I'd also like to make one with only chloroplast sequences. Put I can't find any tools to do this. I don't want to download the entire nt and then use masking. I'd rather be using a smaller db.

Are there tools to do this? update_blastdb.pl does not seem to be able to 'filter' sequences.

I can use Entrez to get a subset nt, but can't see a way to download the output (the download browser buttons don't work for this - the download always freezes).

Thanks,

T.

blast database

Do you know which specific sequences you want in the database (like a list of identifiers)? Do you have fasta files of these sequences?

As Dk says below, please do not post your comments and notes as answers. I have deleted them, since they were detracting from the question.

8 answers

You can use eutils to get the IDs programmatically and download the according Fasta files. However this will take some scripting.

If you have the Fasta file you can use makeblastdb that comes with the local blast installation to make a database out of the Fasta files.

Great - thanks, that seems to work.

T.

Hi,

no, I don't have identifiers. I want to build it by taxonomic subdivision - PLN, or by keywords or other information, since taxon ids are sometimes wrong or absent.

Perhaps I could get a list of identifiers only from Entrez - I have not tried this for large lists.

Thanks.

update: no, there does not seem to be a way to output a list of identifiers from entrez (except one page at a time).

update: Entrez seems to have a 200 item limit on output via the browser

Thanks - I thought there must be a tool like eutils. Looks like I have to learn yet a nother scripting language and set of tools... jeez.

T.

Ok, I have tried to use this script:

http://www.ncbi.nlm.nih.gov/books/NBK25498/#chapter3.Application_3_Retrieving_large

but always get an empty output file and the script runs suspiciously quickly.

Any ideas?

Thanks,

T.

my IT dept. is blocking 'Get' use!! :(

Hi T. Just for future reference, please use the reply button to the comments instead of posting a new answer. We don't want to clutter up the question with a bunch of posts that aren't actually answers.

Log in to answer this question.