Since the returns involve changes in the stock price over consecutive time periods, the answer to your question depends on the joint distribution of the stock price over time. Since you have only specified the marginal distribution of the stock price, you have not given enough information in your question to determine the distribution of interest. Nevertheless, I will give you a useful result here that holds in a simple case, and which can be extended more broadly to other models with some additional algebraic work.
Stationary AR(1) model: Suppose the log-price follows the stationary model of the form:
$$\ln P_t = (1-\phi) \mu + \phi \ln P_{t-1} + \varepsilon_t \quad \quad \quad \varepsilon_t \sim \text{N}( 0, (1-\phi^2) \sigma^2),$$
where $-1 < \phi < 1$ is the auto-correlation parameter of the model. In this case the stationary marginal distribution of the log-price is as you specified in your question (i.e., $\ln P_t \sim \text{N} (\mu, \sigma^2)$). Now, since $\ln P_t = \ln P_{t-1} + \ln (1+R_t)$ the density function of the latter term is:
$$\begin{equation} \begin{aligned}
p(\ln (1+R_t) = r)
&= \int \limits_{-\infty}^\infty p(\ln (1+R_t) = r | \ln P_{t-1} = s) \cdot \text{N}(s|\mu, \sigma^2) ds \\[6pt]
&= \int \limits_{-\infty}^\infty p(\ln P_t = r+s | \ln P_{t-1} = s) \cdot \text{N}(s|\mu, \sigma^2) ds \\[6pt]
&= \int \limits_{-\infty}^\infty \text{N}(r+s | \mu + \phi (s-\mu), (1-\phi^2) \sigma^2) \cdot \text{N}(s|\mu, \sigma^2) ds. \\[6pt]
\end{aligned} \end{equation}$$
To solve this integral equation we need to complete the square in the exponential terms of the normal density functions in the integrand. This would lead to a normal distribution, so that $1+R \sim \text{LogN}$. Obtaining the parameters of the distribution would require you to complete the square in the above integral, and perform the subsequent algebra. I leave this as an exercise.