Hi, and welcome!
Take a look at the pegas package which has an amova function (yes, with an m), which returns
a list with a table of sums of square deviations (SSD), mean square deviations (MSD), and the number of degrees of freedom, and a vector of variance components.
I'm not in your field, but recommend looking at this because the package describes itself
Description Functions for reading, writing, plotting, analysing, and manipulating allelic and haplotypic data, including from VCF files, and for the analysis of population nucleotide sequences and micro-satellites including coalescent analyses, linkage disequilibrium, population structure (Fst, Amova) and equilibrium (HWE), haplotype networks, minimum span- ning tree and network, and median-joining networks.
There's also an MMD, mismatch distribution function that
draws a histogram of the frequencies of pairwise distances from a set of DNA sequences.