5 mL microfuge tube, air dried briefly and suspended in 200 μL TE buffer resulting in a DNA concentration of approximately 1 μg/μL. Sequencing and annotation Random and φ52237-sequence guided φX216 genome fragment clones were constructed by this website restriction digest of purified φX216 genomic DNA with EcoRI, EcoRI + HindIII or AgeI and ligation with EcoRI, EcoRI + HindIII or SmaI digested pUC19 DNA [25], respectively, followed by

transformation of E. coli DH5α or GBE180 [26] using standard transformation protocols [27] and recovery of white colonies on LB plates containing 100 μg/mL ampicillin and 50 μg/mL 5’-bromo-4-chloro-3-indolyl-β-D-galactopyranoside (X-gal). φ52237-sequence-guided PCR amplicons were designed to close gaps and confirm fragment clone borders. Sequencing was accomplished using M13F and M13R primers, as well as φ52237-sequence guided primer walking of fragment clones and PCR amplicons using an ABI 3130xL Genetic Analyzer (Applied Biosystems, Carlsbad, CA) at the Colorado learn more State University Proteomics and Metabolomics Facility. φX216 Illumina sequencing libraries were prepared using the TruSeq DNA Sample Preparation Kit v2, (Illumina, San Diego, CA), following the manufacturer’s instructions. Phage DNA was fragmented to a range of 300–400 bp using a Covaris acoustic shearing device, (Covaris Inc., Woburn, MA) followed by 3′ adenylation and adapter ligation. Ligation products were purified on an agarose gel and the DNA fragments

enriched via PCR. Fragmented Phage DNA was sequenced by high-throughput Protein kinase N1 Illumina parallel sequencing using 100 bp mate-pair Illumina HiSeq 2000 reversible terminator chemistry. The library was run on 15% of a single lane. Reads were trimmed for quality and de novo short-read genome assembly was performed using the Velvet 1.1.05 sequence assembler algorithms with a hash length of

99 and a final graph with 3 nodes and n50 of 37412 nt [28]. Open reading frames were identified with GeneMark gene prediction software using a viral-optimized Heuristic approach [29]. Putative gene identification was conducted by sequence alignment with φ52237 (GenBank:DQ087285.2) [8] and individual open reading frames queried using the NCBI Basic Alignment Search Tool (BLAST). Genome annotation, mapping, sequence alignments, and comparative analyses were conducted using Gene Construction Kit v3.0 and Geneious Pro 5.4.6 bioinformatics software. The annotation map was created using Adobe Illustrator CS5. The final φX216 genome sequence has been deposited in GenBank under accession # JX681814. Acknowledgements Funding was provided by the Defense Threat Reduction Agency grant W81XWH-07-C0061. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Electronic supplementary material Additional file 1: φX216 host range, word document, Host range of φX216. Table of φX216 host range for 72 B. pseudomallei strains and other Burkholderia species.