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.

MOSAIK: A hash-based algorithm for accurate next-generation sequencing read mapping

MaxSSmap: A GPU program for mapping divergent short reads to genomes with the maximum scoring subsequence

kngMap: Sensitive and Fast Mapping Algorithm for Noisy Long Reads Based on the K-Mer Neighborhood Graph

ReadsMap: a new tool for high precision mapping of DNAseq and RNAseq read sequences

Bitmapper: an Efficient All-Mapper Based on Bit-Vector Computing

Real-time mapping of nanopore raw signals

Minimap and miniasm: fast mapping and de novo assembly for noisy long sequences

pathMap: a path-based mapping tool for long noisy reads with high sensitivity

Too many needles in this haystack: algorithms for the analysis of next generation sequence data

Umap: Use Unique Sequence for Alignment of Short Sequence Reads and SNP Detection

AMAS: Optimizing the Partition and Filtration of Adaptive Seeds to Speed Up Read Mapping

smsMap: mapping single molecule sequencing reads by locating the alignment starting positions

Sap-A Sequence Mapping And Analyzing Program For Long Sequence Reads Alignment And Accurate Variants Discovery

Fast and accurate mapping of Complete Genomics reads

An efficient Burrows-Wheeler transform-based aligner for short read mapping

A Re-Sequencing Tool For High Mismatch-Tolerant Short Read Alignment Based On Burrows-Wheeler Transform

RawHash: Enabling Fast and Accurate Real-Time Analysis of Raw Nanopore Signals for Large Genomes

Transforming Genomes Using MOD Files with Applications.

MECAT: an ultra-fast mapping, error correction and<i>de novo</i>assembly tool for single-molecule sequencing reads

ESKEMAP: exact sketch-based read mapping

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