Wrong result of evalf in SymPy but only for some values of the precision

Question

I think I have encountered a bug of Sympy\'s evalf() method with substitutions passed.

By accident, I found an expression that evaluates to a wrong value