INNER CODE UNIT · C++
dups
refresh-bio/FAMSA · src/msa.cpp:545
dups = sequences.size() - mapped_seqs.size();
LOG_VERBOSE << mapped_seqs.size() << "/" << sequences.size() << " sequences retained." << endl;
}
// only one unique sequence - move input to output end exit
if (mapped_seqs.size() == 1) {
final_profile = new CProfile(¶ms);
for (int i = 0; i < (int)sequences.size(); ++i) {
final_profile->AppendRawSequence(CGappedSequence(move(sequences[i])));
}
return true;
}
// store mappings and temporarily reset numerical identifiers (to make medoid trees work)
for (int i = 0; i < (int)mapped_seqs.size(); ++i) {
mapped_seqs[i]->sequence_no = i;
}