Mind the Gaps: Validating species presence under sparse genomic recovery in clinical metagenomics

1Clinical metagenomic sequencing is increasingly used to detect pathogens directly from clinical samples, but sparse and uneven genome recovery caused by low microbial biomass and stochastic fragment loss makes species-level validation difficult. Under these conditions, the same observed coverage may arise either from genuine genome-wide representation or from reads confined to conserved loci shared among related taxa. Current heuristic thresholds based on read counts, breadth or related alignment summaries target the realised coverage rather than the generating process. We show that, under hidden fragment loss, coverage is the wrong inferential object for uniformly valid species-level inference. What remains informative is the spatial dispersion of surviving reads across genomic coordinates. This change of inferential object allows us to ask how much of the candidate genome must be represented in any source capable of reproducing the observed evidence. This yields the Metagenomic Alignment Validation Index (MAVI), a conservative lower bound on the fraction of the candidate genome required to explain the observed read dispersion. MAVI therefore provides an error-controlled, interpretable measure of genomic confidence: the minimum genomic content a competing source would need to share with the candidate to mimic the observed evidence. Evaluated across 643 long-read metagenomic samples, including 280 clinical samples, MAVI reproduced independently supported organism calls, resolved sparse low-coverage ambiguities and exposed workflow-level collapse of stochastic dispersion, without relying on heuristic coverage thresholds.