I use 4 concatenated molecular markers for phylogenetic inference. One of them, ribosomal 16s, shows a large amount of indels in my group of study, and I would like to use that information which is otherwise lost.
Could you point me to a thorough description of the best procedure for coding gaps as an additional partition and incorporating them in RAXml/MrBayes (substitution model, etc.)?