 #define EM 0.33333333333
 #define EI 0.25
 
-#define MAX_BAND_WIDTH 7
-#define MAX_BAND_WIDTH2 (MAX_BAND_WIDTH * 2 + 1)
-#define MIN_BASE_QUAL 4
+#define MIN_BAND_WIDTH 7
+#define MIN_BAND_WIDTH2 (MIN_BAND_WIDTH * 2 + 1)
 
 // all bases with q < minBaseQual are up'd to this value
 #define MIN_BASE_QUAL 4

//        fprintf(stderr, "read %d: _iqual = { % 3d % 3d % 3d % 3d % 3d ... % 3d % 3d % 3d % 3d % 3d }\n", read_index,
//                inputQualities[0], inputQualities[1], inputQualities[2], inputQualities[3], inputQualities[4],
//                inputQualities[queryLen - 5], inputQualities[queryLen - 4], inputQualities[queryLen - 3], inputQualities[queryLen - 2], inputQualities[queryLen - 1]);
//        fprintf(stderr, "read %d: c->bw = %d, bw = %d, l_ref = %d, l_query = %d\n", read_index, MAX_BAND_WIDTH, bandWidth, referenceLength, queryLen);

        /*** forward ***/
        // f[0]

        else
            bandWidth = queryLen;

        if (MIN_BAND_WIDTH < abs(referenceLength - queryLen))
        {
            bandWidth = abs(referenceLength - queryLen) + 3;
        }

        if (bandWidth > MIN_BAND_WIDTH)
            bandWidth = MIN_BAND_WIDTH;

        if (bandWidth < abs(referenceLength - queryLen))
        {
