R script used for the phyloseq analysis of the Poirier et al., food microbiota dataset and comparative analysis between 16S and gyrB amplicon sequencing.