English

Sketching and Sequence Alignment: A Rate-Distortion Perspective

Information Theory 2021-07-12 v1 math.IT

Abstract

Pairwise alignment of DNA sequencing data is a ubiquitous task in bioinformatics and typically represents a heavy computational burden. A standard approach to speed up this task is to compute "sketches" of the DNA reads (typically via hashing-based techniques) that allow the efficient computation of pairwise alignment scores. We propose a rate-distortion framework to study the problem of computing sketches that achieve the optimal tradeoff between sketch size and alignment estimation distortion. We consider the simple setting of i.i.d. error-free sources of length nn and introduce a new sketching algorithm called "locational hashing." While standard approaches in the literature based on min-hashes require B=(1/D)O(logn)B = (1/D) \cdot O\left( \log n \right) bits to achieve a distortion DD, our proposed approach only requires B=log2(1/D)O(1)B = \log^2(1/D) \cdot O(1) bits. This can lead to significant computational savings in pairwise alignment estimation.

Keywords

Cite

@article{arxiv.2107.04202,
  title  = {Sketching and Sequence Alignment: A Rate-Distortion Perspective},
  author = {Ilan Shomorony and Govinda M. Kamath},
  journal= {arXiv preprint arXiv:2107.04202},
  year   = {2021}
}
R2 v1 2026-06-24T04:01:43.114Z