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
|
//![main]
#include <iostream>
#include <seqan/align.h>
using namespace seqan2;
int main()
{
typedef String<Dna> TSequence; // sequence type
typedef StringSet<TSequence> TStringSet; // container for strings
typedef StringSet<TSequence, Dependent<> > TDepStringSet; // dependent string set
typedef Graph<Alignment<TDepStringSet> > TAlignGraph; // alignment graph
//![main]
//![init]
TSequence seq1 = "AAATGACGGATTG";
TSequence seq2 = "TGGGA";
TStringSet sequences;
appendValue(sequences, seq1);
appendValue(sequences, seq2);
TAlignGraph alignG(sequences);
//![init]
//![alignment]
int score = globalAlignment(alignG, Score<int, Simple>(1, -1, -1), AlignConfig<true, false, false, true>());
std::cout << "Score: " << score << std::endl;
std::cout << alignG << std::endl;
return 0;
}
//![alignment]
|
