We study possible effects of supersymmetry (SUSY) in CP asymmetries in non-leptonic B-d,B-s decays in a variety of SUSY flavour models considered in literature, We use both mass insertion and vertex mixing methods to calculate squark-gluino box diagrams contribution to B-d,B- s-(B) over bar(d, s) mixings. With the squark mixing parameter eta = 0.22, and with large new CP phases, it turns out that the CP asymmetries to be measured in upcoming B-factories, HERA-B and LHC-B, can be completely dominated by the SUSY contribution in almost every considered model. Discrimination between the different models can be done by comparing experimental results in different decay modes. In some models squark masses up to similar to 5 TeV can be probed through these experiments provided the SUSY contribution to B-(B) over bar mixing is at 10% level, \M-12(SUSY)/M-12(SM)\ similar to 0.1. This implies that models with heavy squarks have a fair chance to be tested in the future CP experiments before LHC.