INNER CODE UNIT · C++

combined_qual

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

    combined_qual = combined_read->qual;

    combined_read->seq_len = combined_seq_len;
    combined_read->qual_len = combined_seq_len;

    /* Copy the beginning of read 1 (not in the overlapped region).  */
    while (overlap_begin--) {
        *combined_seq++ = *seq_1++;
        *combined_qual++ = *qual_1++;
    }

    /* Copy the overlapped region.  */
    while (overlap_len--) {
        if (*seq_1 == *seq_2) {
            /* Same base in both reads.  Take the higher quality
             * value. */
            *combined_seq = *seq_1;
            *combined_qual = max(*qual_1, *qual_2);

View source record →