INNER CODE UNIT · C++

len

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

      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);

      // Complement the sequence
      std::transform(seq.begin(), seq.end(), seq.begin(), complement);

      // Copy the sequence in reverse
      std::copy(seq.rbegin(), seq.rend(), std::back_inserter(bnd));
      std::copy(find2_it + 1, alt.end(), std::back_inserter(bnd));
    }
  }
  else
  {
    find_it = std::find(alt.begin(), alt.end(), ']');

    if (find_it == alt.end())
      invalid_bnd_alt();

View source record →