Abstract:We investigate the fundamental limits of the recently proposed random access coverage depth problem for DNA data storage. Under this paradigm, it is assumed that the user information consists of $k$ information strands, which are encoded into $n$ strands via a generator matrix $G$. During the sequencing process, the strands are read uniformly at random, as each strand is available in a large number of copies. In this context, the random access coverage depth problem refers to the expected number of reads (i.e., sequenced strands) required to decode a specific information strand requested by the user. This problem heavily depends on the generator matrix $G$, and besides computing the expectation for different choices of $G$, the goal is to construct matrices that minimize the maximum expectation over all possible requested information strands, denoted by $T_{\max}(G)$. In this paper, we introduce new techniques to investigate the random access coverage depth problem, capturing its combinatorial nature and identifying the structural properties of generator matrices that are advantageous. We establish two general formulas to determine $T_{\max}(G)$ for arbitrary generator matrices. The first formula depends on the linear dependencies between columns of $G$, whereas the second formula takes into account recovery sets and their intersection structure. We also introduce the concept of recovery balanced codes and provide three sufficient conditions for a code to be recovery balanced. These conditions can be used to compute $T_{\max}(G)$ for various families of codes, such as MDS, simplex, Hamming, and binary Reed-Muller codes. Additionally, we study the performance of modified systematic MDS and simplex matrices, showing that the best results for $T_{\max}(G)$ are achieved with a specific combination of encoded strands and replication of the information strands.

On the Coverage Required for Diploid Genome Assembly

Hypergraph covering problems motivated by genome assembly questions

Theoretical Bounds on Mate-Pair Information for Accurate Genome Assembly

Reconstruction of Diploid High-Order 3D Genome Interactions from Long Noisy Concatemers

Sequencing depth and coverage: key considerations in genomic analyses

Exploring Genome Characteristics and Sequence Quality Without a Reference

A Combinatorial Perspective on Random Access Efficiency for DNA Storage

Computational Genome Analysis

Covering All Bases: The Next Inning in DNA Sequencing Efficiency

GCI: a continuity inspector for complete genome assembly

Maximum likelihood genome assembly

Sequence coverage required for accurate genotyping by sequencing in polyploid species

Exemplar or Matching: Modeling DCJ Problems with Unequal Content Genome Data

Complex genome assembly based on long-read sequencing

Do Read Errors Matter for Genome Assembly?

Scalable telomere-to-telomere assembly for diploid and polyploid genomes with double graph

New algorithms for accurate and efficient de novo genome assembly from long DNA sequencing reads

DEGAP: Dynamic elongation of a genome assembly path

Comparison of Hi-C-Based Scaffolding Tools on Plant Genomes

Chromosome-scale, haplotype-resolved assembly of human genomes

ComHapDet: A Spatial Community Detection Algorithm for Haplotype Assembly