Permalink
Browse files

Remove quality scores after correction

  • Loading branch information...
1 parent 2dc5e59 commit 604f8856920487b8beee9d89ab8a96d6908f0bc6 @jts committed Jan 14, 2013
Showing with 1 addition and 1 deletion.
  1. +1 −1 src/Algorithm/ErrorCorrectProcess.cpp
@@ -375,7 +375,6 @@ ErrorCorrectResult ErrorCorrectProcess::kmerCorrection(const SequenceWorkItem& w
const double MIN_TRIM_RATIO = 0.65;
result.correctSequence = readSequence.substr(0, lastSolidBase);
double len_ratio = (double)result.correctSequence.length() / workItem.read.seq.length();
-// printf("LSB: %d TRIMMED TO [%zu %zu] = %lf %s\n", lastSolidBase, result.correctSequence.length(), workItem.read.seq.length(), len_ratio, result.correctSequence.toString().c_str());
if(len_ratio >= MIN_TRIM_RATIO)
result.kmerQC = true;
else
@@ -770,6 +769,7 @@ void ErrorCorrectPostProcess::process(const SequenceWorkItem& item, const ErrorC
SeqRecord record = item.read;
record.seq = result.correctSequence;
+ record.qual = "";
if(readQCPass || m_pDiscardWriter == NULL)
{

0 comments on commit 604f885

Please sign in to comment.