pub fn calculate_genome_stats(fasta_path: &str) -> GenomeAssemblyStats