1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
//![start]
#include <seqan3/alphabet/nucleotide/all.hpp>
#include <seqan3/core/debug_stream.hpp>
#include <seqan3/std/ranges>
using seqan3::operator""_dna5;
//![start]
auto my_reverse_complement = std::views::reverse | std::views::transform([] (auto const d)
{
return seqan3::complement(d);
});
//![end]
int main()
{
std::vector<seqan3::dna5> vec{"ACCAGATTA"_dna5};
seqan3::debug_stream << vec << '\n'; // will print "ACCAGATTA"
auto v = vec | my_reverse_complement;
seqan3::debug_stream << v << '\n'; // prints "TAATCTGGT"
}
//![end]
|
