Data faults in sensor networks must be marked to ensure accurate inferences. We introduce a two phase semi-realtime end-to-end Bayesian fault detection system for sensor networks. The first phase selects a subset of agreeing sensors from which a model of expected behavior is derived. The second phase uses this subset to derive and tag questionable sensor data. To accurately model the data, we use a hierarchical Bayesian space-time (HBST) model, as compared to the linear autoregressive modeling used in previous works. Applying this system to simulated and real world data, results are excellent when the phenomenon is well modeled by the HBST model. It achieves high detection rates and almost zero false detection rates. Results also indicate that in cases of critically low spatial sampling density a more accurate model is required.