bool IsOriginalClusterLost(size_t orig_cluster_ind) {
auto reads = original_clusters_.GetReadsByCluster(orig_cluster_ind);
for(auto read = reads.begin(); read != reads.end(); read++)
if(constructed_clusters_.ReadExists(original_reads_[*read]))
return false;
return true;
}
vector<string> GetReadNamesFromCluster(size_t cluster_ind, Clusters &clusters, vector<string> &reads) {
vector<size_t> read_inds = clusters.GetReadsByCluster(cluster_ind);
vector<string> read_names;
for(auto it = read_inds.begin(); it != read_inds.end(); it++)
read_names.push_back(reads[*it]);
return read_names;
}
