Following ITS sequencing, reads were demultiplexed using the 2017.10 release of QIIME2 requiring a barcode match. Because the forward and reverse reads do not overlap for this amplicon, DADA2 was used to trim forward and reverse reads to 200 and 100bp respectively. Cutadapt was used to trim primer sequences from the sequence variants identified by DADA2. Sequences that did not contain these primer sequences were discarded. VSEARCH was used to perform de novo chimera detection and removal. The remaining ITS sequence variants were clustered into operational taxonomic units (OTUs) at 97% similarity using Abundant OTU+. The software package ITSx was then used to filter non-ITS and non-fungal OTUs from our data set
