case 'M': d->c->cap_mapQ = atoi(optarg); break;
case 'd': d->max_depth = atoi(optarg); break;
case 'c': d->format |= BAM_PLF_CNS; break;
case 'M': d->c->cap_mapQ = atoi(optarg); break;
case 'd': d->max_depth = atoi(optarg); break;
case 'c': d->format |= BAM_PLF_CNS; break;