ArrayPlot[Function[u,MapAt[Red&,u[[2]],u[[1,2]]]]/@TuringMachine[2506,{1,{{},0}},50]]