question about emmeans post-hoc analysis

Hi list,

I have a question about post-hoc analyses using emmeans function. The model
I am using has four independent variables, pre-test scores, group (4
levels), time (4 levels), and the interaction between group and time. The
dependent variable is post-test scores. I would like to know the
differences between the four groups in terms of their post-test scores at
each time point while controlling for the effects of pre-test scores. My
model is:

*mod<-lmer(posttest scores~pre-test
scores+group+time+group*time+(1|Subject),data=GJT)*

For the post-hoc analyses, I used emmeans. The R code is:

*emmeans(mod, list(pairwise~group|time), adjust="bonferroni")*

My question is: Do the pair-wise comparisons carried out through emmeans
account for the influence of pre-test scores? Because the omnibus test
shows that the effect of pre-test scores is significant, I would like to
exclude the influence of pre-test scores on pair-wise comparisons.

Thank you in advance for any advice.

Best,

Mengxia