To make this question easier, I'll put the center of the incircle in the origin. Any integral solution can be easily transferred to the location you indicated. Let $M,N,P$ be the corners of the triangle, as used in your question. Let the coordinates of $M$ be $(M_x,M_y)\in\mathbb Z^2$. One important quantity in your scenario is
$$d:=M_x^2+M_y^2-25$$
As you can see, as long as $M$ is outside the circle, that quantity will be positive. Its square root can be used to describe the tangents from $M$ to the circle:
\begin{align*}
t_1 &= \left\{(x,y)\;\middle\vert\;
(5M_x+\sqrt dM_y)x + (5M_y-\sqrt dM_x)y = 5\left(M_x^2+M_y^2\right) \right\} \\
t_2 &= \left\{(x,y)\;\middle\vert\;
(5M_x-\sqrt dM_y)x + (5M_y+\sqrt dM_x)y = 5\left(M_x^2+M_y^2\right) \right\}
\end{align*}
The way I obtained these tangents is by performing the dual of a conic-line intersection, as described in section 11.3 of Perspectives on Projective Geometry by Richter-Gebert.
If these tangents are to cross any points with integer coordinates, their slopes have to be rational (or $\infty$ for vertical lines). These slopes are
$$-\frac{5M_x\pm\sqrt dM_y}{5M_y\mp\sqrt dM_x}$$
There are two conceivable ways these slopes might be rational: either because $\sqrt d$ is rational, or because numerator and denominator are rational multiples of the same irrational number. If one treats $\mathbb Q[\sqrt d]$ as a two-dimensional $\mathbb Q$-vectorspace, then this translates into the requirement that the two vectors have to be linearily dependent, which here means
$$\det\begin{pmatrix}
5M_x & \pm M_y \\
5M_y & \mp M_x
\end{pmatrix}=\mp5\left(M_x^2+M_y^2\right)=0$$
So the only other way would mean $M$ is the origin, which is of course impossible. So we are left with the requirement that $\sqrt d$ be rational, so $d$ must be a square number.
Now one possible way to tackle this problem is iterating over all integral coordinates in a given range, and checking whether the resulting $d$ is square. One can build up a collection of possible corner points, and using the coordinates of the tangents, one can see which corners could belong to the same triangle. The $t_1$ of one point must match the $t_2$ of the other edge, but scalar multiples of the coefficients describe the same line. I've implemented such an enumeration. Many results contain $5$ as one of the coordinates of one of the points. But not all do, so one of the more interesting results is
\begin{align*}
M&=(-11, -2) & N&=( 1, 7) & P&=( 13, -9)
\end{align*}
After moving $C$ back to where you defined it, this becomes
\begin{align*}
M&=(-10, 0) & N&=( 2, 9) & P&=( 14, -7)
\end{align*}
However, the solution above has aright angle at $N$, so as your question asks for a slution without a right angle, this one is not acceptable. Thanks to @coffeemath for making me aware of this condition. It seems that at least if the absolute value of the coordinates is restricted to no more than $1000$, then all solution without a right angle neccessarily have either a horizontal or a vertical tangent. So here is another go, again first in the untranslated situation:
\begin{align*}
M&=(-11, -10) & N&=( -3, 5) & P&=( 25, 5)
\end{align*}
After moving $C$ back to where you defined it, this becomes
\begin{align*}
M&=(-10, -8) & N&=( -2, 7) & P&=( 26, 7)
\end{align*}
I also adapted my program to filter out the right-angled triangles automatically. There are still many possible solutions remaining, but none with all coordinates less than $25$ in absolute value, in the coordinate system where the circle center is the origin.
If you are only given one set of sides and the two points, and don't want to say it is square (which would give you the other set of sides), you don't have enough information to solve the problem. The rectangle could be as you have drawn the left picture, or $L1,L2$ (which are equal) could be shorter and the rectangle more tilted.
If you are given both sides as on the right, you have enough information. Do you need the coordinates of the other corners (which are not the $(X1,Y2)$ you label) or is that just a means to get the angle of the bottom side? You can evaluate the angle of he diagonal $(X1,Y1)$ to $(X2,Y2)$ as Atan2(Y1-Y2,X1-X2). Note that this uses the usual math convention of $+x$ to the right, $+y$ up-you may need to reorient it for your coordinates. The angle from the diagonal to the lower side is $\arctan(\frac {70}{60})$ and the angle of the lower side is the sum of these.
Best Answer
If $\gamma = 0$, vector $\vec{u} = ( w, 0 )$ and vector $\vec{v} = ( 0, b )$.
Rotating a 2D vector $(x , y)$ counterclockwise yields $$\begin{cases} x^, = x \cos\gamma - y \sin\gamma \\ y^, = x \sin\gamma + y \cos\gamma \end{cases} \tag{1}\label{na1}$$ therefore $$\vec{u} = ( w \cos\gamma ,\, w \sin\gamma ) \tag{2}\label{na2}$$ and $$\vec{v} = ( -b \sin\gamma ,\, b \cos\gamma ) \tag{3}\label{na3}$$
If we know $\gamma$, $w$, $b$, and $\vec{z_c}$, then the four vertices of the rotated rectangle are $$\begin{cases} \vec{z_1} = \vec{z_c} - \vec{u} + \vec{v} \\ \vec{z_2} = \vec{z_c} + \vec{u} + \vec{v} \\ \vec{z_3} = \vec{z_c} + \vec{u} - \vec{v} \\ \vec{z_4} = \vec{z_c} - \vec{u} - \vec{v} \end{cases} \tag{4}\label{na4}$$ or equivalently, $$\begin{cases} x_1 = x_c - w \cos\gamma - b \sin\gamma \\ x_2 = x_c + w \cos\gamma - b \sin\gamma \\ x_3 = x_c + w \cos\gamma + b \sin\gamma \\ x_4 = x_c - w \cos\gamma + b \sin\gamma \end{cases}, \qquad \begin{cases} y_1 = y_c - w \sin\gamma + b \cos\gamma \\ y_2 = y_c + w \sin\gamma + b \cos\gamma \\ y_3 = y_c + w \sin\gamma - b \cos\gamma \\ y_4 = y_c - w \sin\gamma - b \cos\gamma \end{cases} \tag{5}\label{na5}$$