\begin{Verbatim}[commandchars=\\\{\}]
\PYG{g+gp}{\PYGZgt{}\PYGZgt{}\PYGZgt{} }\PYG{k+kn}{from} \PYG{n+nn}{Bio.Phylo.PAML} \PYG{k+kn}{import} \PYG{n}{codeml}
\PYG{g+gp}{\PYGZgt{}\PYGZgt{}\PYGZgt{} }\PYG{n}{cml} \PYG{o}{=} \PYG{n}{codeml}\PYG{o}{.}\PYG{n}{Codeml}\PYG{p}{()}
\PYG{g+gp}{\PYGZgt{}\PYGZgt{}\PYGZgt{} }\PYG{n}{cml}\PYG{o}{.}\PYG{n}{alignment} \PYG{o}{=} \PYG{l+s+s2}{\PYGZdq{}Tests/PAML/alignment.phylip\PYGZdq{}}
\PYG{g+gp}{\PYGZgt{}\PYGZgt{}\PYGZgt{} }\PYG{n}{cml}\PYG{o}{.}\PYG{n}{tree} \PYG{o}{=} \PYG{l+s+s2}{\PYGZdq{}Tests/PAML/species.tree\PYGZdq{}}
\PYG{g+gp}{\PYGZgt{}\PYGZgt{}\PYGZgt{} }\PYG{n}{cml}\PYG{o}{.}\PYG{n}{out\PYGZus{}file} \PYG{o}{=} \PYG{l+s+s2}{\PYGZdq{}results.out\PYGZdq{}}
\PYG{g+gp}{\PYGZgt{}\PYGZgt{}\PYGZgt{} }\PYG{n}{cml}\PYG{o}{.}\PYG{n}{working\PYGZus{}dir} \PYG{o}{=} \PYG{l+s+s2}{\PYGZdq{}./scratch\PYGZdq{}}
\PYG{g+gp}{\PYGZgt{}\PYGZgt{}\PYGZgt{} }\PYG{n}{cml}\PYG{o}{.}\PYG{n}{set\PYGZus{}options}\PYG{p}{(}\PYG{n}{seqtype}\PYG{o}{=}\PYG{l+m+mi}{1}\PYG{p}{,}
\PYG{g+gp}{... }        \PYG{n}{verbose}\PYG{o}{=}\PYG{l+m+mi}{0}\PYG{p}{,}
\PYG{g+gp}{... }        \PYG{n}{noisy}\PYG{o}{=}\PYG{l+m+mi}{0}\PYG{p}{,}
\PYG{g+gp}{... }        \PYG{n}{RateAncestor}\PYG{o}{=}\PYG{l+m+mi}{0}\PYG{p}{,}
\PYG{g+gp}{... }        \PYG{n}{model}\PYG{o}{=}\PYG{l+m+mi}{0}\PYG{p}{,}
\PYG{g+gp}{... }        \PYG{n}{NSsites}\PYG{o}{=}\PYG{p}{[}\PYG{l+m+mi}{0}\PYG{p}{,} \PYG{l+m+mi}{1}\PYG{p}{,} \PYG{l+m+mi}{2}\PYG{p}{],}
\PYG{g+gp}{... }        \PYG{n}{CodonFreq}\PYG{o}{=}\PYG{l+m+mi}{2}\PYG{p}{,}
\PYG{g+gp}{... }        \PYG{n}{cleandata}\PYG{o}{=}\PYG{l+m+mi}{1}\PYG{p}{,}
\PYG{g+gp}{... }        \PYG{n}{fix\PYGZus{}alpha}\PYG{o}{=}\PYG{l+m+mi}{1}\PYG{p}{,}
\PYG{g+gp}{... }        \PYG{n}{kappa}\PYG{o}{=}\PYG{l+m+mf}{4.54006}\PYG{p}{)}
\PYG{g+gp}{\PYGZgt{}\PYGZgt{}\PYGZgt{} }\PYG{n}{results} \PYG{o}{=} \PYG{n}{cml}\PYG{o}{.}\PYG{n}{run}\PYG{p}{()}
\PYG{g+gp}{\PYGZgt{}\PYGZgt{}\PYGZgt{} }\PYG{n}{ns\PYGZus{}sites} \PYG{o}{=} \PYG{n}{results}\PYG{o}{.}\PYG{n}{get}\PYG{p}{(}\PYG{l+s+s2}{\PYGZdq{}NSsites\PYGZdq{}}\PYG{p}{)}
\PYG{g+gp}{\PYGZgt{}\PYGZgt{}\PYGZgt{} }\PYG{n}{m0} \PYG{o}{=} \PYG{n}{ns\PYGZus{}sites}\PYG{o}{.}\PYG{n}{get}\PYG{p}{(}\PYG{l+m+mi}{0}\PYG{p}{)}
\PYG{g+gp}{\PYGZgt{}\PYGZgt{}\PYGZgt{} }\PYG{n}{m0\PYGZus{}params} \PYG{o}{=} \PYG{n}{m0}\PYG{o}{.}\PYG{n}{get}\PYG{p}{(}\PYG{l+s+s2}{\PYGZdq{}parameters\PYGZdq{}}\PYG{p}{)}
\PYG{g+gp}{\PYGZgt{}\PYGZgt{}\PYGZgt{} }\PYG{n+nb}{print}\PYG{p}{(}\PYG{n}{m0\PYGZus{}params}\PYG{o}{.}\PYG{n}{get}\PYG{p}{(}\PYG{l+s+s2}{\PYGZdq{}omega\PYGZdq{}}\PYG{p}{))}
\end{Verbatim}
