p value of hypothesis testing for data distribution fit python scipy

Question

I am trying to obtain the data distribution using hypothesis testing. The given data is said to follow Gamma distribution (obtained using JS divergence and MLE tests). However,