INNER CODE UNIT · C++
start_it
DecodeGenetics/graphtyper · src/graph/genomic_region.cpp:214
auto start_it = reference.begin() + (pos - GenomicRegion::begin);
std::size_t const ref_len =
std::min(static_cast<std::size_t>(std::distance(start_it, reference.end())), record.ref.seq.size());
std::vector<char> const reference_ref(start_it, start_it + ref_len);
if (reference_ref.size() > 0 && !prefix_match(record.ref.seq, reference_ref))
{
print_log(log_severity::warning,
__HERE__,
" Record @ position ",
pos,
" (REF = ",
std::string(record.ref.seq.begin(), record.ref.seq.end()),
')',
" did not match the reference genome (",
std::string(reference_ref.begin(), reference_ref.end()),
")");