To analyze pyrosequencing data, raw sequences were first screened to remove redundancy using Galaxy software (http:// usegalaxy.org). The remaining reads were then imported into CLC Genomics Workbench (CLC bio, Aarhus, Denmark), trimmed to remove Nextera-specific transposon sequences as well as short and low quality reads, and assembled using the CLC de novo assembler. Both singleton and assembled contiguous sequences (contigs) were queried against the GenBank database (http:// www.ncbi.nlm.nih.gov/GenBank) using the basic local alignment search tools blastn and blastx [12] , with an e-value cut-off of 10 and word sizes of 11 and 3 for blastn and blastx queries, respectively.