Return the sequence of the specified contig node. The sequence * may be ambiguous or reverse complemented. */
| 189 | * may be ambiguous or reverse complemented. |
| 190 | */ |
| 191 | static 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. */ |
| 205 | static int getDistance(const Graph& g, |
no test coverage detected