Visualisation of ATAC seq data in Python

Hello Guys,

I’ve just started analysing my first sc ATAC seq dataset.
I’m using snapATAC2 and so far I’ve been able to preprocess, batchcorrect and cluster without any major problems.

I’am however having problems creating coverage plots that look good. I’ve tried using seaborn with and mathplotlib but nothing that I create myself looks nearly as good as anything found in published papers.

Can you recommend a library in python that would have this as a built in function?
I would like to create a plot similar to this for the region of my interest:

Very thankful for any sort of advice I could get regarding this.

All the best