I mainly focus on dealing with Reduced Representation Sequence (RRS) data obtained by RAD-seq, MIG-seq and GRAS-Di is assumed (but any other data could be applied).
calc_pairwise_sim.R is the script which reads phylip format generated by Stacks to calculate number of substitutions. You can visualize heat maps as well.
# code block
Rscript xx.R xx.vcf
Rscript xx.R xx.vcf