We analyze the light curve of 1740 flare stars to study the relationship between the magnetic feature characteristics and the identified flare activity. Coverage and stability of magnetic features are inspired by rotational modulation of light curve variations and flare activity of stars are obtained using our automated flare detection algorithm. The results show that (i) Flare time occupation ratio (or flare frequency) and total power of flares increase by increasing relative magnetic feature coverage and contrast in F-M type stars (ii) Magnetic feature stability is highly correlated with the coverage and the contrast of the magnetic structures as this is the case for the Sun (iii) Stability, coverage and contrast of the magnetic features, time occupation ratio and total power of flares increases for G, K and M-type stars by decreasing Rossby number due to the excess of produced magnetic field from dynamo procedure until reaching to saturation level.