Let there be given a probability measure μ on the unit circle T of the complex plane and consider the inner product induced by μ . In this paper we consider the problem of orthogonalizing a sequence of monomials { z r k } k , for a certain order of the r k ∈ Z , by means of the Gram–Schmidt orthogonalization process. This leads to a sequence of orthonormal Laurent polynomials { ψ k } k . We show that the matrix representation with respect to { ψ k } k of the operator of multiplication by z is an infinite unitary or isometric matrix allowing a ‘snake-shaped’ matrix factorization. Here the ‘snake shape’ of the factorization is to be understood in terms of its graphical representation via sequences of little line segments, following an earlier work of S. Delvaux and M. Van Barel. We show that the shape of the snake is determined by the order in which the monomials { z r k } k are orthogonalized, while the ‘segments’ of the snake are canonically determined in terms of the Schur parameters for μ . Isometric Hessenberg matrices and unitary five-diagonal matrices (CMV matrices) follow as a special case of the presented formalism.