This is a test version of Biostars. For the public version, visit https://www.biostars.org.
blastp commandline: against a protein database of a specific organism

Hello,

It is a very basic question regarding the -db option. I am using standalone blast+ and I want to blast my query.fasta ONLY against specific organism (e.g human) protein database only. I know blasting against refseq_protein/ swissprot are the usual choice

I am confused which database to select from ftp://ftp.ncbi.nlm.nih.gov/blast/db/ .

How can I specify the specific organism in command line.

Thanks

blast command-line-blast blastp

I am also search like this, did you get answer for that...

Thanks

2 answers

You might want to check this thread: Vertebrate Subset Nr Database? Build My Own?

There is also information in BLAST Command Line Applications User Manual under the heading "Extract all human sequences from the nr database"

http://www.ncbi.nlm.nih.gov/books/NBK1763/

Thanks Siva, that was very helpful

Another alternative is to download all the peptide sequences using BioMart for your genome of interest.

Choose Ensembl Genes --> Homo Sapiens --> Sequences, then choose peptides, output to FASTA, and save. Then build locally as you would any other blast database. This will obviously miss heavily processed isoforms, so it depends what your application is.

Thanks hurfdurf. I will definitely try this one.

Log in to answer this question.