INNER CODE UNIT · C
strand
dpryan79/MethylDackel · MBias.c:185
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);
meths[strand-1]->unmeth2[(plp[0]+i)->qpos]++;
} else {
assert((meths[strand-1]->unmeth1[(plp[0]+i)->qpos]) < 0xFFFFFFFF);
meths[strand-1]->unmeth1[(plp[0]+i)->qpos]++;
}