INNER CODE UNIT · C++
len
DecodeGenetics/graphtyper · src/graph/constructor.cpp:401
auto const len = EXTRA_SEQUENCE_LENGTH - bnd.size() + 1;
read_reference_seq(bnd, fasta_index, chrom_idx_2, pos, len); // Read mate locus
append_sv_tag_to_node(bnd); // Put SV tag
}
else
{
// Case 2: S [chr:pos[NNNS => Extending reversed sequence left of chr:pos
print_log(log_severity::debug, "[", __HERE__, "] BND variant case 2 @ ", var.pos + 1);
append_sv_tag_to_node(bnd);
auto find2_it = std::find(find_it + 1, alt.end(), '[');
if (find2_it == alt.end())
invalid_bnd_alt();
auto const len = EXTRA_SEQUENCE_LENGTH - std::distance(find2_it, alt.end());
std::vector<char> seq;
read_reference_seq(seq, fasta_index, chrom_idx_2, pos - 1, len);