Abstract
Seed-chain-extend withk-mer seeds is a powerful heuristic technique for sequence alignment used by modern sequence aligners. Although effective in practice for both runtime and accuracy, theoretical guarantees on the resulting alignment do not exist for seed-chain-extend. In this work, we give the first rigorous bounds for the efficacy of seed-chain-extend withk-mersin expectation. Assume we are given a random nucleotide sequence of length ∼nthat is indexed (or seeded) and a mutated substring of length ∼m≤nwith mutation rate θ < 0.206. We prove that we can find ak= Θ(logn) for thek-mer size such that the expected runtime of seed-chain-extend under optimal linear-gap cost chaining and quadratic time gap extension isO(mnf(θ)logn), wheref(θ) < 2.43 · θ holds as a loose bound. The alignment also turns out to be good; we prove that more than1−O(1m)fraction of the homologous bases isrecoverableunder an optimal chain. We also show that our bounds work whenk-mers aresketched, that is, only a subset of allk-mers is selected, and that sketching reduces chaining time without increasing alignment time or decreasing accuracy too much, justifying the effectiveness of sketching as a practical speedup in sequence alignment. We verify our results in simulation and on real noisy long-read data and show that our theoretical runtimes can predict real runtimes accurately. We conjecture that our bounds can be improved further, and in particular,f(θ) can be further reduced.
Funder
Natural Sciences and Engineering Research Council of Canada
Publisher
Cold Spring Harbor Laboratory
Subject
Genetics (clinical),Genetics
Cited by
2 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献