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
|
#include <iostream>
#include <seqan/sequence.h>
#include <seqan/seq_io.h>
using namespace seqan2;
int main(int argc, char const ** argv)
{
if (argc != 2)
{
std::cerr << "USAGE: build_fai FILE.fa\n";
return 0;
}
FaiIndex faiIndex;
if (!build(faiIndex, argv[1]))
{
std::cerr << "ERROR: Could not build FAI index for file " << argv[1] << ".\n";
return 0;
}
CharString faiFilename = argv[1];
append(faiFilename, ".fai");
if (!save(faiIndex, toCString(faiFilename)))
{
std::cerr << "ERROR: Could not write the index to file!\n";
return 0;
}
std::cout << "Index file " << faiFilename << " was successfully created.\n";
return 0;
}
|
