I need to compute pairwise sequence divergence between two species. I have the fasta files of the two sequences aligned to an outgroup. I was wondering if you can suggest a program for doing this. It is important if the method is able to handle missing data - aka ignore missing data. I can comfortable with coding in perl and can implement something myself but I think parsing two fasta files will be quite slow. Any suggestions?