Virulence and resistance genes carried on mobile DNAs such as genomic islands (GIs) and plasmids promote bacterial pathogen emergence. An early step in the mobilization of GIs is their excision, which produces both a circular form of the GI and a deletion site in the chromosome
circular forms have also been described for some bacterial insertion sequences (ISs). We demonstrate that the recombinant sequence produced at the junction of such circles, and their corresponding deletion sites, can be detected sensitively in high throughput sequencing data, using new computational methods that enable empirical discovery of new mobile DNAs. Applied to the rich mobilome of a single strain (Kpn2146) of the emerging multidrug-resistant pathogen Klebsiella pneumoniae, our approach detected circular junctions for six GIs and seven IS types (several of the latter not previously known to circularize). Our methods further revealed differential biology of multiple mobile DNAs, imprecision of integrases and transposases, and differential activity among identical IS copies for IS26, ISKpn18 and ISKpn21. Exonuclease was used to enrich for circular dsDNA molecules, and internal calibration with the native Kpn2146 plasmids showed that not all molecules bearing GI and IS circular junctions were circular dsDNAs. Transposition events were also detected, revealing replicon preference (ISKpn18 preferring a conjugative IncA/C2 plasmid), local action (IS26), regional preferences, selection (against capsule synthesis), and left-right IS end swapping. Efficient discovery and global characterization of numerous mobile elements per experiment will allow detailed accounting of bacterial evolution, explaining the new gene combinations that arise in emerging pathogens.