Landau-Khalatnikov-Fradkin transformations (LKFTs) yield the gauge dependence of correlation functions within the class of linear covariant gauges. We derive the LKFT for the quark propagator and explicitly evaluate it up to the two loop level in the chiral limit. Although the number of diagrams to be evaluated is significantly larger than with the conventional computational scheme, the diagrams are simpler in nature, thereby leading to a considerably faster evaluation of the gauge dependent part than naively expected. Finally, we also resum the LKFT generated terms and compare our results with earlier work in the literature.