revcomp for uracil. fixes issue-682
@@ -39,6 +39,10 @@ void reverseComplement(string &seq) {
case 'T':
seq[i] = 'A';
break;
+ // issue 682
+ case 'U':
+ seq[i] = 'A';
+ break;
case 'a':
seq[i] = 't';
break;
@@ -51,6 +55,10 @@ void reverseComplement(string &seq) {
case 't':
seq[i] = 'a';
break;
+ // issue 682
+ case 'u':
+ seq[i] = 'a';
+ break;
// unknown (N)
case 'N':
@@ -185,6 +193,8 @@ void getDnaContent(const string &seq, int &a, int &c, int &g, int &t, int &n, in
break;
case 'T':
case 't':
+ case 'U':
+ case 'u':
t++;
break;
case 'N':
