Request uncompressed lengths for read length histograms. Issue #1906.

marbl · Feb 26, 2021 · f1eadb3 · f1eadb3
1 parent 2f73439
commit f1eadb3
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/pipelines/canu/SequenceStore.pm b/src/pipelines/canu/SequenceStore.pm
@@ -418,7 +418,7 @@ sub generateReadLengthHistogram ($$) {
         $cmd .= "  -raw \\\n"                  if ($tag eq "cor");
         $cmd .= "  -corrected \\\n"            if ($tag eq "obt");
         $cmd .= "  -corrected -trimmed \\\n"   if ($tag eq "utg");
-        $cmd .= "  -histogram \\\n";
+        $cmd .= "  -uncompressed \\\n";
         $cmd .= "  -lengths \\\n";
         $cmd .= "> ./$asm.seqStore/readlengths-$tag.dat \\\n";
         $cmd .= "2> ./$asm.seqStore/readlengths-$tag.err \n";
@@ -462,6 +462,7 @@ sub generateReadLengthHistogram ($$) {
         $cmd .= "  -raw \\\n"                  if ($tag eq "cor");
         $cmd .= "  -corrected \\\n"            if ($tag eq "obt");
         $cmd .= "  -corrected -trimmed \\\n"   if ($tag eq "utg");
+        $cmd .= "  -uncompressed \\\n";
         $cmd .= "  -histogram \\\n";
         $cmd .= "> ./$asm.seqStore/readlengths-$tag.txt \\\n";
         $cmd .= "2> ./$asm.seqStore/readlengths-$tag.err \n";