From 8e13416bd89df5a0d841f1c64b9c9319deaf90d5 Mon Sep 17 00:00:00 2001 From: Adam Novak Date: Fri, 31 May 2024 14:33:24 -0700 Subject: [PATCH] Note when parts of a dumped chain go off of all the reference paths --- src/minimizer_mapper_from_chains.cpp | 15 +++++++++++++++ 1 file changed, 15 insertions(+) diff --git a/src/minimizer_mapper_from_chains.cpp b/src/minimizer_mapper_from_chains.cpp index a84ded02058..b8d2dc77f51 100644 --- a/src/minimizer_mapper_from_chains.cpp +++ b/src/minimizer_mapper_from_chains.cpp @@ -238,6 +238,21 @@ void MinimizerMapper::dump_debug_dotplot(const std::string& name, const VectorVi exp.field(minimizers[seed.source].pin_offset()); } } + if (offsets.empty()) { + // Note that we don't actually have a position + exp.line(); + if (!marker.empty()) { + // Sentinel and a marker and a subscript + exp.field("NO_PATH-" + marker + "-" + std::to_string(run_number)); + } else { + // Sentinel alone + exp.field("NO_PATH"); + } + // Put it at 0 on no path + exp.field(0); + // Offset in read *of the pin point* (not of the forward-strand start of the minimizer) + exp.field(minimizers[seed.source].pin_offset()); + } } }