English

Finding Maximal Exact Matches in Graphs

Data Structures and Algorithms 2023-07-04 v2

Abstract

We study the problem of finding maximal exact matches (MEMs) between a query string QQ and a labeled graph GG. MEMs are an important class of seeds, often used in seed-chain-extend type of practical alignment methods because of their strong connections to classical metrics. A principled way to speed up chaining is to limit the number of MEMs by considering only MEMs of length at least κ\kappa (κ\kappa-MEMs). However, on arbitrary input graphs, the problem of finding MEMs cannot be solved in truly sub-quadratic time under SETH (Equi et al., ICALP 2019) even on acyclic graphs. In this paper we show an O(nLdL1+m+Mκ,L)O(n\cdot L \cdot d^{L-1} + m + M_{\kappa,L})-time algorithm finding all κ\kappa-MEMs between QQ and GG spanning exactly LL nodes in GG, where nn is the total length of node labels, dd is the maximum degree of a node in GG, m=Qm = |Q|, and Mκ,LM_{\kappa,L} is the number of output MEMs. We use this algorithm to develop a κ\kappa-MEM finding solution on indexable Elastic Founder Graphs (Equi et al., Algorithmica 2022) running in time O(nH2+m+Mκ)O(nH^2 + m + M_\kappa), where HH is the maximum number of nodes in a block, and MκM_\kappa is the total number of κ\kappa-MEMs. Our results generalize to the analysis of multiple query strings (MEMs between GG and any of the strings). Additionally, we provide some preliminary experimental results showing that the number of graph MEMs is an order of magnitude smaller than the number of string MEMs of the corresponding concatenated collection.

Keywords

Cite

@article{arxiv.2305.09752,
  title  = {Finding Maximal Exact Matches in Graphs},
  author = {Nicola Rizzo and Manuel Cáceres and Veli Mäkinen},
  journal= {arXiv preprint arXiv:2305.09752},
  year   = {2023}
}

Comments

21 pages, 3 figures. To be published in the proceedings of WABI 2023. This article supersedes part of arXiv:2302.01748

R2 v1 2026-06-28T10:36:22.762Z