X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=trimseqscommand.cpp;h=2dc07794be5abb2f239e4b228c3aacbf51c51b6e;hb=d04f948b1a2a1a2984fc4a45d04403b8c121c5bc;hp=d5df4730c9f118966de5f3514b4ec276faffce05;hpb=5b3bedc6fbfed63854cc6a7d47454580b2734c08;p=mothur.git diff --git a/trimseqscommand.cpp b/trimseqscommand.cpp index d5df473..2dc0779 100644 --- a/trimseqscommand.cpp +++ b/trimseqscommand.cpp @@ -556,16 +556,16 @@ int TrimSeqsCommand::driverCreateTrim(string filename, string qFileName, string if(qFileName != ""){ + int origLength = currSeq.getNumBases(); if(qThreshold != 0) { success = currQual.stripQualThreshold(currSeq, qThreshold); } else if(qAverage != 0) { success = currQual.cullQualAverage(currSeq, qAverage); } else if(qRollAverage != 0) { success = currQual.stripQualRollingAverage(currSeq, qRollAverage); } else if(qWindowAverage != 0){ success = currQual.stripQualWindowAverage(currSeq, qWindowStep, qWindowSize, qWindowAverage); } else { success = 1; } - -// if (qtrim == 1 && (origSeq.length() != currSeq.getUnaligned().length())) { -// success = 0; //if you don't want to trim and the sequence does not meet quality requirements, move to scrap -// } + + //you don't want to trim, if it fails above then scrap it + if ((!qtrim) && (origLength != currSeq.getNumBases())) { success = 0; } if(!success) { trashCode += 'q'; } }