Let $${\mathbb K} $$ denote a field. Let it d denote a nonnegative integer and consider a sequence p=( $$\theta_i, \theta^*_i,i=0...d; \varphi_j, \phi_j,j=1...{\it d})$$ consisting of scalars taken from $${\mathbb K} $$ . We call p a parameter array whenever: (PA1) $$\theta_i \not=\theta_j, \; \theta^*_i\not=\theta^*_j$$ if $$i\not=j$, $(0 \leq i, j\leq d)$; (PA2) $ \varphi_i\not=0$, $(1 \leq i \leq d)$; (PA3) $\varphi_i = \phi_1 \sum_{h=0}^{i-1} ({\theta_h-\theta_{d-h}})/({\theta_0-\theta_d}) + (\theta^*_i-\theta^*_0)(\theta_{i-1}-\theta_d)$ $(1 \leq i \leq d)$; (PA4) $\phi_i = \varphi_1 \sum_{h=0}^{i-1} ({\theta_h-\theta_{d-h}})/({\theta_0-\theta_d}) + (\theta^*_i-\theta^*_0)(\theta_{d-i+1}-\theta_0)$ $(1 \leq i \leq d)$; (PA5) $(\theta_{i-2}-\theta_{i+1})(\theta_{i-1}-\theta_i)^{-1}$, $(\theta^*_{i-2}-\theta^*_{i+1})(\theta^*_{i-1}-\theta^*_i)^{-1}$$ are equal and independent of i for $$2 \leq i \leq d-1$$ . In Terwilliger, J. Terwilliger, Linear Algebra Appl., Vol. 330(2001) p. 155 we showed the parameter arrays are in bijection with the isomorphism classes of Leonard systems. Using this bijection we obtain the following two characterizations of parameter arrays. Assume p satisfies PA1 and PA2. Let A, B,A^*, B^* denote the matrices in $${Mat}_{{\it d}+1}$$ ( $${\mathbb K} $$ ) which have entries A ii =θ i , B ii =θ d-i , A * ii =θ* i , B * ii =θ* i (0 ≤ i ≤ d), A i,i-1=1, B i,i-1=1, A * i-1,i =φ i , B * i-1,i =ϕ i (1 ≤ i ≤ d), and all other entries 0. We show the following are equivalent: (i) p satisfies PA3–PA5; (ii) there exists an invertible G ∈Mat d+1( $${\mathbb K} $$ ) such that G −1 AG=B and G −1 A * G=B *; (iii) for 0 ≤ i ≤ d the polynomial $$ \sum_{n=0}^i \frac{ (\lambda-\theta_0) (\lambda-\theta_1) \cdots (\lambda-\theta_{n-1}) (\theta^*_i-\theta^*_0) (\theta^*_i-\theta^*_1) \cdots (\theta^*_i-\theta^*_{n-1}) } {\varphi_1\varphi_2\cdots \varphi_n}$$ is a scalar multiple of the polynomial $$\sum_{n=0}^i \frac{ (\lambda-\theta_d) (\lambda-\theta_{d-1}) \cdots (\lambda-\theta_{d-n+1}) (\theta^*_i-\theta^*_0) (\theta^*_i-\theta^*_1) \cdots (\theta^*_i-\theta^*_{n-1}) } {\phi_1\phi_2\cdots \phi_n}.$$ We display all the parameter arrays in parametric form. For each array we compute the above polynomials. The resulting polynomials form a class consisting of the q-Racah, q-Hahn, dual q-Hahn, q-Krawtchouk, dual q-Krawtchouk, quantum q-Krawtchouk, affine q-Krawtchouk, Racah, Hahn, dual-Hahn, Krawtchouk, Bannai/Ito, and Orphan polynomials. The Bannai/Ito polynomials can be obtained from the q-Racah polynomials by letting q tend to −1. The Orphan polynomials have maximal degree 3 and exist for ( $${\mathbb K} $$ )=2 only. For each of the polynomials listed above we give the orthogonality, 3-term recurrence, and difference equation in terms of the parameter array.