A family of non-Hermitian real and tridiagonal-matrix candidates H (N ) for a hiddenly Hermitian (a.k.a. quasi-Hermitian) quantum Hamiltonian is proposed and studied. Fairly weak assumptions are imposed upon the unperturbed matrix (the square-well-simulating spectrum of H (N ) 0is not assumed equidistant) as well as upon its maximally non-Hermitian N−parametric antisymmetric-matrix perturbations (matrix W (N ) (λ) is not even required to be PT −symmetric). In spite of that, the "physical" parametric domain D [N ] is (constructively!) shown to exist, guaranteeing that in its interior the spectrum remains real and non-degenerate, rendering the quantum evolution unitary. Among the non-Hermitian degeneracies occurring at the boundary ∂D [N ] of the domain of stability our main attention is paid to their extreme version corresponding to the Kato's exceptional point of order N (EPN). The localization of the EPNs and, in their vicinity, of the quantum-phase-transition boundaries ∂D [N ] is found feasible, at the not too large N, using computer-assisted symbolic manipulations including, in particular, the Gröbner basis elimination and the high-precision arithmetics.