The flexibility in supporting heterogeneous services with vastly different technical requirements is one of the distinguishing characteristics of the fifth generation (5G) communication systems and beyond. One viable solution is to divide the system bandwidth into several bandwidth parts (BWPs), each having a distinct numerology optimized for a particular service. However, multiplexing of mixed numerologies over a unified physical infrastructure comes at the cost of induced interference. In this paper, we develop an analytical system model for inter-numerology interference (InterNI) analysis in orthogonal frequency-division multiplexing (OFDM) systems with and without filter processing in the presence of mixed numerologies. With the analytical model, the level of InterNI is quantified by the developed analytical metric, which is expressed as a function of several system parameters. This leads to an analysis and evaluation of these parameters for meeting a given distortion target. Moreover, a case study on power allocation utilizing the derived analysis is presented, where an optimization problem of maximizing the sum rate is formulated, and a solution is also provided. It is also demonstrated that a filtered-OFDM system better accommodates the coexistence of mixed numerologies. The proposed model provides an accurate analytical guidance for the multi-service design in 5G and beyond systems.