INNER CODE UNIT · C++

end

soedinglab/plass · lib/flash/combine_reads.cpp:288

    end = read_1->seq_len - min_overlap + 1;
    for (int i = start; i < end; i++) {
        unsigned num_mismatches;
        unsigned mismatch_qual_total;
        int overlap_len = read_1->seq_len - i;

        compute_mismatch_stats(read_1->seq + i,
                               read_2->seq,
                               read_1->qual + i,
                               read_2->qual,
                               haveN,
                               &overlap_len,
                               &num_mismatches,
                               &mismatch_qual_total);

        if (overlap_len >= min_overlap) {
            float score_len = (float)min(overlap_len, max_overlap);
            float qual_score = mismatch_qual_total / score_len;

View source record →