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
|
#include <iostream>
#include <seqan/modifier.h>
#include <seqan/stream.h>
using namespace seqan2;
int main()
{
// A nested modifier.
typedef ModifiedString<String<Dna>, ModComplementDna> TMyComplement;
typedef ModifiedString<TMyComplement, ModReverse> TMyReverseComplement;
// The original string.
String<Dna> myString = "attacgg";
// A reverse complemented string.
TMyReverseComplement myReverseComplement(myString);
std::cout << myString << "\n"
<< myReverseComplement << "\n";
replace(myString, 1, 1, "cgt");
std::cout << myString << "\n"
<< myReverseComplement << "\n"
<< DnaStringReverseComplement(myString) << "\n";
return 0;
}
|
