Abstract:A key step in sequence similarity search is to identify seeds that are found in both the query and the reference sequence. A seed is a shorter substring (e.g., a k-mer) or pattern (e.g., a spaced k-mer) constructed from the sequences. A well-known trade-off in applications such as read mapping is that longer seeds offer fast searches through fewer spurious matches but lower sensitivity in variable regions as longer seeds are more likely to harbor mutations. Some recent developments on seed constructs have considered approximate (or fuzzy) seeds such as k-min-mers, strobemers, BLEND, SubSeqHash, TensorSketch, and more, that can match over smaller mutations and, thus, suffer less from sensitivity issues in variable regions. Nevertheless, the sensitivity-to-speed trade-off still exists for such constructs. In other applications, such as genome assembly, using multiple sizes of k-mers is effective. While this can be achieved in read mapping through, e.g., MEM construction from an FM-index, such seed constructs are typically much slower than hash-based constructs. To this end, we introduce multi-context seeds (MCS). In brief, MCS are strobemers where the hashes of individual strobes are partitioned in the hash value representing the seed. Such partitioning enables a cache-friendly approach to search for both full and partial matches of a subset of strobes. For example, both the full strobemer and the first strobe (a k-mer) can be queried. We demonstrate that MCS improves sequence matching statistics over standard strobemers and k-mers without compromising seed uniqueness. We demonstrate the practical applicability of MCS by implementing them in strobealign. Strobealign with MCS comes at no cost in memory and only little cost in runtime while offering increased mapping accuracy over default strobealign using simulated Illumina reads across genomes of various complexity. We also show that strobealign with MCS outperforms minimap2 in short-read mapping and is comparable to BWA-MEM in accuracy in high-variability sequences. MCS provides a fast seed alternative that addresses the trade-offs between seed length and alignment accuracy.

Seeding with Minimized Subsequence.

Efficient Seeding for Error-Prone Sequences with SubseqHash2

Efficient Algorithms for Finding a Longest Common Increasing Subsequence

MISSH: Fast Hashing of Multiple Spaced Seeds

CHTKC: a robust and efficient k-mer counting algorithm based on a lock-free chaining hash table

A Novel Key Point Based MLCS Algorithm for Big Sequences Mining

Mining Multiple Longest Common Subsequence in Big Sequences

Short Read Alignment Based On Maximal Approximate Match Seeds

Efficient Approximate Subsequence Matching Using Hybrid Signatures

Towards a Better Solution to the Shortest Common Supersequence Problem: the Deposition and Reduction Algorithm

Multi-context seeds enable fast and high-accuracy read mapping

A compressive seeding algorithm in conjunction with reordering-based compression

An Efficient Mining Algorithm for Key Segment from DNA Sequences

A simple refined DNA minimizer operator enables 2-fold faster computation

A Space-Saving Based MLCS Algorithm

A Branch Elimination-Based Efficient Algorithm for Large-Scale Multiple Longest Common Subsequence Problem

Overlap Digraph： an Effective Model for Finding Good Spaced Seeds for Biological Sequence Local Alignment

Distortion-Resistant Hashing for rapid search of similar DNA subsequence

An Efficient Filtration Method Based on Variable-Length Seeds for Sequence Alignment.

Fundamental Limits of Multiple Sequence Reconstruction from Substrings

Maximum Match Subsequence Alignment Algorithm Finely Grained (MMSAA FG)