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
|
#include <alignment/datastructures/alignmentset/SAMQVConversion.hpp>
#include <pbdata/FASTQSequence.hpp>
#include <cassert>
#include <cstddef>
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];
}
}
|