1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
#include <seqan/stream.h>
#include <seqan/sequence.h>
#include <seqan/file.h>
using namespace seqan2;
int main()
{
String<Dna5> nucleotides = "AGTCGTGNNANCT";
String<Dna5> lesser;
String<Dna5> greater;
for (unsigned i = 0; i < length(nucleotides); ++i){
if (nucleotides[i] < 'G')
appendValue(lesser, nucleotides[i]);
else if (nucleotides[i] > 'G')
appendValue(greater, nucleotides[i]);
}
std::cout << "Lesser nucleotides: " << lesser << std::endl;
std::cout << "Greater nucleotides: " << greater << std::endl;
}
|
