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
|
#include <iostream>
#include <seqan/basic.h>
#ifndef STDLIB_VS
#include <seqan/blast.h>
using namespace seqan2;
int main()
{
std::string inPath = getAbsolutePath("/tests/blast/plus_comments_defaults.m9");
std::ifstream fin(toCString(inPath), std::ios_base::in | std::ios_base::binary);
auto fit = directionIterator(fin, Input());
typedef std::pair<std::string, std::string> THsp;
std::vector<THsp> hsps;
while (!atEnd(fit))
{
// skip any comment lines
if (!onMatch(fit, BlastTabularLL()))
{
skipUntilMatch(fit, BlastTabularLL());
if (atEnd(fit))
break;
}
// resize output list
resize(hsps, length(hsps)+1);
// read only the first two fields into our variables
readMatch(fit, BlastTabularLL(), back(hsps).first, back(hsps).second);
}
std::sort(std::begin(hsps), std::end(hsps));
auto last = std::unique(std::begin(hsps), std::end(hsps));
hsps.erase(last, hsps.end());
for (THsp const & hsp : hsps)
std::cout << '(' << hsp.first << ", " << hsp.second << ")\n";
return 0;
}
#else
int main()
{
std::cerr << "Demo not run, because of a bug in Microsoft Visual Studio 2015.\n";
return 0;
}
#endif
|
