Question: How to blast mutiple sequences to local database?
0
gravatar for tmms
2.5 years ago by
tmms0
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 st.ph.n2.5k • written 2.5 years ago by tmms0
0
gravatar for st.ph.n
2.5 years ago by
st.ph.n2.5k
Philadelphia, PA
st.ph.n2.5k 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 st.ph.n2.5k
Please log in to add an answer.

Help
Access

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