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
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