1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
#include <iostream>
#include <seqan/sequence.h>
#include <seqan/stream.h>
#include <seqan/store.h>
using namespace seqan;
int main()
{
CharString fileName = getAbsolutePath("demos/tutorial/genome_annotations/assignment_annotations.gtf");
GffFileIn file(toCString(fileName));
FragmentStore<> store;
readRecords(store, file);
// Create AnnotationTree iterator
Iterator<FragmentStore<>, AnnotationTree<> >::Type it;
it = begin(store, AnnotationTree<>());
// Move iterator one node down
goDown(it);
std::cout << "Is leaf: " << isLeaf(it) << std::endl;
return 0;
}
|
