We re-examine electron and hole transport in pure and stabilized amorphous selenium (a-Se) and attempt to construct a DOS distribution in the mobility gap below E c and above E v based on time-of-flight (TOF) transient photoconductivity measurements. First, we review the current status of a-Se, its recent use in commercial X-ray detectors, and the scientific information that is available in terms of its density of states (DOS). We review and describe a convenient multiple trapping transport mechanism to calculate the TOF transient photocurrents in a system with a distributed DOS in the mobility gap and invoke a number of assumptions that allow the photocurrent to be calculated as an inverse Laplace transform. The assumptions behind the model are critically examined. Then, by comparing the calculated TOF photocurrent shapes directly with experimental waveforms, the DOS distribution has been extracted in a-Se in the vicinity of conduction and valence bands (below E c and above E v ). Below the conduction band, the DOS distribution seems to demonstrate three peaks at around E c -0.30 eV, E c -0.45 eV and below E c -0.55 eV. First principles atomistic modeling provided insight into the origins of these peaks and the peaks have been assigned to certain valence-alteration pairs (VAPs) appearing in different configurations. In contrast to the conduction band vicinity, the DOS distribution above the valence band seems to be almost totally featureless. Within the range of experimental uncertainties, a featureless DOS below E v down to 0.55 eV seems to be able to explain a vast range of data, that is, it can explain hole transport not only as a function of the electric field and temperature down to 123 K, but also the thickness dependence of the mobility over two orders of magnitude. At room temperature, the hole transport is nondispersive but at temperatures below *200 K, it becomes dispersive exactly along the lines argued by Pfister and Scher (Adv Phys 27:747, 1978), and the mobility evinces a clear thickness dependence. The first principles calculations yield two defect levels related to VAP type defects in this region of the mobility gap but the valence band tail is so broad that it masks these peaks. This seeming asymmetry between the conduction and valence bands may be related to their physical origins: the conduction band arises from antibonding states and valence band is due to interacting lone pairs.