Let (X, ρ, µ) be a space of homogeneous type in the sense of Coifman and Weiss, and Y(X) a ball quasi-Banach function space on X, which supports a Fefferman-Stein vectorvalued maximal inequality, and the boundedness of the powered Hardy-Littlewood maximal operator on its associate space. The authors first introduce the Hardy space H * Y (X), associated with Y(X), via the grand maximal function, and then establish its various real-variable characterizations, respectively, in terms of radial or non-tangential maximal functions, atoms or finite atoms, and molecules. As an application, the authors give the dual space of H * Y (X), which proves to be a ball Campanato-type function space associated with Y(X). All these results have a wide range of generality and, particularly, even when they are applied to variable Hardy spaces, the obtained results are also new. The major novelties of this article exist in that, to escape the reverse doubling condition of µ and the triangle inequality of ρ, the authors cleverly construct admissible sequences of balls, and fully use the geometrical properties of X expressed by dyadic reference points or dyadic cubes and, to overcome the difficulty caused by the lack of the good dense subset of H * Y (X), the authors further prove that Y(X) can be embedded into the weighted Lebesgue space with certain special weight, and then can fully use the known results of the weighted Lebesgue space.