pub fn get_peaks(x: &[Homology]) -> Vec<usize>