Abstract
Shotgun metagenomic sequencing can determine both the taxonomic and functional content of microbiomes. However, functional classification for metagenomic reads remains highly challenging as protein mapping tools require substantial computational resources and yield ambiguous classifications when short reads map to homologous proteins originating from different bacteria. Here we introduce
k
Mermaid for the purpose of uniquely mapping bacterial short reads to taxa-agnostic clusters of homologous proteins, which can then be used for downstream analysis tasks such as read quantification and pathway or global functional analysis. Using a nested hash map containing amino acid
k
-mer profiles as a model for protein assignment,
k
Mermaid achieves the sensitivity of popular existing protein mapping tools while remaining highly resource efficient. We evaluate
k
Mermaid on simulated data and data from human fecal samples as well as demonstrate the utility of
k
Mermaid for classifying reads originating from new, unseen proteins.
k
Mermaid allows for highly accurate, unambiguous and ultrafast metagenomic read assignment into protein clusters, with a fixed memory usage, and can easily be employed on a typical computer.