# How to create figure 4

# Need to do barcode error correction

# Then need to run calculate_total_model_snp_snp_associations.py (loop over species)
# Then need to run calculate_cross_species_snp_barcode_statistics.py (will create an intermediate file)
# The plot_cross_species_gametes.py
# Can also run plot_sixpanel_traficplot.py example_fourgamete_figure.settings.txt
# Can also run plot_vulgatus_figure.py

# Also need to run calculate_raw_total_snp_snp_associations.py 
# for the shared vs distance plot (loop over species)
# plot_distance_counts.py