The phase diagram of the Bose-Hubbard model in the presence of off-diagonal disorder is determined using Quantum Monte Carlo simulations. A sequence of quantum glass phases intervene at the interface between the Mott insulating and the Superfluid phases of the clean system. In addition to the standard Bose glass phase, the coexistence of gapless and gapped regions close to the Mott insulating phase leads to a novel Mott glass regime which is incompressible yet gapless. Numerical evidence for the properties of these phases is given in terms of global (compressibility, superfluid stiffness) and local (compressibility, momentum distribution) observables.PACS numbers: 03.75. Lm, 03.75.Ss, 05.30.Jp, 32.80.Pj The competition between disorder, interactions and commensurability in quantum many-body systems is known to produce novel quantum glassy phases, characterized by a gapless spectrum and by the absence of a global order parameter [1].4 He adsorbed on porous media[2], granular superconductors[3], disordered magnets [4] are but a few manifestations of localization effects due to random potentials in interacting bosons. In recent years, ultracold atomic gases in magneto-optical traps have opened a new frontier in the study of strongly correlated systems, as unprecedented control over experimental parameters in these systems makes them ideally suited for studying many-body phenomena. Disorder can be generated in optical lattices by exposure to speckle lasers [5,6], incommensurate lattice-forming lasers [7,8,9], and by other means [10]. The interplay between disorder and interactions in trapped Bose-Einstein condensates has recently been explored experimentally in 87 Rb, both in the continuum [6,11,12] and in an optical lattice [9,13].Theoretically, the effects of random potentials on interacting bosons in periodic lattices have been studied using analytic [1,14,15,16] and numerical techniques [17]. It is now well established that even infinitesimally small potential disorder can destroy the direct superfluid (SF) to Mott insulator (MI) transition in one dimension by introducing an intervening insulating, but compressible, Bose glass (BG) phase [1,18]. Surprisingly, while the effects of potential disorder have been widely investigated, other kinds of disorder, e.g. hopping or interaction strengths, have remained largely unexplored until recently, when it was demonstrated that these can lead to qualitatively quite different phenomena [19,20,21]. In the quantum rotor model with off-diagonal disorder, the SF to MI transition takes place via an intermediate Mott glass (MG) phase [20,21]-a unique incompressible, yet gapless, glassy regime that was first reported in disordered fermions with extended range interactions [22]. It is thus of great interest to investigate whether such a phase also appears in the Bose-Hubbard model (BHM), and if particle-hole symmetry is essential for the stabilization of such a MG.We use large-scale QMC simulations to study the effects of off-diagonal disorder in the Bose-Hubbard model on a on...