I'm working on an interespecific phylogeny of psbA-trnH haplotypes and I would like to run Bayesian and Maximum Likelihood analyses considering gaps as fifth base. There are some indels blocks specific to each species, then I think this information is important in the separation of the three species involved in the analysis. Could anyone tell me which software can I use to run a Maximum Likelihood analysis considering gaps as fifth base?