-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy pathtypes.h
More file actions
105 lines (86 loc) · 2.63 KB
/
Copy pathtypes.h
File metadata and controls
105 lines (86 loc) · 2.63 KB
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
// -*- c++ -*-
/**
*
* LightStringGraph
*
* Lightweight String Graph Construction.
*
* Copyright (C) 2013, 2014 Stefano Beretta, Yuri Pirola, Marco Previtali
*
* Distributed under the terms of the GNU General Public License (or the Lesser
* GPL).
*
* This file is part of LightStringGraph.
*
* LighStringGraph is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* LightStringGraph is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with LightStringGraph. If not, see <http://www.gnu.org/licenses/>.
*
**/
#ifndef TYPES_H
#define TYPES_H
#include <vector>
#include <cstdint>
using std::vector;
#define QIntervalManager IntervalManager< QInterval >
#define SameLengthArcIntervalManager IntervalManager< ArcInterval >
#define SameLengthEdgeLabelIntervalManager IntervalManagerRLE< EdgeLabelInterval >
#define BasicArcIntervalManager MultiIntervalManager<ArcInterval>
// As defined in BEETL/src/shared/Tools.hh (ElementType)
// Type to represent: Number of sequences
// below limits to 4 billion reads max - change to uint64_t for more
typedef uint32_t SequenceNumber;
// Type to represent: Sequence length (in biologic case 100)
typedef uint8_t SequenceLength;
struct GSAEntry
{
SequenceLength sa; // suffix array position
// SequenceNumber numSeq; // sequence number
};
// END of definition taken from BEETL
#define PRINT_SL( x ) static_cast<uint32_t>( (x) )
// Position on a BWT
typedef uint64_t BWTPosition;
// Possible value of LCP
typedef SequenceLength LCPValue;
// Edge Length
typedef SequenceLength EdgeLength; // should be enough for short reads
// Nucleotide counter
typedef uint64_t NucleoCounter;
// Possible Nucleotide in sequence
enum Nucleotide
{
BASE_$,
BASE_A,
BASE_C,
BASE_G,
BASE_N,
BASE_T,
BASE_Z,
// NOT_IN_ALPHABET,
ALPHABET_SIZE
};
#define CHAR_BASE_$ '$'
#define CHAR_BASE_A 'A'
#define CHAR_BASE_C 'C'
#define CHAR_BASE_G 'G'
#define CHAR_BASE_N 'N'
#define CHAR_BASE_T 'T'
#define CHAR_BASE_Z 'Z'
#define CHAR_ALPHABET_SIZE '\0'
#define CHAR_BASE_A_LOW 'a'
#define CHAR_BASE_C_LOW 'c'
#define CHAR_BASE_G_LOW 'g'
#define CHAR_BASE_N_LOW 'n'
#define CHAR_BASE_T_LOW 't'
#define CHAR_BASE_Z_LOW 'z'
#endif