The lattice Boltzmann method (LBM) has two key steps: collision and streaming. In a conventional LBM, the streaming is exact, where each distribution function is perfectly shifted to the neighbor node on the uniform mesh arrangement. This advantage may curtail the applicability of the method to problems with complex geometries. To overcome this issue, a high-order meshless interpolation-based approach is proposed to handle the streaming step. Owing to its high accuracy, the radial basis function (RBF) is one of the popular methods used for interpolation. In general, RBF-based approaches suffer from some stability issues, where their stability strongly depends on the shape parameter of the RBF. In the current work, a stabilized RBF approach is used to handle the streaming. The stabilized RBF approach has a weak dependency on the shape parameter, which improves the stability of the method and reduces the dependency of the shape parameter. Both the stabilized RBF method and the streaming of the LBM are used for solving three benchmark problems. The results of the stabilized method and the perfect streaming LBM are compared with analytical solutions or published results. Excellent agreements are observed, with a little advantage for the stabilized approach. Additionally, the computational cost is compared, where a marginal difference is observed in the favor of the streaming of the LBM. In conclusion, one could report that the stabilized method is a viable alternative to the streaming of the LBM in handling both simple and complex geometries.