The generalized Chaplygin gas (GCG), with the equation of state p = -A/rho(alpha), was recently proposed as a candidate for dark energy in the universe. In this paper we confront the GCG with Type Ia supernova (SN Ia) data using available samples. Specifically, we have tested the GCG cosmology in three different classes of models with (1) Omega(m) = 0.3 and Omega(Ch) = 0.7, (2) Omega(m) = 0.05 and Omega(Ch) = 0.95, and (3) Omega(m) = 0 and Omega(Ch) = 1, as well as a model without prior assumptions on Omega(m). The best-fit models are obtained by minimizing the chi(2) function. We supplement our analysis with confidence intervals in the (A(0), alpha)-plane by marginalizing the probability density functions (pdf's) over the remaining parameters assuming uniform priors. We have also derived one-dimensional pdf's for Omega(Ch) obtained from joint marginalization over alpha and A(0). The maximum value of such a pdf provides the most probable value of Omega(Ch) within the full class of GCG models. The general conclusion is that SN Ia data give support to the Chaplygin gas (with alpha = 1). However, a noticeable preference for A(0)-values close to 1 means that the alpha dependence becomes insignificant. This is reflected in one-dimensional pdf's for alpha that turned out to be flat, meaning that the power of the present supernova data to discriminate between various GCG models (differing by alpha) is weak. Extending our analysis by relaxing the prior assumption of the flatness of the universe leads to the result that even though the best-fit values of Omega(k) are formally nonzero, they are still close to the flat case. Our results show clearly that in GCG cosmology, distant (i.e., z > 1) supernovae should be brighter than in the Lambda CDM model. Therefore, one can expect that future supernova experiments (e.g., SNAP) having access to higher redshifts will eventually resolve the issue of whether the dark energy content of the universe could be described as a Chaplygin gas. Moreover, it would be possible to differentiate between models with various values of the alpha-parameter and/or discriminate between GCG, Cardassian, and Lambda CDM models. This discriminative power of the forthcoming mission has been demonstrated on simulated SNAP data.