Trustworthy Multimodal Regression with Mixture of Normal-inverse Gamma Distributions