Bar graph and box-and-whisker plot sorted over binary variable


I have a dataset that looks like that below.

I want to create graphs like the example below:

The variable names var1-var4 in my dataset would replace the variable names on the x axis in the graph. My variable "status" is what is labeled Unstable/Stable in the graph. Pairid is a variable that describes which pairid each row belongs to, as this is a matched case-control study with two controls per case.

Have you got tips on how to create such a graph?
Thank you in advance.