diamond icon indicating copy to clipboard operation
diamond copied to clipboard

i have question for diamond blastp

Open gyudae opened this issue 3 years ago • 1 comments

hi, Im student who studying bioinformatics.

i tried to annotate my metagenome data using diamond tools, and i want to get only one best alignment result for each CDS reads. I found the --max-target-seqs option, but i worry whether i can set this parameter to 1. Because they said in the explanation "Note that this parameter does not only affect the reporting, but also the algorithm as it is taken into account for heuristics that eliminate hits prior to full gapped extension."

Please tell me how can i get only one best alignment without affecting the algorithm. Or if i dont have to worry about the note, then please explain about that

thanks for reading my question

GyuDae Lee

gyudae avatar Nov 30 '22 03:11 gyudae

You can use --max-target-seqs 1 without worry normally, the effect on the accuracy is small. To not affect the algorithm you can add the --no-ranking option.

bbuchfink avatar Dec 25 '22 10:12 bbuchfink