INNER CODE UNIT · C
rb
dpryan79/MethylDackel · common.c:145
rb = bounds[4*strand+3];
} else {
lb = bounds[4*strand];
rb = bounds[4*strand+1];
}
lb = (lb<b->core.l_qseq) ? lb : b->core.l_qseq;
//trim on the left
if(lb) {
for(i=0; i<lb; i++) {
qual[i] = 0;
if(i&1) seq[i>>1] |= 0xf;
else seq[i>>1] |= 0xf0;
}
}
//trim on the right