A Monte Carlo technique was used to simulate primary recrystallization in IF steels. In order to consider anisotropic properties of grain boundary energy and grain boundary mobility, functions of boundary misorientation and coincident site lattice (CSL) were introduced. In order to consider the energy stored within grain, a parameter Q describing the pattern quality in Electron Backscatter Diffraction (EBSD) measurement was used. The EBSD measurement result of local area was analyzed to examine the possibility of the model to simulate the evolution of microstructure and texture during recrystallization. The nucleation sites of recrystallization were high angle grain boundaries (HAGBs) and shear bands. The main texture components of simulation were quite similar with those obtained experimentally.