diff --git a/src/readdepth.cpp b/src/readdepth.cpp index 0176ca8..b637092 100644 --- 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)) } }