From f1eadb317cc2731fc4a2a10c222c684d830919cc Mon Sep 17 00:00:00 2001 From: "Brian P. Walenz" Date: Thu, 25 Feb 2021 19:57:44 -0500 Subject: [PATCH] Request uncompressed lengths for read length histograms. Issue #1906. --- src/pipelines/canu/SequenceStore.pm | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/pipelines/canu/SequenceStore.pm b/src/pipelines/canu/SequenceStore.pm index 6c18ce506..0f4e85698 100644 --- 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";