  Static and dynamic aspects of the fission process of $^{226}$Th are analyzed
in a self-consistent framework based on relativistic energy density
functionals. Constrained relativistic mean-field (RMF) calculations in the
collective space of axially symmetric quadrupole and octupole deformations,
based on the energy density functional PC-PK1 and a $\delta$-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\"odinger 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 $^{226}$Th. 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 pre-neutron
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.
