This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Error using BLAST+

I am trying to use BLASTx to query a nucleotide sequence against a protein database and I wanted to add the encoded protein information to the output dataset. I however keep getting this error

blastx -db UniprotConus -query /Users/evanclark/Desktop/Cone_Snail_Assembly_split/aaaa-Cone_Snail_Assembly.txt -outfmt qseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore sseq qseq 6 -out result1.txt
Error: Too many positional arguments (1), the offending value: sseqid
Error:  (CArgException::eSynopsis) Too many positional arguments (1), the offending value: sseqid
fasta uniprot blast protein
• 5,916 views • link
updated by Ram Ram • written by onspotproductions onspotproductions

2 answers

You may want to put " or ' around your formatting arguments:

blastx -db UniprotConus -query /Users/evanclark/Desktop/Cone_Snail_Assembly_split/aaaa-Cone_Snail_Assembly.txt -outfmt "qseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore sseq qseq 6" -out result1.txt

If that doesn't work, put the 6 before the other arguments:

blastx -db UniprotConus -query /Users/evanclark/Desktop/Cone_Snail_Assembly_split/aaaa-Cone_Snail_Assembly.txt -outfmt "6 qseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore sseq qseq" -out result1.txt
• 0 views • link
updated by Ram Ram • written by Philipp Bayer Philipp Bayer

more concise with:

blastx -db UniprotConus -query /Users/evanclark/Desktop/Cone_Snail_Assembly_split/aaaa-Cone_Snail_Assembly.txt -outfmt "6 std sseq qseq" -out result1.txt
• 0 views • link
written by Yanpeng Chen • Yanpeng Chen

Log in to answer this question.

More posts like this