A unified numerically solvable framework for dispersion relations with an arbitrary number of species drifting at arbitrary directions and with Krook collision is derived in linear uniform/homogenous kinetic plasma model, which greatly extends the standard one given by Stix (1992). The purpose of this work is to provide a kinetic plasma dispersion relation tool not only in the physical model but also the numerical approach as general/powerful as possible. As a general application, we give the dispersion relations which assume the equilibrium distribution function to be bi-Maxwellian including parallel drift, two directions of perpendicular drift (i.e., drift across magnetic field), ring beam and loss-cone. Both electromagnetic and electrostatic versions are provided, with also the Darwin (a.k.a., magnetoinductive or magnetostatic) version. The species can be treated either magnetized or unmagnetized. Later, the equations are transformed to the matrix form that can be solved by using the powerful matrix algorithm (Xie and Xiao, 2016), which is the first approach can give all the important solutions of a linear kinetic plasma system without requiring initial guess for root finding. To the best of our knowledge, the present model is the most comprehensive model in literature for the distribution function constructed from Maxwellian distributions, which thus can be applied widely for studying waves and instabilities in space, astrophysics, fusion and laser plasmas. We limit the present work to non-relativistic case. Program summaryProgram Title: BOProgram Files doi:http://dx.doi.org/10.17632/cvbrftzfy5.1Licensing provisions: BSD 3-clauseProgramming language: MatlabNature of problem: The linear fluid and kinetic waves and instabilities in plasma can be described by dispersion relations. The challenges are to provide a dispersion relation as general as possible and to obtain all the solutions of it, which is the goal of BO tool. The kinetic version of BO tool provides a unified numerically solvable framework for kinetic dispersion relations, which greatly extends the standard one (Stix, 1992), with an arbitrary number of species. It contains many new features, including anisotropic temperature/loss cone/drift in arbitrary direction/ring beam/collision, unmagnetized/magnetized, electrostatic/electromagnetic/Darwin, can support parallel wave vector kz≤0, etc.Solution method: Approximating the plasma dispersion function in the kinetic dispersion relation with J-pole Padé expansion. Transforming the dispersion relation to an equivalent matrix eigenvalue problem and find all the solutions using standard matrix eigenvalue library function.Additional comments including restrictions and unusual features: Kinetic relativistic effects are not included in the present version yet.