Consider the following toy example:

sol = NDSolve[{ x'[t] == {-1, -2} x[t], x[0] == {2, 2} }, x, {t, 0, 4} ]; Plot[x[t] /. sol, {t, 0, 4}, PlotRange -> All] 

As seen above, both plots are drawn in the same colour.

The usual trick of adding Evaluate does not work here, I guess because x[t] /. sol evaluates to a list with a single element:

and it's the InterpolatingFunction that later evaluates to a list.

Manually specifying the PlotStyle also doesn't seem to work.

How can I ensure that the different plots are assigned different colours/styles?