Voglio analizzare delle sequenze in cui penso ci possa essere della selezione positiva usando PAML. Ho piu' di 100 sequenze e mi e' stato consigliato di usare MrBayes per creare l'albero filogenetico.
Ho un paio di domande:
1) Ho letto che PAML considera le indels nell'allineamento nucleotidico come caratteri ambigui e che e' meglio togliere tutte le posizioni dove le indels sono presenti o lasciare solo quelle dove quasi tutte le sequenze hanno dati.
Per questo ho usato in PAML:
cleandata = 1 *remove sites with ambiguity data (1:YES, 0:no)
Per l'albero filogenetico (MrBayes) volevo convertire le indels in linguaggio binario e usarle come una "separate binary data partition" per sfruttare al massimo l'informazione filogenetica dei miei dati.
Vorrei sapere se usuare questo albero in PAML potrebbe essere un problema dato che alcune delle posizioni usate in MrBayes non saranno usate da PAML.
2)Posso usare le branch lengths ottenute in MrBayes e usare: