2.3 Fixation index (FST) calculation
For each dataset, we used ANGSD to generate unbiased allele frequency
spectra of our populations (AFS) and from that, calculated the weightedF ST for each gene (Korneliussen et al., 2014).
While assigning SNP loci to genes, we extended the annotated gene
boundaries one-thousand bases at the start and end of the gene, to
capture the upstream regulatory region and the downstream untranslated
RNA region that might be missing in the gene annotation.