1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
|
#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);
// Iterate over all leafs, count and print the result
Iterator<FragmentStore<>, AnnotationTree<> >::Type it;
it = begin(store, AnnotationTree<>());
unsigned count = 0;
std::cout << "Number of children for each mRNA: " << std::endl;
// Go down to the first leaf (first child of the first mRNA)
while (goDown(it))
{}
while (!atEnd(it))
{
++count;
// Iterate over all siblings and count
while (goRight(it))
++count;
std::cout << count << std::endl;
count = 0;
// Jump to the next mRNA or gene, go down to its first leaf and count it
if (!atEnd(it))
{
goNext(it);
if (!atEnd(it))
while (goDown(it))
{}
}
}
return 0;
}
|
