ó ¼Tc@s¦ddlZddlZddlZddlZddlZddlZddlZddl Z ddl Z d„Z d„Z dd„Zddd„Zd„ZdS( iÿÿÿÿNc Csìd}t|dƒ.}x$|D]}|jƒ}|d7}qWWdQXtj||fdtjƒ}d}t|dƒn}xd|D]\}|jƒjdƒ}x4tt|ƒƒD] }t||ƒ|||f;ss.{0}/armatus -i {1} -g {2} -o {3}/{4} -k {5} -mtgammaiÿÿÿÿs{0}/{1}t.s.gamma.iis rm -rf {0}(tostpathtexiststmakedirstformattsystemtmyutilt listfilestfindRRtjoinRt setdefaultRt itertoolst combinationstlisttAssertionErrortkeys(R tmaxscalet outprefixt armaoutfoldert ARMATUSPATHttopcountt intersecttcodet scale2domst tarmaoutfiletoutpathstoutpathRtscaletdom1tdom2((s../sedfm/lib/ArmatusLib.pyt runArmatus/s&  C  6%)gš™™™™™©?c CsÐd„}tjj|ƒs+tj|ƒndj|||||||ƒ}tj|ƒdj||ƒ} t| ƒ} x<tjt | ƒdƒD]"\} } || | ƒ st ‚qWtjdj|ƒƒ| S(sáRuns Armatus only for consensus Args: filepath: maxscale: outprefix: armaoutfolder: ARMATUSPATH: topcount: stepsize: Returns: consdomains: consensus domains cSs8|\}}|\}}||ks0||kr4tStS(N(RR(RRR R!R"R#((s../sedfm/lib/ArmatusLib.pyR$]ss5{0}/armatus -i {1} -g {2} -o {3}/{4} -m -k {5} -s {6}s{0}/{1}.consensus.txtis rm -rf {0}( R'R(R)R*R+R,RR2R3R4R5( R R7R8R9R:R;tstepsizeR<R=RAt consdomainsRCRD((s../sedfm/lib/ArmatusLib.pytrunArmatusConsensusPs !  %cCsZd„}tjj|ƒs+tj|ƒndj||||||ƒ}tj|ƒgtj|ƒD]-}|jdƒdkrfdj||ƒ^qf} i} xž| D]–} t | ƒ} t dj | j |dƒdj dƒdd !ƒƒ} t | j |dƒdj dƒd ƒ}| j| iƒt| ƒ| | |s            !