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
|
//![main]
#include <iostream>
#include <seqan/align.h>
using namespace seqan;
int main()
{
//![main]
//![init]
Align<String<AminoAcid> > ali;
resize(rows(ali), 2);
assignSource(row(ali, 0), "PNCFDAKQRTASRPL");
assignSource(row(ali, 1), "CFDKQKNNRTATRDTA");
//![init]
//![ali]
Score<int> sc(3, -2, -1, -5);
unsigned count = 0;
LocalAlignmentEnumerator<Score<int>, Unbanded> enumerator(sc);
while (nextLocalAlignment(ali, enumerator) && count < 3)
{
std::cout << "Score = " << getScore(enumerator) << std::endl;
std::cout << ali;
++count;
}
return 0;
}
//![ali]
|
