Fix bug with FASTQ indexing

morispi · Oct 20, 2020 · d9df3d9 · d9df3d9
1 parent 05adb0e
commit d9df3d9
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/src/utils.cpp b/src/utils.cpp
@@ -172,7 +172,7 @@ void indexReads(robin_hood::unordered_map<std::string, std::vector<bool>>& index
 	while (header.length() > 0) {
 		// Get header
 		header.erase(0, 1);
-		// header = splitString(header, " ")[0];
+		header = splitString(header, " ")[0];
 
 		// Get sequence, watching out for multiline FASTA/FASTQ
 		getline(f, seq);
@@ -188,7 +188,6 @@ void indexReads(robin_hood::unordered_map<std::string, std::vector<bool>>& index
 		// Index header/sequence pair
 		std::transform(sequence.begin(), sequence.end(), sequence.begin(), ::toupper);
 		index[header] = fullstr2num(sequence);
-		// cerr << "Indexed " << header << endl;
 
 		// Skip remaining lines if FASTQ
 		if (seq[0] == '+') {