The process $e^+e^-\rightarrow q\bar{q}$ with $q\bar{q}=s\bar{s},c\bar{c},b\bar{b},t\bar{t}$ plays a central role in the physics programs of high energy electron-positron colliders operating from $\mathcal{O}$ (100~GeV) to $\mathcal{O}$ (1~TeV) center of mass energies. Furthermore, polarized beams as available at the International Linear Collider (ILC) are an essential input for the complete measurement of the helicity amplitudes that govern the production cross section. Quarks, especially the heaviers, are likely messengers to new physics and at the same time they are ideal benchmark processes for detector optimization. All four processes call for superb primary and secondary vertex measurements, a high tracking efficiency to correctly measure the vertex charge and excellent hadron identification capabilities. Strange, charm and bottom production are already available below the $t\bar{t}$ threshold. We will show with detailed detector simulations of the International Large Detector (ILD) that production rate and the forward backward asymmetries of the the different processes can be measured at the 0.1% to 0.7% level and how systematic errors can be controlled to reach this level of accuracy. The importance of operating at different center of mass energies and the discovery potential in terms of Randall-Sundrum models with warped extra dimensions will be outlined.