In this article, we explore the famous Selkov-Schnakenberg (SS) system of coupled nonlinear partial differential equations (PDEs) for Lie symmetry analysis (LSA), Self Adjointness (SA) and Conservation laws (CLs). Moreover miscellaneous soliton solutions like dark, bright, periodic, rational, Jacobian elliptic function (JEF), Weierstrass elliptic and hyperbolic solutions of SS system will be achieved by a well-known technique called Sub-ODE. All these results are displayed graphically by 3D, 2D and contour plots.