In recent years, lepton flavour violating (LFV) decays are one of the most trending topics to probe new physics (NP). The latest results of LHCb have motivated us to study the LFV decays through $$b\rightarrow s$$
b
→
s
transition. The branching ratios, forward backward asymmetries and longitudinal polarization fractions of $$B(B_{s})\rightarrow K^{*}(\varphi )l_{1}l_{2}$$
B
(
B
s
)
→
K
∗
(
φ
)
l
1
l
2
decays are studied in non-universal $$Z^{'}$$
Z
′
model. Here, we have structured the four-fold angular distribution of the decays in terms of transversity amplitudes. The variation of the observables in whole kinematic region shows the sensitivity of NP. The observables estimated in this work are very intriguing and might lead a new track towards NP in near future.