MCPcopy Create free account
hub / github.com/BirolLab/abyss / getSequence

Function getSequence

MergePaths/PathConsensus.cpp:191–202  ·  view source on GitHub ↗

Return the sequence of the specified contig node. The sequence * may be ambiguous or reverse complemented. */

Source from the content-addressed store, hash-verified

189 * may be ambiguous or reverse complemented.
190 */
191static const Sequence getSequence(ContigNode id)
192{
193 if (id.ambiguous()) {
194 string s(id.ambiguousSequence());
195 if (s.length() < opt::k)
196 transform(s.begin(), s.end(), s.begin(), ::tolower);
197 return string(opt::k - 1, 'N') + s;
198 } else {
199 string seq(g_contigs[id.id()]);
200 return id.sense() ? reverseComplement(seq) : seq;
201 }
202}
203
204/** Return the distance from vertex u to v. */
205static int getDistance(const Graph& g,

Callers 2

mergePathFunction · 0.70
alignMultiFunction · 0.70

Calls 8

ambiguousMethod · 0.80
ambiguousSequenceMethod · 0.80
idMethod · 0.80
senseMethod · 0.80
reverseComplementFunction · 0.50
lengthMethod · 0.45
beginMethod · 0.45
endMethod · 0.45

Tested by

no test coverage detected