What tool can compute SNPs from a set of assemblies and output the SNPs into a dataframe?

Question

[O.rka]

I have several assemblies all from the same bacterial strain. I want to do a SNPs analysis but I'm having trouble finding tools that will do this and output the data into a dataframe of SNPs by position.

I have the fasta files for each assembly along with the shotgun genomic reads that went into (some of) them.