Observations suggest that twisted magnetic flux tubes are ubiquitous in the Sun's atmosphere. The main aim of this work is to advance the study of axisymmetric modes of magnetic flux tubes by modeling both twisted internal and external magnetic field, when the magnetic twist is weak. In this work, we solve the derived wave equations numerically assuming twist outside the tube is inversely proportional to the distance from its boundary. We also study the case of constant magnetic twist outside the tube and solve these equations analytically. We show that the solution for a constant twist outside the tube is a good approximation to the case where the magnetic twist is proportional to $1/r$, namely the error is in all cases less than $5.4\%$.The solution is in excellent agreement with solutions to simpler models of twisted magnetic flux tubes, i.e. without external magnetic twist. It is shown that axisymmetric Alfv\'{e}n waves are naturally coupled with magnetic twist as the azimuthal component of the velocity perturbation is nonzero. We compared our theoretical results with observations and comment on what the Doppler signature of these modes is expected to be. Lastly, we argue that the character of axisymmetric waves in twisted magnetic flux tubes can lead to false positives in identifying observations with axisymmetric Alfv\'{e}n waves.