if other genomes are needed use https://github.com/BushmanLab/intSiteCaller/blob/master/programFlow.R get_reference_genome <- function(reference_genome)