The comprehensive characterization of recombinant adeno-associated viral (rAAV) integration frequency and persistence for assessing rAAV vector biosafety in gene therapy is severely limited due to the predominance of episomal rAAV vector genomes maintained in vivo. Introducing rAAV insertional standards (rAIS), we show that linear amplification-mediated (LAM)-PCR and deep sequencing can be used for validated measurement of rAAV integration frequencies. Integration of rAAV2/1 or rAAV2/8, following intramuscular (IM) or regional intravenous (RI) administration of therapeutically relevant vector doses in nine adult non-human primates (NHP), occurs at low frequency between 10(-4) and 10(-5) both in NHP liver and muscle, but with no preference for specific genomic loci. High resolution mapping of inverted terminal repeat (ITR) breakpoints in concatemeric and integrated vector genomes reveals distinct vector recombination hotspots, including large deletions of up to 3 kb. Moreover, retrieval of integrated rAAV genomes indicated approximately threefold increase in liver compared to muscle. This molecular analysis of rAAV persistence in NHP provides a promising basis for a reliable genotoxic risk assessment of rAAV in clinical trials.