pub fn consalign<T, U>(
    fasta_records: &FastaRecords,
    align_prob_mats_with_rna_id_pairs: &SparseProbMatsWithRnaIdPairs<T>,
    bpp_mats: &SparseProbMats<T>,
    feature_scores: &FeatureCountsPosterior,
    insert_prob_set_pairs_with_rna_id_pairs: &ProbSetPairsWithRnaIdPairs
) -> MeaStructAlign<T, U>where
    T: Unsigned + PrimInt + Hash + FromPrimitive + Integer + Ord + Sync + Send + Display,
    U: Unsigned + PrimInt + Hash + FromPrimitive + Integer + Ord + Sync + Send + Display,