Let $k>2$ and $\epsilon>0$ be fixed real numbers. Let $x$ be a large parameter. It suffices to show that there exists a prime pair $(p,q)$ such that $x<p<2x$ and $p^k<q<p^k(1+\epsilon)$ and $\left(\frac{p}{q}\right)=1$. Note that for primes $p\equiv 1\pmod{4}$, the last condition is equivalent to $\left(\frac{q}{p}\right)=1$. Let us now call a prime $p\equiv 1\pmod{4}$ bad when there is no prime $q$ such that $p^k<q<p^k(1+\epsilon)$ and $\left(\frac{q}{p}\right)=1$. It suffices to show that the number of bad primes $p\in(x,2x)$ is $\ll x/(\log x)^2$.

For a prime $p$, let $\chi_p$ (resp. $1_p$) be the quadratic (resp. trivial) Dirichlet character mod $p$. If $p$ is bad, then with the usual notations we have
$$\sum_{\substack{p^k<n\leq p^k(1+\epsilon)\\\left(\frac{q}{p}\right)=1}}\Lambda(n)\ll p^{k/2}=o(p^k),$$
hence also
$$\left(\psi(p^k(1+\epsilon),1_p)-\psi(p^k,1_p)\right)+\left(\psi(p^k(1+\epsilon),\chi_p)-\psi(p^k,\chi_p)\right)=o(p^k).$$
The first difference is asymptotically $\epsilon p^k$, hence for a bad prime $p$ we have
$$\left|\psi(p^k(1+\epsilon),\chi_p)-\psi(p^k,\chi_p)\right|\gg p^k,$$
hence also
$$\max_{y<(3x)^k}\left|\psi(y,\chi_p)\right|\gg p^k.$$
On the other hand, by the usual proof of the Bombieri-Vinogradov theorem (see e.g. Chapter 24 in Huxley: The distribution of prime numbers) we can see that that
$$\sum_{x<p<2x}\max_{y<(3x)^k}\left|\psi(y,\chi_p)\right|\ll\frac{x^{k+1}}{(\log x)^2},$$
using that $x$ is large and $k>2$. As we observed above, for a bad prime $p\in(x,2x)$ the corresponding term on the left hand side is $\gg x^k$, whence the number of such primes is $\ll x/(\log x)^2$. The proof is complete.