This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Taxonomic analysis of metagenomic data

How can I quickly calculate the taxonomic distribution of a metagenome (assembled or unassembled data)?

sequence
• 2,756 views • link
updated by Shicheng Guo Shicheng Guo • written by bird77 • bird77

4 answers

You can do this very quickly (in a few seconds) with BBMap's Sketch tool, which will compare the data to RefSeq:

sendsketch.sh in=reads.fq reads=4m depth

or

sendsketch.sh in=contigs.fa
• 0 views • link
written by Brian Bushnell Brian Bushnell

There are a number of k-mer based taxonomic tools available: https://omictools.com/taxonomy-dependent2-category

• 0 views • link
written by Sej Modha Sej Modha

What do you mean by Taxonomic analysis? Do you mean taxonomic binning of the sequences or profiling of the sequences?

I wrote a tool in graduate school that uses k-mers and an optimization method for taxonomic profiling of a metagenome dataset in seconds. The tool is named FOCUS and you can learn more about it here. The page also teaches more details about binning and profiling in case you are interested,

• 0 views • link
written by onestop_data onestop_data

How about fastq_screen and Diamond?, I prefer Diamond with my experience.

fastq_screen : https://www.bioinformatics.babraham.ac.uk/projects/fastq_screen/_build/html/index.html

Diamond: https://github.com/bbuchfink/diamond and DIAMOND_analysis_counter.py (SAMSA2)

• 0 views • link
written by Shicheng Guo Shicheng Guo

DIAMOND is simply an aligner. If you don't refer to a tool that does not have a database associated to it, it does no good. Maybe you can add some databases and tools which take the DIAMOND output such as MEGAN which takes the DIAMOND output when aligning against the NR/NT database and can give you back the taxonomic and functional analysis.

I talk here more about DIAMOND and Rapsearch2.

• 0 views • link
written by onestop_data onestop_data

Log in to answer this question.

More posts like this