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$
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$ – Chris_Rands Aug 12 '20 at 21:04
  • 1
    $\begingroup$ It does work well, but it's less accurate at lower taxonomic ranks (as most tools are) $\endgroup$ – aLbAc Aug 14 '20 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, privacy policy and cookie policy

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