16S rRNA gene variable region amplicon data analysis
Reads were generated from specimens according to the routine protocol of the University of Michigan Microbiome Core. They were amplified with Earth Microbiome Project primers targeting the V4 variable region of the gene (intended to generate full overlap in the paired-end reads) with an error-correcting polymerase, followed by sequencing with the Illumina MiSeq platform. Processing of the raw FASTQ files was as per Golob et al.