This paper considers a linear quasi-cyclic product code of two given quasi-cyclic codes of relatively prime lengths over finite fields. We give the spectral analysis of a quasi-cyclic product code in terms of the spectral analysis of the row- and the column-code. Moreover, we provide a new lower bound on the minimum Hamming distance of a given quasi-cyclic code and present a new algebraic decoding algorithm.More specifically, we prove an explicit (unreduced) basis of an l\_a l\_b-quasi-cyclic product code in terms of the generator matrix in reduced Gr{ö}bner basis with respect to the position-over-term order (RGB/POT) form of the l\_a-quasi-cyclic row- and the l\_b-quasi-cyclic column-code, respectively. This generalizes the work of Burton and Weldon for the generator polynomial of a cyclic product code (where l\_a =l\_b=1). Furthermore, we derive the generator matrix in Pre-RGB/POT form of an l\_a l\_b-quasi-cyclic product code for two special cases: (i) for l\_a=2 and l\_b=1, and (ii) if the row-code is a 1-level l\_a-quasi-cyclic code (for arbitrary l\_a) and l\_b=1.For arbitrary l\_a and l\_b, the Pre-RGB/POT form of the generator matrix of an l\_a l\_b-quasi-cyclic product code is conjectured.The spectral analysis is applied to the generator matrix of the product of an l-quasi-cyclic and a cyclic code, and we propose a new lower bound on the minimum Hamming distance of a given l-quasi-cyclic code. In addition, we develop an efficient syndrome-based decoding algorithm for l-phased burst errors with guaranteed decoding radius.