Return whether the read seq is in colour space. */
| 120 | |
| 121 | /** Return whether the read seq is in colour space. */ |
| 122 | static bool isColourSpace(const string& seq) |
| 123 | { |
| 124 | assert(!seq.empty()); |
| 125 | size_t i = seq.find_first_of("ACGTacgt0123", 1); |
| 126 | return i != string::npos && isdigit(seq[i]); |
| 127 | } |
| 128 | |
| 129 | /** Read a single record. */ |
| 130 | Sequence FastaReader::read(string& id, string& comment, |