da_out1 <- differential_abundance_wrapper(counts, md = md, 
  formula = formula_glmer_binomial1, K = K)

apply(da_out1$adjp < FDR_cutoff, 2, table)

adjp_BCRXLvsRef
FALSE               5
TRUE                3

da_out2 <- differential_abundance_wrapper(counts, md = md, 
  formula = formula_glmer_binomial2, K = K)

apply(da_out2$adjp < FDR_cutoff, 2, table)

adjp_BCRXLvsRef
FALSE               2
TRUE                6