We study the first curvature correction to the string amplitude of four Kaluza--Klein (KK) modes on
AdS3×S3×M4, with
M4=K3 or
T4, in type IIB string theory, which is holographically dual to the four--point correlator
⟨Op1Op2Op3Op4⟩ of certain half--BPS operators in the boundary D1--D5 CFT. The result takes the form of an integral over the Riemann sphere, analogous to the flat-space Virasoro--Shapiro amplitude, but with insertions of single-valued multiple polylogarithms of weight three. Our results are obtained in two steps. First, we derive the
AdS3×S3 Virasoro--Shapiro amplitude in the special case
⟨OpOpO1O1⟩, by matching the CFT block expansion with an ansatz based on single-valued multiple polylogarithms. We then employ the
AdS×S Mellin formalism to generalize the result to the general case of four arbitrary KK modes
⟨Op1Op2Op3Op4⟩. Our analysis yields an infinite set of results for operator anomalous dimensions and OPE data in D1--D5 CFT at strong coupling. In particular, the resulting scaling dimensions of certain operators are shown to be consistent with classical string theory computations.