The SVs were called by svtools as described in https://github.com/hall-lab/svtools/blob/master/Tutorial.md. The R custom codes providing here were used to
- convert sv vcf format into 0 (absent) and 1 (present) SVs via "modify_lumpy.R"
- calculate 500-kb windowed SNP and SV diversity via "sv_snp_diversity.R"
- plot a genome-wide circos via "8k_circos_5group.R"
- plot global site frequency spectrum via "site_frequency_spectrum"
- plot linkage disequilibrium decay via "linkage.R"
- perform k-mean clustering via "k-clustering.R"
- plot the k-mean clusters with geological location via "map_distribution.R"
- plot ADMIXTURE result with k-mean clustering via "plot_admixture_347.R"
- determine cluster-specific SVs via "unique_cluster_snpeff.R"
- perform GO enrichment analysis via "GO_cluster_specific_based_on_snpeff"