Function getTreeSequences

RResolver/SequenceTree.cpp:172–205

}

std::vector<Sequence>
getTreeSequences(
    const ContigNode& start,
    const int overlap,
    const int maxLength,
    const bool forward,
    const int maxPaths)
{
	assert(overlap >= 0);
	assert(maxLength > 0);
	assert(maxPaths >= 1);

	auto traces = getTreeTraces(start, overlap, maxLength, forward, maxPaths);

	std::vector<Sequence> sequences;
	sequences.reserve(traces.size());

	for (const auto& trace : traces) {
		Sequence sequence;
		sequence.reserve(trace.size() * EXPECTED_BASES_PER_NODE);
		if (forward) {
			for (Trace::const_iterator it = trace.begin(); !(it == trace.end()); it++) {
				sequence += it->treeSequence();
			}
		} else {
			for (Trace::const_reverse_iterator it = trace.rbegin(); !(it == trace.rend()); it++) {
				sequence += it->treeSequence();
			}
		}
		sequences.push_back(sequence);
	}

	return sequences;
}