///|
/// Count normalized bases in the ORIGIN sequence.
pub fn Record::sequence_stats(self : Record) -> SequenceStats {
  let mut a = 0
  let mut c = 0
  let mut g = 0
  let mut t = 0
  let mut n = 0
  for ch in self.origin {
    match ch {
      'A' => a = a + 1
      'C' => c = c + 1
      'G' => g = g + 1
      'T' | 'U' => t = t + 1
      _ => n = n + 1
    }
  }
  let length = a + c + g + t + n
  let gc_percent = if length == 0 {
    0.0
  } else {
    (g + c).to_double() * 100.0 / length.to_double()
  }
  { length, a, c, g, t, n, gc_percent }
}