These scripts run in R or python3

1.virome_data_run.py: The script was used to remove false positive viral 
genomes identified from vibrant.

2.skip.py: The script was used to filter viral genomes greater or less the
settled length threshold.

3.microdiversity/classified_population.R: The script was used to classify the viral populations as multi-zonal, regional or local viral populations according to their distribution in soils with various land use type.

4.microdiversity/microdiveristy.R: The script was used to calculate the 
average microdiversity of multi-zonal, regional or local viral populations
in each samples.

5.microdiveristy/microdiveristy_multi.R: The script was used to calculate
the average microdiversity of multi-zonal viral populations in each samples.

6.microdiveristy/pNpS.R: The script was used to calculate the average pNpS value
of multi-zonal, regional or local viral populations in each samples.

7.microdiveristy/all_positive_gene.R:  The script was used to calculate the positive selected genes