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
|
#include <iostream>
#include <seqan/stream.h>
#include <seqan/modifier.h>
using namespace seqan2;
struct ConvertCT :
public std::function<Dna(Dna)>
{
inline Dna operator()(Dna x) const
{
if (x == 'C') return 'T';
return x;
}
};
int main()
{
typedef String<Dna> TSequence;
TSequence seq1 = "CCCGGCATCATCC";
TSequence seq2 = "CTTGGCATTATTC";
std::cout << seq1 << std::endl;
std::cout << seq2 << std::endl;
std::cout << std::endl;
typedef ModifiedString<TSequence, ModView<ConvertCT> > TModCT;
TModCT modCT1(seq1);
TModCT modCT2(seq2);
std::cout << modCT1 << std::endl;
std::cout << modCT2 << std::endl;
return 0;
}
|
