How to blast multiple sequences against one database at once? (command line)
1
0
Entering edit mode
13 months ago
M. ▴ 30

Hi,

I'm very new to the command line and I am trying to use blast. I have a set of sequences in one fasta file that I need to blast against a database that is also is a fasta file. Is there a blast option or can I do a for loop for that?

Thanks in advance!

unix blast • 523 views
ADD COMMENT
2
Entering edit mode
13 months ago
GenoMax 142k

Unless you want to get separate result files for each sequence you can cat seq1.fa seq2.fa ... seqN.fa > total.fa. And then use total.fa as -query total.fa for your blast. If you do want separate result files then use a loop/xargs construct to feed blast with one file at a time.

ADD COMMENT
0
Entering edit mode

I want seperate results, I'm trying to do a loop. Thank you!

ADD REPLY

Login before adding your answer.

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