Automated interpretable discovery of heterogeneous treatment effectiveness: A COVID-19 case study