Exact four-body equations of Alt, Grassberger and Sandhas are solved for neutron-3 He and proton-3 H scattering in the energy regime above the four-nucleon breakup threshold. Cross sections and spin observables for elastic, transfer, charge-exchange, and breakup reactions are calculated using realistic nucleon-nucleon interaction models, including the one with effective many-nucleon forces due to explicit ∆-isobar excitation. The experimental data are described reasonably well with only few exceptions such as vector analyzing powers.PACS numbers: 21.45.-v, 25.10.+s, 21.30.-x, 24.70.+s Collisions and reactions are among the most important processes used to study various quantum systems ranging from ultracold atoms to nuclear and particle physics. However, reliable information about their properties can be extracted from experimental data only when an accurate theoretical description of the scattering process is available which is much more complicated to obtain as compared to bound systems. With the three-particle system already well under control the next step is the study of collisions and reactions involving four-particles.Due to its inherent complexity, rich structure of resonances, and multitude of channels the four-nucleon (4N ) system constitutes a highly challenging but also promising theoretical laboratory to test the nucleon-nucleon (N N ) force models. But for that to be possible one needs to be able to solve numerically, over a broad range of energy, the corresponding scattering equations in momentum or coordinate space. Work on this problem evolved slowly over the years but took a fast leap forward in the last ten years.Four-nucleon scattering results with realistic force models emerged first through coordinate space calculations but were limited to single channel n-3 H and p-3 He [1-3] and p-3 H [4] reactions below inelastic threshold. In that region the 4N system exhibits a rich structure of resonances [5] in different partial waves that have been well identified in the literature and whose understanding in terms of the underlying force models constitutes a major yet unresolved challenge for theory. More recent results show that adding a three-nucleon (3N ) force does not necessarily improve the agreement with the experimental data [3,6,7], unless particular 3N force models are used [8]. As in the three-nucleon system, complex scaling methods are now being used to calculate single channel reactions above breakup threshold, however, with semirealistic S-wave interactions so far [9].Given that the treatment of the Coulomb interaction between protons became possible in momentum-space calculations by using the method of screening and renor- and their respective time reversal. In this energy range calculations were done using realaxis integration with subtraction method, spline interpolation, and Padé summation [15] for matrix elements of the transition operators [12]. This approach was extended to allow for an explicit ∆-isobar excitation yielding effective 3N and 4N forces [16]. Unlike coordinate s...