For decades, in situ hybridization methods have been essential tools for studies of vertebrate development and disease, as they enable qualitative analyses of mRNA expression in an anatomical context. Quantitative mRNA analyses typically sacrifice the anatomy, relying on embryo microdissection, dissociation, cell sorting, and/or homogenization. Here, we eliminate the tradeo↵ between quantitation and anatomical context, using multiplexed in situ hybridization chain reaction (HCR) to perform accurate and precise relative quantitation of mRNA expression with subcellular resolution within whole-mount vertebrate embryos. Gene expression can be queried in two directions: read-out from anatomical space to expression space reveals co-expression relationships in selected regions of the specimen; conversely, read-in from multidimensional expression space to anatomical space reveals those anatomical locations in which selected gene co-expression relationships occur. As we demonstrate by examining gene circuits underlying somitogenesis, quantitative read-out and read-in analyses provide the strengths of flow cytometry expression analyses, but by preserving subcellular anatomical context, they enable iterative bi-directional queries that open a new era for in situ hybridization. Traditional in situ hybridization approaches based on catalytic reporter deposition (CARD) yield high-contrast images of mRNA expression domains within whole-mount vertebrate embryos