On Julia 1.11 and recent versions of AlgebraicSolving.jl, the following example

using AlgebraicSolving

R, (x, y, z) = polynomial_ring(QQ, [:x, :y, :z], internal_ordering=:degrevlex)

I = Ideal([
    x + 2 * y + 2 * z - 1,
    x^2 + 2 * y^2 + 2 * z^2 - x,
    2 * x * y + 2 * y * z - y
])

println(groebner_basis(I, eliminate=2))

unexpectedly returns QQMPolyRingElem[84*z^4 + y - 39], while msolve 0.9.1 correctly returns [84*z^4-40*z^3+z^2+z]: with -e 2 for the same input.

x,y,z
0
x+2*y+2*z-1,
x^2+2*y^2+2*z^2-x,
2*x*y+2*y*z-y

One may suspect a parsing error at the interface level, as the coefficients seem quite similar.