Python Script for Calculating Joint Maximum likelihood estimate of nucleotide polymorphism (theta = (4Nu)) based on number of segregating sites