INNER CODE UNIT · C++

is_similar

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

      auto is_similar = [](std::vector<char> const & seq1, std::vector<char> const & seq2) -> bool
      {
        seqan::Dna5String seqan_seq1;
        seqan::Dna5String seqan_seq2;

        std::size_t constexpr MAX_SIZE = 1000;

        // At most align MAX_SIZE bp
        if (seq1.size() > MAX_SIZE && seq2.size() > MAX_SIZE)
        {
          std::string tmp1(begin(seq1), begin(seq1) + MAX_SIZE);
          seqan_seq1 = tmp1.c_str();

          std::string tmp2(begin(seq2), begin(seq2) + MAX_SIZE);
          seqan_seq2 = tmp2.c_str();
        }
        else
        {

View source record →