  In the previous article Eur. Phys. J. Plus 136, 590 (2021) (arXiv:
2110.15396) we have proposed a reconstruction method for the kinetic gravity
braiding theory in the framework of the flat Friedman-Robertson-Walker
spacetime. Here we develop this method in the Bianchi I spacetime model for a
subclass of the Horndeski theory: $G_5 \sim \phi $, $G_2(X)\neq0$,
$G_3(X)\neq0$. The Hubble parameter $H(t)$ and the canonical kinetic term
$X(t)$ are set a priori. The choice of the function $X(t)$ determines the
anisotropic properties of the Universe. This makes it possible to provide
believable anisotropy. The presented method allows for a realistic model of the
Universe to simply reconstruct some scalar field theory. Reconstruction example
is given for anisotropic model of a post-inflationary transition to the
radiation-dominated phase. The model is investigated for the absence ghosts and
Laplacian instabilities.
