The eigenvalue distribution of a Wilson loop operator of fixed shape undergoes a transition under scaling at infinite N. We derive a large N scaling function in a double scaling limit of the average characteristic polynomial associated with the Wilson loop operator in two dimensional QCD. We hypothesize that the transition in three and four dimensional large N QCD are also in the same universality class and provide a numerical test for our hypothesis in three dimensions.