INNER CODE UNIT · C++
duplicate
refresh-bio/FAMSA · src/msa.cpp:643
CGappedSequence* duplicate = new CGappedSequence(*representative);
duplicate->id = std::move(current.id);
duplicate->original_no = current.original_no;
ordered_alignment[current.original_no] = duplicate;
}
}
final_profile->data = std::move(ordered_alignment);
}
// store stats
if (params.areStatsStored()) {
// this is time consuming - only in very verbose
if (params.very_verbose_mode) {
int64_t sackin = tree.calculateSackinIndex();
statistics.put("guide_tree.sackin", sackin);
statistics.put("guide_tree.sackin_norm", sackin / (double)sequences.size());