Protocol Online logo
Top : New Forum Archives (2009-): : Bioinformatics and Biostatistics

Stand Alone Blast - Creating a Blast Database (May/03/2011 )

I basically have these 2 files, the

Input file = sequences.txt
Database file = homologene_result.txt

Isn't "blastall" an older version command ? That is why i was using makeblastdb

This is what i typed

172-30-8-231:bin sukeerthiteja$ makeblastdb -in homologene_result.txt -dbtype prot -parse_seqids -out teja1 blastp -query sequences.txt -db teja2

AND i also tried the code using .txt next to teja2.txt, but in both cases i got this error.

Error: (106.18)
NCBI C++ Exception:
Error: (CArgException::eSynopsis) Too many positional arguments (1), the offending value: blastp
Error: (CArgException::eSynopsis) Application's initialization failed

Please help me understand where exactly i am going wrong.

Thanking you
Teja

-Sukeerthi Teja-

First, you must turn your FASTA formatted sequences into a BLAST database:

formatdb -i input_sequences.fasta -p T -n input_sequences_database_name.db


Then you can BLAST your query sequence(s) against the database you made:

blastall -p blastp -d input_sequences_database_name.db -F F -i query_sequence.fasta -o blast_output_file.txt

-HomeBrew-