English

On the number of $k$-mers admitting a given lexicographical minimizer

Data Structures and Algorithms 2024-12-25 v2 Discrete Mathematics Combinatorics

Abstract

The minimizer of a word of size kk (a kk-mer) is defined as its smallest substring of size mm (with mkm\leq k), according to some ordering on mm-mers. minimizers have been used in bioinformatics -- notably -- to partition sequencing datasets, binning together kk-mers that share the same minimizer. It is folklore that using the lexicographical order lead to very unbalanced partitions, resulting in an abundant literature devoted to devising alternative orders for achieving better balanced partitions. To the best of our knowledge, the unbalanced-ness of lexicographical-based minimizer partitions has never been investigated from a theoretical point of view. In this article, we aim to fill this gap and determine, for a given minimizer, how many kk-mers would admit the chosen minimizer -- i.e. what would be the size of the bucket associated to the chosen minimizer in the worst case, where all kk-mers would be seen in the data. We show that this number can be computed in O(km)O(km) space and O(km2)O(km^2) time. We further introduce approximations that can be computed in O(k)O(k) space and O(km)O(km) time. We also show on genomic datasets that the practical number of kk-mers associated to a minimizer are closely correlated to the theoretical expected number. We introduce two conjectures that could help closely approximating the total number of kk-mers sharing a minimizer. We believe that characterising the distribution of the number of kk-mers per minimizer will help devise efficient lexicographic-based minimizer bucketting.

Keywords

Cite

@article{arxiv.2412.17492,
  title  = {On the number of $k$-mers admitting a given lexicographical minimizer},
  author = {Florian Ingels and Camille Marchet and Mikaël Salson},
  journal= {arXiv preprint arXiv:2412.17492},
  year   = {2024}
}