1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
//![includes]
#include <iostream>
#include <seqan/index.h>
using namespace seqan2;
//![includes]
//![initialization]
int main()
{
Index<DnaString, IndexQGram<OneGappedShape> > index("CATGATTACATA");
stringToShape(indexShape(index), "1101");
//![initialization]
//![output]
hash(indexShape(index), "ATCA");
for (unsigned i = 0; i < length(getOccurrences(index, indexShape(index))); ++i)
std::cout << getOccurrences(index, indexShape(index))[i] << std::endl;
return 0;
}
//![output]
|
