Near-Optimal Average-Case Approximate Trace Reconstruction from Few Traces

Abstract: In the standard trace reconstruction problem, the goal is to \emph{exactly} reconstruct an unknown source string $\mathsf{x} \in {0,1}<sup>n$ from independent "traces", which are copies of $\mathsf{x}$ that have been corrupted by a $\delta$-deletion channel which independently deletes each bit of $\mathsf{x}$ with probability $\delta$ and concatenates the surviving bits. We study the \emph{approximate} trace reconstruction problem, in which the goal is only to obtain a high-accuracy approximation of $\mathsf{x}$ rather than an exact reconstruction. We give an efficient algorithm, and a near-matching lower bound, for approximate reconstruction of a random source string $\mathsf{x} \in {0,1}<sup>n$ from few traces. Our main algorithmic result is a polynomial-time algorithm with the following property: for any deletion rate $0 &lt; \delta &lt; 1$ (which may depend on $n$), for almost every source string $\mathsf{x} \in {0,1}<sup>n$, given any number $M \leq \Theta(1/\delta)$ of traces from $\mathrm{Del}<em>\delta(\mathsf{x})$, the algorithm constructs a hypothesis string $\widehat{\mathsf{x}}$ that has edit distance at most $n \cdot (\delta M)<sup>{\Omega(M)}$ from $\mathsf{x}$. We also prove a near-matching information-theoretic lower bound showing that given $M \leq \Theta(1/\delta)$ traces from $\mathrm{Del}</em>\delta(\mathsf{x})$ for a random $n$-bit string $\mathsf{x}$, the smallest possible expected edit distance that any algorithm can achieve, regardless of its running time, is $n \cdot (\delta M)<sup>{O(M)}$.