Display spectrum of region of interest with Rasterio and Geopandas

Question

I believe this is a very simple question, but for some reason I really struggle with finding the way to do that.

I have plot a satellite image that has 13 bands and