INNER CODE UNIT · C++
padding_size
DecodeGenetics/graphtyper · src/graph/constructor.cpp:792
std::size_t padding_size = EXTRA_SEQUENCE_LENGTH - dup.size();
/// Read the beginning of the duplicated sequence
std::copy(dup.begin(), dup.end(), std::back_inserter(dup_begin));
read_reference_seq(dup_begin, fasta_index, chrom_idx, var.pos + 1, padding_size);
append_sv_tag_to_node(dup_begin);
sv.model = "BREAKPOINT1";
sv.related_sv = static_cast<int>(graph.SVs.size()) + 1; // Next SV
graph.SVs.push_back(sv);
/// Read the end of the duplicated sequence
// Make sure we don't read before the chromosome starts!
padding_size = std::min(padding_size, static_cast<std::size_t>(var2.pos));
append_sv_tag_to_node(dup_end);
read_reference_seq(dup_end, fasta_index, chrom_idx, var2.pos - padding_size + 1, padding_size);
std::move(dup.begin(), dup.end(), std::back_inserter(dup_end));
sv.related_sv = static_cast<int>(graph.SVs.size()) - 1;