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
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
|
#include "cleanup.hpp"
seqan3::cleanup index_file{"index.file"};
#if SEQAN3_WITH_CEREAL
#include <fstream>
#include <cereal/archives/binary.hpp>
#endif //SEQAN3_WITH_CEREAL
#include <seqan3/core/debug_stream.hpp>
#include <seqan3/search/search.hpp>
#include <seqan3/search/fm_index/bi_fm_index.hpp>
#include <seqan3/search/fm_index/fm_index.hpp>
#include <seqan3/std/filesystem>
using namespace std::string_literals;
int main()
{
{
//![text_collection]
std::vector<std::string> text{{"Garfield the fat cat without a hat."},
{"He is infinite, he is eternal."},
{"Yet another text I have to think about."}};
seqan3::fm_index index{text};
seqan3::bi_fm_index bi_index{text};
//![text_collection]
}
#if SEQAN3_WITH_CEREAL
{
//![store]
#include <fstream> // for writing/reading files
#include <cereal/archives/binary.hpp> // for storing/loading indices via cereal
std::string text{"Garfield the fat cat without a hat."};
seqan3::fm_index index{text};
{
std::ofstream os{"index.file", std::ios::binary};
cereal::BinaryOutputArchive oarchive{os};
oarchive(index);
}
//![store]
}
{
//![load]
// we need to tell the index that we work on a single text and a `char` alphabet before loading
seqan3::fm_index<char, seqan3::text_layout::single> index;
{
std::ifstream is{"index.file", std::ios::binary};
cereal::BinaryInputArchive iarchive{is};
iarchive(index);
}
//![load]
}
#endif //SEQAN3_WITH_CEREAL
{
//![error_search]
std::string text{"Garfield the fat cat without a hat."};
seqan3::fm_index index{text};
seqan3::configuration const cfg = seqan3::search_cfg::max_error_total{seqan3::search_cfg::error_count{1}} |
seqan3::search_cfg::max_error_substitution{seqan3::search_cfg::error_count{0}} |
seqan3::search_cfg::max_error_insertion{seqan3::search_cfg::error_count{1}} |
seqan3::search_cfg::max_error_deletion{seqan3::search_cfg::error_count{1}};
seqan3::debug_stream << search("cat"s, index, cfg) << '\n';
// prints: [<query_id:0, reference_id:0, reference_pos:14>,
// <query_id:0, reference_id:0, reference_pos:17>,
// <query_id:0, reference_id:0, reference_pos:18>,
// <query_id:0, reference_id:0, reference_pos:32>]
//![error_search]
}
{
//![multiple_queries]
std::string text{"Garfield the fat cat without a hat."};
seqan3::fm_index index{text};
std::vector<std::string> query{"cat"s, "hat"s};
seqan3::debug_stream << search(query, index) << '\n';
// prints: [<query_id:0, reference_id:0, reference_pos:17>,
// <query_id:1, reference_id:0, reference_pos:31>]
//![multiple_queries]
}
{
//![error_sum]
seqan3::configuration const cfg = seqan3::search_cfg::max_error_total{seqan3::search_cfg::error_count{2}} |
seqan3::search_cfg::max_error_substitution{seqan3::search_cfg::error_count{2}} |
seqan3::search_cfg::max_error_insertion{seqan3::search_cfg::error_count{1}} |
seqan3::search_cfg::max_error_deletion{seqan3::search_cfg::error_count{1}};
//![error_sum]
}
{
//![hit_best]
seqan3::configuration const cfg = seqan3::search_cfg::max_error_total{seqan3::search_cfg::error_count{1}} |
seqan3::search_cfg::max_error_substitution{seqan3::search_cfg::error_count{0}} |
seqan3::search_cfg::max_error_insertion{seqan3::search_cfg::error_count{1}} |
seqan3::search_cfg::max_error_deletion{seqan3::search_cfg::error_count{1}} |
seqan3::search_cfg::hit_single_best{};
//![hit_best]
}
{
//![hit_strata]
seqan3::configuration cfg = seqan3::search_cfg::max_error_total{seqan3::search_cfg::error_count{2}} |
seqan3::search_cfg::max_error_substitution{seqan3::search_cfg::error_count{0}} |
seqan3::search_cfg::max_error_insertion{seqan3::search_cfg::error_count{1}} |
seqan3::search_cfg::max_error_deletion{seqan3::search_cfg::error_count{1}} |
seqan3::search_cfg::hit_strata{2};
using seqan3::get; // Required to be able to find the correct get function.
get<seqan3::search_cfg::hit_strata>(cfg).stratum = 1; // The stratum is now 1 and not 2 anymore.
//![hit_strata]
}
//![hit_dynamic]
seqan3::search_cfg::hit hit_dynamic{seqan3::search_cfg::hit_all{}}; // Initialise with hit_all configuration.
bool hit_with_strata = static_cast<bool>(std::rand() & 1); // Either false or true.
if (hit_with_strata)
hit_dynamic = seqan3::search_cfg::hit_strata{2}; // Search instead with strata mode.
seqan3::configuration const cfg = seqan3::search_cfg::max_error_total{seqan3::search_cfg::error_count{2}} |
seqan3::search_cfg::max_error_substitution{seqan3::search_cfg::error_count{0}} |
seqan3::search_cfg::max_error_insertion{seqan3::search_cfg::error_count{1}} |
seqan3::search_cfg::max_error_deletion{seqan3::search_cfg::error_count{1}} |
hit_dynamic; // Build the configuration by adding the dynamic hit configuration.
//![hit_dynamic]
}
|
