INNER CODE UNIT · C++

padding_length

DecodeGenetics/graphtyper · src/graph/constructor.cpp:565

      long padding_length = EXTRA_SEQUENCE_LENGTH - sv.seq.size();
      assert(padding_length > 0);
      std::copy(sv.seq.begin(), sv.seq.end(), std::back_inserter(alt1));

      read_reference_seq(alt1, fasta_index, chrom_idx, var.pos + 1, padding_length);

      append_sv_tag_to_node(alt1);
      sv.related_sv = static_cast<int>(graph.SVs.size()) + 1; // Next SV is related
      sv.model = "BREAKPOINT1";
      graph.SVs.push_back(sv);
      append_sv_tag_to_node(alt2);

      read_reference_seq(alt2, fasta_index, chrom_idx, var.pos - padding_length, padding_length + 1);

      std::copy(sv.seq.begin(), sv.seq.end(), std::back_inserter(alt2));
      sv.related_sv = static_cast<int>(graph.SVs.size()) - 1; // Previous SV is related
      sv.model = "BREAKPOINT2";
      graph.SVs.push_back(std::move(sv));

View source record →

📰 Research Paper
Loading…
⏳ Fetching content…