2
$\begingroup$

I have a couple of hundred bacterial sequences of 2-30 genes of interest each, recovered from metagenomics. None of them encode rRNA. Normally I'd just BLAST the one gene I already know to be reliable phylogenetic marker, where present, but the scale of my current dataset makes it a bit daunting.

Is there any software for taxonomic assignment that you could recommend? It's fine if it requires command line use. There are so many options listed online that I didn't even know where to start... Please and thank you.

$\endgroup$

2 Answers 2

4
$\begingroup$

You could take a look at kraken2 It should fit your purpose.

$\endgroup$
2
  • 1
    $\begingroup$ yes, kraken2 performed well in an independant benchmark lemmi.ezlab.org $\endgroup$ Commented Aug 12, 2020 at 21:04
  • 1
    $\begingroup$ It does work well, but it's less accurate at lower taxonomic ranks (as most tools are) $\endgroup$
    – aLbAc
    Commented Aug 14, 2020 at 9:52
3
$\begingroup$

If you like blast, maybe you can try DIAMOND - the authors claim it's much faster than BLAST.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.