Studies on jet substructure have evolved significantly in recent years. Jet substructure is essentially determined by QCD radiations and non-perturbative effects. Predictions of jet substructure are usually different among Monte Carlo event generators, and are governed by the parton shower algorithm implemented. For leading logarithmic parton shower, even though one of the core variables is the evolution variable, its choice is not unique. We examine evolution variable dependence of the jet substructure by developing a parton shower generator that interpolates between different evolution variables using a parameter α. Jet shape variables and associated jet rates for quark and gluon jets are used to demonstrate the α-dependence of the jet substructure. We find angular ordered shower predicts wider jets, while relative transverse momentum (pT) ordered shower predicts narrower jets. This is qualitatively in agreement with the missing phase space of pT ordered showers. Such difference can be reduced by tuning other parameters of the showering algorithm, especially in the low energy region, while the difference tends to increase for high energy jets.