Identify amino acid mutation spread in large population

[Biophylo2001]

Is there any way to Identify amino acid mutations spread in large population with .pb fileand the -matutils -summary command? I can only get a tsv file with nucleotide mutation with their number of occurrences. So, I was wondering if we could get the same for amino acid mutation/mutation in covid-19 genes (ORF1ab, S,E etc etc )
Context: I am studying the Covid-19 mutation in a specific country with around 3500 samples to identify which mutation has spread in more population or is spreading more rapidly.