Covariance-based method for eigenstate factorization and generalized singlets