:align intron1.fasta
:align intron2.fasta
:align intron3.fasta
:align coding.fasta

:smodel 1,2,3:GTR +> ASRV.Free(n=4)
:imodel 1,2,3:RS07

:smodel 4:|w:GTR +> x3 +> MutSel +> dNdS(w)| +> M3
:imodel 4:none

