search

bbuchfink / diamond

Accelerated BLAST compatible local sequence aligner.
GNU General Public License v3.0
992 stars 182 forks source link

Command line options #777

Open saassli21 opened 5 months ago

saassli21 commented 5 months ago

I blastp on my data and it's really much faster than I exprected. i used some filters based on the documentation in "Wiki" section ( https://github.com/bbuchfink/diamond/wiki/3.-Command-line-options ) in order to decrease my dataset but I can't do the link between the option and header of my output file :

" ./diamond blastp --min-score 500 --evalue 0.0001 --query-cover 90 --approx-id 70 -d /lab/saassli/blast/db_swissprot/swissprot -q /lab/saassli/blast/input.fa --header -o /lab/saassli/blast/output_swissprot/matches_filter_evalue.tsv "

for example I know that for --evalue I must survey the e-value column in my output file --min-score it's bit-score column in my output file

but for --approx-id and --query-cover I don't know which column I should survey and how they are calculated Thanks, Saassli,

bbuchfink commented 5 months ago

These columns are not reported by default, you can see them by changing the format e.g. -f 6 qcovhsp approx_pident.