scipy fit for t distribution seems broken for bi-modal data
The t-test has many assumptions. That dataset violates several of them:


Data should be sufficiently large (>30 independent points)
Data should be approximately normally distributed


Given that the assumptions are violated, you can not expect to valid results.