fn make_ref_to_aln_map(ref_seq: &[Nuc]) -> Vec<usize> {
  let mut coord_map = Vec::<usize>::with_capacity(ref_seq.len());

  for (i, nuc) in ref_seq.iter().enumerate() {
    if !nuc.is_gap() {
      coord_map.push(i);
    }
  }

  coord_map.shrink_to_fit();
  coord_map
}
/// Makes the "reference to alignment" coordinate map: from reference coordinates to alignment coordinates.
/// Given a position of a letter in the reference sequence, the "reference to alignment" coordinate map allows to
/// lookup the position of the corresponding letter in the aligned sequence.
///
Makes the "reference to alignment" coordinate map: from reference coordinates to alignment coordinates.
Given a position of a letter in the reference sequence, the "reference to alignment" coordinate map allows to
lookup the position of the corresponding letter in the aligned sequence.
Makes the " reference to alignment " coordinate map : from reference coordinates to alignment coordinates . Given a position of a letter in the reference sequence the " reference to alignment " coordinate map allows to lookup the position of the corresponding letter in the aligned sequence .
fn make_ref_to_aln_map(ref_seq: &[Nuc]) -> Vec<usize> {
  let mut coord_map = Vec::<usize>::with_capacity(ref_seq.len());
  for (i, nuc) in ref_seq.iter().enumerate() {
    if !nuc.is_gap() {
      coord_map.push(i);
    }
  }
  coord_map.shrink_to_fit();
  coord_map
}
fn make_ref_to_aln_map ( ref_seq : & [ Nuc ] ) -> Vec < usize > { let mut coord_map = Vec :: < usize > :: with_capacity ( ref_seq . len ( ) ) ; for ( i , nuc ) in ref_seq . iter ( ) . enumerate ( ) { if ! nuc . is_gap ( ) { coord_map . push ( i ) ; } } coord_map . shrink_to_fit ( ) ; coord_map }
Makes the "reference to alignment" coordinate map: from reference coordinates to alignment coordinates.
Makes the " reference to alignment " coordinate map : from reference coordinates to alignment coordinates .
ref_seq
&[Nuc]
ref_seq
&[Nuc]