blastp commandline: against a protein database of a specific organism
2
0
Entering edit mode
6.9 years ago
Anushka ▴ 20

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 blastp db command line blast • 7.4k views
ADD COMMENT
0
Entering edit mode

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

Thanks

ADD REPLY
1
Entering edit mode
6.9 years ago
Siva ★ 1.8k

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/

 

ADD COMMENT
0
Entering edit mode

Thanks Siva, that was very helpful

ADD REPLY
1
Entering edit mode
6.9 years ago
hurfdurf ▴ 490

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.

ADD COMMENT
0
Entering edit mode

Thanks hurfdurf. I will definitely try this one.

ADD REPLY

Login before adding your answer.

Traffic: 1604 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