Static and dynamic aspects of the fission process of 226Th are analyzed in a self-consistent framework based on relativistic energy density functionals. Constrained relativistic mean-field calculations in the collective space of axially symmetric quadrupole and octupole deformations, based on the energy density functional PC-PK1 and a δ- force pairing, are performed to determine the potential energy surface of the fissioning nucleus, the scission line, the single-nucleon wave functions, energies, and occupation probabilities, as functions of deformation parameters. Induced fission dynamics is described using the time-dependent generator coordinate method in the Gaussian overlap approximation. A collective Schrödinger equation, determined entirely by the microscopic single-nucleon degrees of freedom, propagates adiabatically in time the initial wave packet built by boosting the ground-state solution of the collective Hamiltonian for 226Th. The position of the scission line and the microscopic input for the collective Hamiltonian are analyzed as functions of the strength of the pairing interaction. The effect of static pairing correlations on the preneutron emission charge yields and total kinetic energy of fission fragments is examined in comparison with available data, and the distribution of fission fragments is analyzed for different values of the initial excitation energy.