outfmt 6 BLAST column names
2
0
Entering edit mode
23 months ago
Paula ▴ 60

Hi!

I am using the command line BLAST tool to search sequences against a custom database. I ran the following command:

blastp -query . /s3_protein_Sequences_sol_1_10.faa -db ./database/nr -outfmt 6 -num_threads 16 -max_target_seqs 5 -evalue 1e-30 -out s3_protein_Sequequences_taxonomy.txt

Here is the first line of the output table :

NODE_101_length_82643_cov_14.288202_80  / seq_85 /  98.925  / 93    / 1 /   0 / 146 / 238 / 1   93/     1.72E-63 /  191

As you can see, the names of the columns are not included. Do you know what each number means? I have searched in the documentation but I cannot find the information.

Thanks a lot!

outfmt6 blast output • 3.9k views
ADD COMMENT
2
Entering edit mode
23 months ago

yes, use outfmt 7 , that will include header lines.

or, you can have a look in the blast -help where it is explained what each of the default output outfmt 6 columns means.

few other points:

  • use that much num_threads will likely not speed up your analysis much, at that number you are already on the plateau of speed gain (moreover only limited parts of the blast analysis is multithreaded anyway)
  • be careful when using the max_target_seq parameter, it might not return what you think it will (for details google or look around here on biostars)
ADD COMMENT
2
Entering edit mode
23 months ago
GenoMax 141k

https://www.metagenomics.wiki/tools/blast/blastn-output-format-6

https://www.ncbi.nlm.nih.gov/books/NBK279684/table/appendices.T.options_common_to_all_blast/

When not provided, the default value is:

'qseqid sseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore', which is equivalent to the keyword 'std'

ADD COMMENT

Login before adding your answer.

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