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
|
#include <iostream>
#include <seqan/store.h>
#include <seqan/misc/interval_tree.h>
#include <seqan/parallel.h>
using namespace seqan2;
// define used types
typedef FragmentStore<> TStore;
//
// 2. Load annotations and alignments from files
//
bool loadFiles(TStore & store, std::string const & annotationFileName, std::string const & alignmentFileName)
{
(void) store;
(void) annotationFileName;
(void) alignmentFileName;
// INSERT YOUR CODE HERE ...
//
return true;
}
//![main]
int main(int, char const **)
{
TStore store;
std::string annotationFileName = getAbsolutePath("demos/tutorial/simple_rna_seq/example.gtf");
std::string alignmentFileName = getAbsolutePath("demos/tutorial/simple_rna_seq/example.sam");
if (!loadFiles(store, annotationFileName, annotationFileName))
return 1;
return 0;
}
//![main]
|
