I have two spatial objects (coordinates in .shp files) that I want to draw together in R, it is using the same CRS (i.e. to map these using an identical, appropriate, projection)

(The R code to define a projection is similar to: SpatialPoints(coords, proj4string=CRS('+init=epsg:4326')))

First shape file has an associated .proj that contains: GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]

Second shape file has an associated .proj that contains: PROJCS["WGS_1984_Web_Mercator_Auxiliary_Sphere",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator_Auxiliary_Sphere"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",0.0],PARAMETER["Standard_Parallel_1",0.0],PARAMETER["Auxiliary_Sphere_Type",0.0],UNIT["Meter",1.0]]

I tried to do a search using key terms in this .proj such as https://spatialreference.org/ref/?search=WGS_1984 but don't know which of the several projections to pick.

Any help to find the proper argument to CRS() would be very useful.