“…To our knowledge, the Bethe ansatz solutions were known only for some QES difference equations in connection with U q (sl(2)) [17], and for exactly solvable difference equations [18,19,20] of the elliptic Rusijsenaars-Schneider system. Here we apply the Bethe ansatz formulation to several explicit examples of QES difference equations [10,11] as deformation of exactly solvable 'discrete' quantum mechanics [2,3], which are difference analogues of the well-known quasi exactly solvable systems, the harmonic oscillator (with/without the centrifugal potential) deformed by a sextic potential and the 1/ sin 2 x potential deformed by a cos 2x potential. As will be shown explicitly in the main text, these Bethe ansatz equations (2.11), (3.11), (3.15), (4.9), (4.17), and (5.5) can be considered as deformations of the equations (2.16), (4.14), (5.9) determining the roots of the corresponding (q-)Askey scheme of hypergeometric orthogonal polynomials [9] (the MeixnerPollaczek, continuous Hahn, continuous dual Hahn, Wilson and Askey-Wilson polynomials and their restrictions) which constitute the eigenfunctions of the undeformed exactly solvable quantum systems.…”