Fix in edge subseq extraction

AntonBankevich · Nov 25, 2024 · d329d64 · d329d64
1 parent 41c4cd4
commit d329d64
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/projects/assembly_graph/assembly_graph_base.hpp b/src/projects/assembly_graph/assembly_graph_base.hpp
@@ -231,7 +231,7 @@ namespace ag {
         const Sequence &truncSeq() const { return seq; }
         const Sequence &getCode() const {return edge_code;}
         void resetEdgeCode() {edge_code = seq.Subseq(0, std::min<size_t>(1, seq.size()));}
-        Sequence kmerSeq(size_t pos) const {return fullSubseq(pos, pos);}
+        Sequence kmerSeq(size_t pos) const {return fullSubseq(pos, pos + getStartSize());}
 
 
         const Vertex &getFinish() const {return *finish;}
@@ -310,9 +310,9 @@ namespace ag {
             VERIFY(from <= start->size() + truncSize());
             VERIFY(to <= start->size() + truncSize());
             if (from >= start->size()) {
-                return truncSeq().Subseq(from - start->size(), to);
+                return truncSeq().Subseq(from - start->size(), to - start ->size());
             } else {
-                return getStart().getSeq().Subseq(from) + truncSeq().Subseq(0, to);
+                return getStart().getSeq().Subseq(from) + truncSeq().Subseq(0, to - start->size());
             }
         }