This paper is concerned with a non-homogeneous in space and non-local in time random walk model for anomalous subdiffusive transport of cells. Starting with a Markov model involving a structured probability density function, we derive the non-local in time master equation and fractional equation for the probability of cell position. We derive the fractional Fokker-Planck equation for the density of cells and apply this equation to the anomalous chemotaxis problem. We show the structural instability of fractional subdiffusive equation with respect to the partial variations of anomalous exponent. We find the criteria under which the anomalous aggregation of cells takes place in the semi-infinite domain.