Print depth median for long reads.

seqan · Nov 29, 2019 · b67a6d9 · b67a6d9
1 parent 13d81af
commit b67a6d9
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/readdepth.cpp b/src/readdepth.cpp
@@ -89,6 +89,8 @@ void ReadDepth::prepAfterHeaderParsing(seqan::BamHeader& header, seqan::BamFileI
         std::sort(backgroundDepth.begin(), backgroundDepth.end(), [](auto &left, auto &right) {return left < right;} );
         double medianDepth = MID_ELEMENT(backgroundDepth);
         this->setMedianDepth( medianDepth );
+
+        printMessage("Depth median for long reads: " + std::to_string(medianDepth))
     }
 }