pub const fn calculate_aa_alignment_params(qry_gaps: &GapCounts, ref_gaps: &GapCounts) -> PeptideAlignmentParams {
  const BASE_BAND_WIDTH: usize = 5; // An arbitrary magic number to give some additional room for alignment
  let band_width = (qry_gaps.internal + ref_gaps.internal) / 3 + BASE_BAND_WIDTH;
  let mean_shift =
    (qry_gaps.leading as i32 - ref_gaps.leading as i32) / 3 + (qry_gaps.internal as i32 - ref_gaps.internal as i32) / 6;
  PeptideAlignmentParams { band_width, mean_shift }
}
/// Deduces aminoacid alignment params from nucleotide alignment.
/// Gap counts from nucleotide alignment are used to estimate rough parameters. This is by contrast to nucleotide
/// alignment, where these parameters are estimated using seed matching.
Deduces aminoacid alignment params from nucleotide alignment.
Gap counts from nucleotide alignment are used to estimate rough parameters. This is by contrast to nucleotide
alignment, where these parameters are estimated using seed matching.
Deduces aminoacid alignment params from nucleotide alignment . Gap counts from nucleotide alignment are used to estimate rough parameters . This is by contrast to nucleotide alignment where these parameters are estimated using seed matching .
pub const fn calculate_aa_alignment_params(qry_gaps: &GapCounts, ref_gaps: &GapCounts) -> PeptideAlignmentParams {
  const BASE_BAND_WIDTH: usize = 5; // An arbitrary magic number to give some additional room for alignment
  let band_width = (qry_gaps.internal + ref_gaps.internal) / 3 + BASE_BAND_WIDTH;
  let mean_shift =
    (qry_gaps.leading as i32 - ref_gaps.leading as i32) / 3 + (qry_gaps.internal as i32 - ref_gaps.internal as i32) / 6;
  PeptideAlignmentParams { band_width, mean_shift }
}
pub const fn calculate_aa_alignment_params ( qry_gaps : & GapCounts , ref_gaps : & GapCounts ) -> PeptideAlignmentParams { const BASE_BAND_WIDTH : usize = 5 ; // An arbitrary magic number to give some additional room for alignment let band_width = ( qry_gaps . internal + ref_gaps . internal ) / 3 + BASE_BAND_WIDTH ; let mean_shift = ( qry_gaps . leading as i32 - ref_gaps . leading as i32 ) / 3 + ( qry_gaps . internal as i32 - ref_gaps . internal as i32 ) / 6 ; PeptideAlignmentParams { band_width , mean_shift } }
Deduces aminoacid alignment params from nucleotide alignment.
Deduces aminoacid alignment params from nucleotide alignment .
qry_gaps
&GapCounts
ref_gaps
&GapCounts
qry_gaps
&GapCounts
ref_gaps
&GapCounts