Display hifiasm coverage in a user-friendly way

hifiasm shows coverage as in rd:i:n tag where coverage is n+1.
asl · Jan 22, 2025 · 051621a · 051621a
1 parent bbb9ecb
commit 051621a
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/graph/assemblygraphbuilder.cpp b/graph/assemblygraphbuilder.cpp
@@ -425,7 +425,7 @@ namespace io {
                 nodeDepth = *kaTag;
             } else if (auto rdTag = gfa::getTag<int64_t>("rd", record.tags)) {
                 graph.m_depthTag = "rd";
-                nodeDepth = *rdTag;
+                nodeDepth = *rdTag + 1;
             } else if (auto kcTag = gfa::getTag<int64_t>("KC", record.tags)) {
                 graph.m_depthTag = "KC";
                 nodeDepth = double(*kcTag) / double(length);

diff --git a/graph/gfawriter.cpp b/graph/gfawriter.cpp
@@ -72,7 +72,7 @@ namespace gfa {
         if (depthTag == "DP" || depthTag == "dp")
             gfaSegmentLine += "\tDP:f:" + QString::number(node->getDepth()).toLatin1();
         else if (depthTag == "rd")
-            gfaSegmentLine += "\trd:i:" + QString::number(node->getDepth()).toLatin1();
+            gfaSegmentLine += "\trd:i:" + QString::number(node->getDepth()-1).toLatin1();
         else if (depthTag == "KC" || depthTag == "RC" || depthTag == "FC")
             gfaSegmentLine += "\t" + depthTag.toLatin1() + ":i:" +
                               QString::number(int(node->getDepth() * gfaSequence.length() + 0.5)).toLatin1();