MetaStrainer: Accurate reconstruction of bacterial strain genotypes from short-read metagenomic samples.

SummaryMetagenomics provides broad insights from microbial communities, but more biological relevant phenotypes are attributed to subtle changes at the strain-level rather than species. Despite development of several tools using different algorithms, resolving individual strains from short-read pair-end sequencing data remains challenging. We developed MetaStrainer, a tool capable of reconstructing strain genotypes from metagenomic data. Compared with existing approaches, MetaStrainer substantially increases genotype accuracy, correctly identifies the number of strains, and accurately estimates their relative abundances. Accuracy of reconstructed genotypes is robust to choice of mapping reference. Availability and implementationMetaStrainer is implemented in Python 3. Source code and instructions are available on GitHub at https://www.github.com/lbobay/MetaStrainer and on Zenodo: https://doi.org/10.5281/zenodo.17872331 Contactljbobay@ncsu.edu Supplementary InformationSupplementary data is available at Bioinformatics online.