This python script takes a no-gap amino acid aligned file of genes. Each gene should start with a ">". It then finds all the fourfold degenerate codons shared by any two sequences and gives the number of SNP differences at the third codon of any fourfold degenerate site.