Convert a path in the de Bruijn graph to a sequence */
| 20 | |
| 21 | /** Convert a path in the de Bruijn graph to a sequence */ |
| 22 | TEST(BloomDBG, pathToSeq) |
| 23 | { |
| 24 | const string inputSeq = "ACGTAC"; |
| 25 | const string spacedSeed = "10001"; |
| 26 | const unsigned k = 5; |
| 27 | const unsigned numHashes = 2; |
| 28 | |
| 29 | MaskedKmer::setLength(k); |
| 30 | MaskedKmer::setMask(spacedSeed); |
| 31 | |
| 32 | Path<BloomDBG::Vertex> path = |
| 33 | BloomDBG::seqToPath(inputSeq, k, numHashes); |
| 34 | ASSERT_EQ(2U, path.size()); |
| 35 | |
| 36 | string outputSeq = BloomDBG::pathToSeq(path, k); |
| 37 | ASSERT_EQ("ACNNAC", outputSeq); |
| 38 | } |