VCF collection for phylogenies and comparative genomics

In order to compare our novel isolates with previously described members of various S. cerevisiae clades, including animal husbandry isolates and a Hungarian baker’s yeast, raw sequencing files from literature [Offei, Vandecruys, De Graeve, Foulquié-Moreno, and Thevelein 2019; Peter et al. 2018; Duan, Han, Wang, Liu, Shi, Li, Zhang, and Bai 2018; Imre et al. 2022; Rácz et al. 2021] were downloaded and included in our pipeline. The complete list of genomes used from literature is in Table S1. Using BAM files, local realignment around indels and joint variant calling and filtering for the isolates were performed with GATK 4.1.9.0. [Poplin et al. 2018; Van der Auwera et al. 2013] with regions annotated in the S288c reference as centromeric regions, telomeric regions, or LTRs excluded. First, genomic VCF files were obtained with the Haplotype Caller, and joint genotyping of the gVCF files was applied. Using this initial VCF, we applied base quality score recalibration using GATK and called the BAM files again in gVCF mode. After joint calling, in the resulting VCF files, only SNPs or only INDELS were selected. SNPs were filtered according to the parameters [Fay, Liu, Ong, Dunham, Cromie, Jeffery, Ludlow, and Dudley 2019]: QD 3.0; FS > 60.0; MQ 60.0; ReadPosRankSum