  We develop an algorithm for recursively constructing orthogonal multiplet
bases for the color space of QCD, for any order of partons and any Nc. This
recipe is then applied for explicitly constructing some of these bases. Using
the bases, a corresponding set of Wigner 6j coefficients are calculated. The
Wigner coefficients offer a method of using multiplet bases without resorting
to the explicit expressions of the basis vectors, which lead to a significant
speed-up compared to other methods of treating full color structure.
