Quantum stabilizer codes often face the challenge of syndrome errors due to error-prone measurements and multiple rounds of syndrome extraction are typically employed. In this paper, we consider phenomenological decoding problems, where data qubit errors may occur between two syndrome extractions, and each syndrome measurement can be faulty. To handle these diverse error sources, we define a generalized check matrix over mixed quaternary and binary alphabets to characterize their error syndromes. This generalized check matrix leads to the creation of a Tanner graph comprising quaternary and binary variable nodes, which facilitates the development of belief propagation (BP) decoding algorithms to tackle phenomenological errors. Additionally, our BP decoders are applicable to general sparse quantum codes.