7 Convert reference sequences to RSEM format
11 std::string convert(const std::string& rawseq) {
12 int size = (int)rawseq.size();
13 std::string seq = rawseq;
14 for (int i = 0; i < size; i++) {
15 seq[i] = toupper(rawseq[i]);
16 if (seq[i] != 'A' && seq[i] != 'C' && seq[i] != 'G' && seq[i] != 'T') seq[i] = 'N';