                 reference_window_clipped.x, reference_window_clipped.y);
 
     uint16 err = ctx.num_errors[read_index];
-    fprintf(stderr, "    number of errors            = [ % 3d ]\n", err);
+    fprintf(stderr, "    errors in clipped region    = [ % 3d ]\n", err);
 
     fprintf(stderr, "\n");
 }

                    ctx.cigar.is_insertion[idx.read_start + current_bp_idx + batch.cigars[event].len - 1] = 1;
                }

                num_errors++;
                break;

            case cigar_event::D:

        const uint32 align_offset = batch.alignment_start[read_index];
        const uint32 reference_alignment_start = seq_base + align_offset;

        const auto& read_window_clipped = ctx.cigar.read_window_clipped[read_index];

        uint16 current_bp_idx = 0;
        uint16 num_errors = 0;

