ó
4ukOc           @   s;  d  d l  Z  d  d l Z  d  d l Z d  d l Z d  d l Z d  d l Z d  d l Z e j d Z e j d Z	 e j d Z
 e j d Z e j d Z e j d Z e j d Z e j d	 Z e j d
 Z e j d Z e j d Z e j d Z e j d Z e j d Z e j d Z e j d Z e  j ƒ  Z d e e  j ƒ  ƒ Z e  j e ƒ e j e e  j j  e d ƒ ƒ e  j! e ƒ e" d d ƒ Z# e# j$ d ƒ e# j$ d ƒ e d k  sÀe d k rÉd Z n  e
 d k rõe# j$ d ƒ e j d d ƒ n  e
 dF k r#e j% e
 e e d e# d d  ƒ n. e
 dG k rQe j& e
 e e d e# d d  ƒ n  e j' d$ d g ƒ Z( e" e d% ƒ Z) e j* ƒ  Z+ e# d& e e+ ƒ IJe  j, d	 ƒ e j ƒ  Z- e j. d' g d( e j/ ƒZ0 e
 d k rÞd) Z n  e d) k re j1 e0 j2 e# e ƒ n  e
 dH k r`e0 j2 j$ d* ƒ e0 j2 j$ d+ ƒ e0 j2 j$ e d ƒ e0 j2 j$ e e+ ƒ d ƒ n] e
 dI k r½e0 j2 j$ d* ƒ e0 j2 j$ d, ƒ e0 j2 j$ e d ƒ e0 j2 j$ e e+ ƒ d ƒ n  e	 d- k rÙe# j$ d. ƒ nx e	 d/ k re# j$ d0 ƒ e0 j2 j$ d1 ƒ e0 j2 j$ d2 ƒ n< e	 d3 k rQe# j$ d4 ƒ e0 j2 j$ d1 ƒ e0 j2 j$ d5 ƒ n  e j3 e e4 e0 j2 e# ƒ Z e d6 k rße5 e ƒ d k  r“d Z n  e0 j2 j$ e d ƒ e0 j2 j$ e e ƒ d ƒ e# j$ d7 e e ƒ d ƒ n  e d6 k rþe0 j2 j$ d2 ƒ n  e d8 k re0 j2 j$ d9 ƒ n  e
 dJ k r<e0 j2 j$ d: ƒ n  e0 j2 j$ d5 ƒ e0 j2 j6 ƒ  e0 j7 ƒ  e) j6 ƒ  e# j6 ƒ  e8 e- Z8 e  j, d; ƒ e" d< d ƒ Z9 e9 j$ d= e  j: ƒ  d d> e e8 ƒ d? ƒ e9 j6 ƒ  e
 d k se
 d k se
 d! k se
 d" k se
 d# k rÐe" d< d ƒ Z9 e9 j$ d ƒ e9 j$ d@ ƒ e9 j6 ƒ  e j' dA dB dC d e e ƒ d dD e  j j  e e ƒ g ƒ e j; e  j j  e e ƒ e< e< ƒ e= d8 k re j> dB e  j j  e e? ƒ ƒ qn> e j> dB e  j j  e e ƒ ƒ e j; e  j j  e e ƒ e@ ƒ eA j6 ƒ  e  j! e ƒ e jB e ƒ dE GHd S(K   iÿÿÿÿNi   i   i   i   i   i   i   i   i	   i
   i   i   i   i   i   i   s	   DISCPARS.s   infile.tempt   MSGFILEt   asM   ---------------------  DNAPARS - DNA PARSIMONY METHODS ---------------------
s   
id   t   ns    
t   infilet   bt   dt   1t   st   it   pst   pot   pwt   neditt   ws   DNAPARS:  SEED= t   dnaparst   stdint   Js   m
s   w
s   d
t   ms   Search option: THOROUGH
t   ls   Search option: LESS THOROUGH
s   s
s   n
t   rs/   Search option: REARRANGE ON JUST ONE BEST TREE
s   y
t   ts&   Using threshold parsimony, THRESHOLD= t   ys   1
s   3
i    t   outfiles   Elapsed time on s   : s    secondss,   -------------------------------------------
s   consense.pyt   outtreet   es   outfile.consenses   DNAPARS Completed(   R   R   (   s   psR
   R   (   R   R   (   s   psR
   R   (   R   R   s   psR
   R   (C   t   ost   os.patht   phylipt   shutilt
   subprocesst   syst   timet   argvt   INFILEt   SOPTIONt   METHODt
   REPLICATESt   PERCENTt	   BLOCKSIZEt   JUMBLEt   NUMJUMt   OUTGROUPt   PARSTYPEt	   THRESHOLDt   TRANSVERSIONt   TERMOUTt	   PRINTDATAt   OUTFILEt   TREEFILEt   getcwdt   STARTDIRt   strt   getpidt   TEMPDIRt   mkdirt   copyfilet   patht   joint   chdirt   opent	   msgfile_ht   writet   stdresamplet   weightless_resamplet   callt   p_testinfilet	   termout_ht   phylip_randomt   RSEEDt   nicet
   start_timet   Popent   PIPEt	   p_PROGRAMt   jumbleR   t   do_outgroupt   NUMSEQt   intt   closet   waitt   time_elapsedt	   outfile_ht   unamet	   merge_msgt   Truet   ALLTREESt   movet   ALLTREEFILEt   Falset
   OUTfile2_ht   rmtree(    (    (    s
   dnapars.pyt   <module>   sØ   	""			



-
<
:"%
