7173ec16f870dfd0a5ee590354fc604ecfcf1bc9 hiram Thu May 16 11:30:42 2019 -0700 correctly dumping guides no redmine diff --git src/utils/crisprKmers/crisprKmers.c src/utils/crisprKmers/crisprKmers.c index 959d532..f1ff7ba 100644 --- src/utils/crisprKmers/crisprKmers.c +++ src/utils/crisprKmers/crisprKmers.c @@ -1313,31 +1313,34 @@ long long vmPeak = currentVmPeak(); verbose(1, "# vmPeak after scanSequence: %lld kB\n", vmPeak); /* processing those guides */ if (verboseLevel() > 1) { if (offTargets) offFile = mustOpen(offTargets, "w"); /* if ranges have been specified, construct queryList of kmers to measure */ if (rangesHash) { /* result here is two exclusive sets: query, and allGuides * the query guides have been extracted from the allGuides */ queryGuides = rangeExtraction(& allGuides); if (dumpGuides) + { writeGuides(queryGuides, dumpGuides); + return; + } } if (queryGuides) copyToArray(queryGuides); if (allGuides) copyToArray(allGuides); vmPeak = currentVmPeak(); verbose(1, "# vmPeak after copyToArray: %lld kB\n", vmPeak); /* larger example: 62646196 kB */ if (threads > 1) { int gB = vmPeak >> 20; /* convert kB to Gb */ threadCount = threads; verbose(1, "# at %d Gb (%lld kB), running %d threads\n", gB, vmPeak, threadCount); } if (queryGuides) // when range selected some query sequences