This paper investigates the FIR systems identification with quantized output observations and a large class of quantized inputs. The limit inferior of the regressors’ frequencies of occurrences is employed to characterize the input’s persistent excitation, under which the strong convergence and the convergence rate of the two-step estimation algorithm are given. As for the asymptotical efficiency, with a suitable selection of the weighting matrix in the algorithm, even though the limit of the product of the Cram´er-Rao (CR) lower bound and the data length does not exist as the data length goes to infinity, the estimates still can be asymptotically efficient in the sense of CR lower bound. A numerical example is given to demonstrate the effectiveness and the asymptotic efficiency of the algorithm.