INNER CODE UNIT · C
base
dpryan79/MethylDackel · MBias.c:180
base = *(seq+pos-localPos);
for(i=0; i<n_plp; i++) {
if(plp[0][i].is_del) continue;
if(plp[0][i].is_refskip) continue;
if(config->bed) if(!readStrandOverlapsBED(plp[0][i].b, config->bed->region[bedIdx])) continue;
strand = getStrand((plp[0]+i)->b);
if(strand & 1) {
if(base != 'C' && base != 'c') continue;
} else {
if(base != 'G' && base != 'g') continue;
}
rv = updateMetrics(config, plp[0]+i);
if(rv != 0) {
if((plp[0]+i)->qpos >= meths[strand-1]->m)
meths[strand-1] = growStrandMeth(meths[strand-1], (plp[0]+i)->qpos);
if(rv < 0) {
if((plp[0]+i)->b->core.flag & BAM_FREAD2) {
assert((meths[strand-1]->unmeth2[(plp[0]+i)->qpos]) < 0xFFFFFFFF);