  Predictions are made for the jet substructure of one-jet events produced in
electron-proton neutral current deep inelastic scattering at the future
Electron-Ion Collider for exchanged four-momentum squared, $Q^2 > 125$ GeV$^2$.
Data are simulated using Monte Carlo event generators PYTHIA 8.304 and RAPGAP
3.308 at the center of mass energies $\sqrt{s}$ = 63.2, 104.9 and 141 GeV. Jets
and subjets are produced by longitudinally invariant $k_T$ cluster algorithm.
The subjet multiplicity distributions and differential jet shapes are measured
for different jet sizes and varying jet-resolution parameter. A comparison is
presented between the $k_T$ and anti-$k_T$ cluster algorithms for the study of
jets and subjets using the simulated data and the HERA data.
