Make a custom track in the UCSC Genome Browser, setting bigWig-sourced signal tracks in "dense" mode, and arrange or sort tracks in your desired ordering. Signal for each row is taken from your bigWig files.
The "dense" mode will give you a figure very heatmap-like in appearance, e.g.:
You can modify the custom track to add color to rows on a per-signal basis. You might assign colors to groups or categories of sample types/conditions/treatments/etc.
If there are annotation tracks you want to add, you can put them into the browser.
Once you have your browser and custom track file arranged the way you want, you can use a tool like soda.py to generate a gallery of browser snapshots, taken at each position in an input BED file.
This can be a single region, per your question, or multiple regions. For example, to make the figure above, I simply echo'ed an ad-hoc BED interval to a text file called
adHoc.bed, and used that as my input:
$ echo -e "chr9\t133252000\t133280861" > adHoc.bed
If you have multiple regions of interest, this input BED file does not need to be in
sort-bed order. Often it is useful to sort regions by other criteria, such as ID value in the fourth column, or score or signal value in the fifth column. This determines the order of snapshots in the gallery.