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