Bacterial biofilms can form in porous media that are of interest in industrial applications ranging from medical implants to biofilters as well as in environmental applications such as in situ groundwater remediation, where they can be critical locations for biogeochemical reactions. The presence of biofilms modifies porous media topology and hydrodynamics by clogging pores and consequently solutes transport and reactions kinetics. The interplay between highly heterogeneous flow fields found in porous media and microbial behavior, including biofilm growth, results in a spatially heterogeneous biofilm distribution in the porous media as well as internal heterogeneity across the thickness of the biofilm. Our study leverages highly resolved three-dimensional X-ray computed microtomography images of bacterial biofilms in a tubular reactor to numerically compute pore-scale fluid flow and solute transport by considering multiple equivalent stochastically generated internal permeability fields for the biofilm. We show that the internal heterogeneous permeability mainly impacts intermediate velocities when compared with homogeneous biofilm permeability. While the equivalent internal permeability fields of the biofilm do not impact fluid-fluid mixing, they significantly control a fast reaction. For biologically driven reactions such as nutrient or contaminant uptake by the biofilm, its internal permeability field controls the efficiency of the process. This study highlights the importance of considering the internal heterogeneity of biofilms to better predict reactivity in industrial and environmental bioclogged porous systems.