Remove read kmer databases when we're done with them.

fabiogentilini · May 31, 2019 · ec37af0 · ec37af0
1 parent 9eec1f6
commit ec37af0
Showing 1 changed file with 13 additions and 0 deletions.
diff --git a/src/pipelines/canu/HaplotypeReads.pm b/src/pipelines/canu/HaplotypeReads.pm
@@ -931,6 +931,19 @@ sub haplotypeSubtractCheck ($@) {
     generateReport($asm);
     resetIteration("haplotype-merylSubtractCheck");
 
+    #  Remove databases.
+
+    print STDERR "--\n";
+
+    foreach my $hap (@haplotypes) {
+        if (getGlobal("saveMerCounts") == 0) {
+            print STDERR "-- Removing meryl database '$path/reads-$hap.meryl'.\n";
+            remove_tree("$path/reads-$hap.meryl")   
+        } else {
+            print STDERR "-- Meryl database '$path/reads-$hap.meryl' saved because 'saveMerCounts=true'.\n";
+        }
+    }
+
   allDone:
     stopAfter("meryl-subtract");
 }