Average var gene expression by obs metadata

Hi all,

I have a single cell expression dataset with clinical metadata. Outside of checking .obs cluster names for my gene of interest, I want to check the other metadata (tumor status, tissue site, etc.) for the expression levels of a gene of interest. I can use Vln plots, but they’re hard to read. So I wanted to get the mean expression of a gene between different .obs. Is there a quick way to calculate mean gene expression between .obs and graph it in a barplot?

Thanks!