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
|
#include <cstddef>
#include <assert.h>
#include "../../../pbdata/FASTQSequence.hpp"
#include "SAMQVConversion.hpp"
void QualityVectorToPrintable(unsigned char *data, int length) {
if (data == NULL) {
return;
}
for (int i = 0; i < length; i++) {
data[i] = (((unsigned char)data[i]) == MAX_STORED_QUALITY) ? MAX_PRINTED_QUALITY : (unsigned char)data[i];
data[i] = (((unsigned char)data[i]) == SENTINAL) ? MAP_SENTINAL : (unsigned char)data[i];
assert(data[i] != 255);
}
}
void QualityStringToStored(unsigned char *data, int length) {
if (data == NULL) {
return;
}
for (int i = 0; i < length; i++) {
data[i] = data[i] - FASTQSequence::charToQuality;
data[i] = ((unsigned char)data[i]) == MAX_PRINTED_QUALITY ? MAX_STORED_QUALITY : (unsigned char)data[i];
data[i] = ((unsigned char)data[i]) == MAP_SENTINAL ? SENTINAL : (unsigned char)data[i];
}
}
|
