[suggestion] what would be the appropriate pipeline to perform joint embedding of GEX and ATAC?

I have processed my snRNA of multiome using scanpy and snATAC using snapatac2. Both have been annotated as well.
I would like to perform joint embedding of rna and atac but the var attributes are different and so are the number of cells in both the modalities.
The overarching goal after embedding is obtain cell_data that will be used as an input for pycistopic->scenic+.