Question: How to blast mutiple sequences to local database?
gravatar for tmms
2.5 years ago by
tmms0 wrote:

Dear all

I am fairly new to using blast, especially at the command line. In the nex paragraph, I will sketch my problem.

I have two files with sequences (formated as a fasta-file: >header \n sequence ...). I already translated them into proteins. I already made a local database of one of the files with the "makeblastdb" command. Now I would like to blast (blastp) every sequence from the other file to that database, but I have no idea how to do this. In the manual (Blast+) I could not find an answer.

Thanks in advance.

ADD COMMENTlink modified 2.5 years ago by • written 2.5 years ago by tmms0
gravatar for
2.5 years ago by
Philadelphia, PA wrote:
/path/to/blastp -query /path/to/inputfile.fasta -db /path/to/dbfile.fasta -num_threads N -outfmt 6 -out blastp.out.txt

where: -query is your fasta that you're blasting -db is the name of the database from makeblastdb -num_threads is the number of processors to use -outfmt 6 is tab-delimited output (my preference for parsing later) See here for other options. -out is the name of the output file

Use /path/to/blastp -h for more options.

ADD COMMENTlink written 2.5 years ago by
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1702 users visited in the last hour