Too long for a comment, but not a full answer.
This post suggests that the norm of a convolution operator $f\mapsto f*g$ on $\mathbb R$ is $\|g\|_{L^1}$ if $g\geq 0$ (in general, it is the $L^\infty$ norm of the Fourier transform of $g$). The main intuition behind it is that, by scaling symmetry, you can think of $g$ as being vey close to a Dirac delta multiplied by $\|g\|_{L^1}$. Most likely, the same is not true here, due to the boundedness of the domain which does not allow for arbitrary scaling of functions.
One estimate I can think of on the operator norm of $T$ is the following. By Riesz-Thorin theorem, one has the bound
$$ \|T\|_{L^2\to L^2}\leq \|T\|_{L^1\to L^1}^{1/2}\|T\|_{L^\infty\to L^\infty}^{1/2}, $$
where $\|T\|_{L^p\to L^p}$ is the operator norm of the operator on $L^p(0,1)$. The norm $\|T\|_{L^1\to L^1}$ is $\|g\|_{L^1(-1/2,1/2)}$, where $g(x):=|\sin(x)|^{-\alpha}$. The same is true for the norm $\|T\|_{L^\infty\to L^\infty}$. So, by the above bound, one has
$$ \|T\|_{L^2\to L^2}\leq \int_{-1/2}^{1/2}|\sin(x)|^{-\alpha}dx, $$
which is quite close to what you have done to show the boundedness of $T$, if not the same. The above estimate is likely not optimal, and I don't know if there is an easy way of finding a better estimate or compute the exact operator norm.
Another way of estimating the norm $\|T\|_{L^2\to L^2}$ is by using Schur's test, but the standard one with $q=p\equiv 1$ gives the exact same bound I wrote above.
Edit 1. This operator is an integral Hilbert-Schmidt operator if $\alpha<1/2$. For this kind of operator, one also has the bound
$$ \|T\|_{L^2\to L^2}\leq \int_{[0,1]\times[0,1]}|\sin(x-y)|^{-2\alpha}dxdy, $$
see e.g. this pdf, Theorem 8.2.1.
The first estimate I gave sounds better, but maybe the second one is better for some values of $\alpha$.
Edit 2. The operator $T$ is a self-adjoint, non-negative compact operator on $L^2(0,1)$ for any $0<\alpha<1$. By general facts, this implies that all the eigenvalues of $T$ are non-negative real numbers, and $\|T\|_{L^2\to L^2}$ coincides with the largest eigenvalue of $T$. Therefore, if one has a way of finding or estimating the eigenvalues of $T$, this in principle gives an estimate on the operator norm. I don't see how to find the eigenvalues at the moment, but this has better hope of giving the sharp constant.
Edit 3. I don't have time to write a proof of $\|T\|_{L^1\to L^1}=\|g\|_{L^1(-1/2,1/2)}$. I can give some heuristics though. For a function $f$ with $\|f\|_{L^1(0,1)}=1$, it is evident that in order to maximize $\|Tf\|_{L^1(0,1)}$, we can assume $f\geq 0$ a.e.. Now think of $f$ as an infinite linear combinations of Dirac deltas $\delta_{x_0}$ centered at $x_0\in [0,1]$. The operator $T$ applied to $\delta_{x_0}$ formally gives
$$ T\delta_{x_0}(x)=|\sin(x-x_0)|^{-\alpha} $$
The norm $\|T\delta_{x_0}\|_{L^1(0,1)}$ is maximized by choosing $x_0=1/2$, and in that case we have
$$ \|T\delta_{1/2}\|_{L^1(0,1)}=\|g\|_{L^1(-1/2,1/2)} $$
by a simple change of variables.
Coming back to $f$, since we want to maximize $\|Tf\|_{L^1(0,1)}$, the optimal way to do that is to choose $f$ in whose linear expansion in terms of Deltas, $\delta_{1/2}$ gives the main contribution (that is, $f$ approximates $\delta_{1/2}$ in some sense). This is because all the contributions coming from $\delta_{x_0}$ with $x_0$ far from $1/2$ will give a sub-optimal contribution to $\|Tf\|_{L^1}$. Then, by choosing $f_n(x)=\chi_{[1/2-1/n,1/2+1/n]}$, one verifies that
$$ \|Tf_n\|_{L^1(0,1)}\to \|g\|_{L^1(-1/2,1/2)} $$
as $n\to\infty$. This suggests that my claim is correct.