1d76bf881d3e18580739e44698a0d1a499c2f9f1 hiram Fri Oct 2 10:31:39 2015 -0700 fixup gcc warnings for -Wunused-but-set-variable refs #16121 diff --git src/lib/dnautil.c src/lib/dnautil.c index 6b08ccd..e12ad44 100644 --- src/lib/dnautil.c +++ src/lib/dnautil.c @@ -984,60 +984,58 @@ * trimmed too. Returns number of bases trimmed. Leaves very last a. */ { return findTailPolyAMaybeMask(dna, size, TRUE, FALSE); } static int findHeadPolyTMaybeMask(DNA *dna, int size, boolean doMask, boolean loose) /* Return size of PolyT at start (if present); mask to 'n' if specified. * This allows a few non-T's as noise to be trimmed too, but skips last * two tt for revcomp'd taa stop codon. */ { int i; int score = 10; int bestScore = 10; int bestPos = -1; -int pastPoly = 0; int trimSize = 0; for (i=0; i<size; ++i) { DNA b = dna[i]; if (b == 'n' || b == 'N') continue; if (score > 20) score = 20; if (b == 't' || b == 'T') { score += 1; if (score >= bestScore) { bestScore = score; bestPos = i; } else if (loose && score >= (bestScore - 8)) { /* If loose, keep extending even if score isn't back up to best. */ bestPos = i; } } else { score -= 10; } if (score < 0) { - pastPoly = i; break; } } if (bestPos >= 0) { trimSize = bestPos+1 - 2; // Leave two for aa in taa stop codon if (trimSize > 0) { if (doMask) memset(dna, 'n', trimSize); } else trimSize = 0; } return trimSize;