The QCD transition temperature is determined using Symanzik improved gauge and stout-link improved staggered fermionic lattice simulations with physical masses for light and strange quarks. Four lattice spacings (Nt=4, 6, 8, 10) were used for continuum extrapolation, with Nt=6, 8, and 10 in the scaling region. The QCD transition is a non-singular cross-over, leading to different Tc values from different observables. The peak of the renormalized chiral susceptibility predicts Tc=151(3)(3) MeV, while the strange quark number susceptibility and Polyakov loop give higher values of 24(4) MeV and 25(4) MeV, respectively. The non-vanishing width of the peaks is also determined. These results are considered the full result for the T≠0 transition, though other lattice formulations are needed for cross-checking. The study highlights the importance of using physical masses and proper scaling to determine Tc accurately. The results show that lattice spacings larger than ~0.20 fm are not in the scaling region, leading to unreliable continuum extrapolations. The study also emphasizes the ambiguity in setting the overall scale and the need for consistent scale-setting methods. The results are consistent with the continuum limit and show that different observables yield different Tc values, reflecting the cross-over nature of the QCD transition.The QCD transition temperature is determined using Symanzik improved gauge and stout-link improved staggered fermionic lattice simulations with physical masses for light and strange quarks. Four lattice spacings (Nt=4, 6, 8, 10) were used for continuum extrapolation, with Nt=6, 8, and 10 in the scaling region. The QCD transition is a non-singular cross-over, leading to different Tc values from different observables. The peak of the renormalized chiral susceptibility predicts Tc=151(3)(3) MeV, while the strange quark number susceptibility and Polyakov loop give higher values of 24(4) MeV and 25(4) MeV, respectively. The non-vanishing width of the peaks is also determined. These results are considered the full result for the T≠0 transition, though other lattice formulations are needed for cross-checking. The study highlights the importance of using physical masses and proper scaling to determine Tc accurately. The results show that lattice spacings larger than ~0.20 fm are not in the scaling region, leading to unreliable continuum extrapolations. The study also emphasizes the ambiguity in setting the overall scale and the need for consistent scale-setting methods. The results are consistent with the continuum limit and show that different observables yield different Tc values, reflecting the cross-over nature of the QCD transition.