ó ›¢yZc@sóddlZddlZddlZddlmZmZmZmZddlm Z ddlm Z m Z ddl m Z mZddlmZeejƒZdd d gZd „Zd „Zejj ejd ƒZed„Zd„ZdS(iÿÿÿÿN(t reportingtqconfigtqutilstca_utils(tcompile_aligneri(tget_path_to_programt is_python2(tjointabspath(t get_loggertUtilst SizeFastat GetFastaStatscCs>x7tD]/}tjjtjj||dƒƒstSqWtS(Ns.class(trequired_java_fnamestostpathtisfileRtFalsetTrue(tdirpatht required_name((s>/home/birch/BIRCH/local/install/quast-4.6.3/quast_libs/gage.pytall_required_java_classes_exists %c CsŸtj|ƒ}tj|ƒ}tjdtj|ƒ|dƒtjj|d|dƒ}tjj|d|dƒ} tjdtj|ƒdtjj |ƒdtjj | ƒdƒt |dƒ} t | dƒ} tj d |t t jjƒ|||ttjƒgd | d | d dtj|ƒd tƒ} | dkritjdtj|ƒdƒntjdtj|ƒdƒ| jƒ| jƒ| S(Ns s...tgage_s.stdouts.stderrsLogging to files s and twtshtstdouttstderrtindentt only_if_debugisFailed.sDone.(Rtlabel_from_fpathtlabel_from_fpath_for_fnametloggertinfot index_to_strRRRtbasenametopentcall_subprocessRRtmisctcontig_aligner_dirpathtstrRt min_contigRtclose( tit contigs_fpathtgage_results_dirpathtgage_tool_patht referencettmp_dirtassembly_labeltcorr_assembly_labelt log_out_fpatht log_err_fpatht log_out_ft log_err_ft return_code((s>/home/birch/BIRCH/local/install/quast-4.6.3/quast_libs/gage.pytrun_gages,"@  !  tgagec CsÖ|rVxItD]A}tjjt|dƒ}tjj|ƒr tj|ƒq q WtStdƒ}|dkr‡t j dtdƒt Stj ƒ}tjtƒt jdtjjtdƒdƒy{gtD]j}tjdtjjt|dƒgd ttjjtdƒd ƒd ttjjtd ƒd ƒƒ^qÎ}Wn4tk rut jd tdƒtj|ƒt SXtj|ƒtd„|Dƒƒs¦ttƒ rÒt j dtdtjsÆdndƒt StS(Ns.classtjavacsZJava compiler not found (javac)! Please install it or compile GAGE java classes manually (s /*.java)!s'Compiling JAVA classes (details are in smake.logs and make.err)s.javaRRRsmake.errsPermission denied accessing s. Did you forget sudo?css|]}|dkVqdS(iN((t.0R6((s>/home/birch/BIRCH/local/install/quast-4.6.3/quast_libs/gage.pys _ss3Error occurred during compilation of java classes (s&/*.java)! Try to compile it manually. sDYou can restart Quast with the --debug flag to see the command line.t(R RRRt gage_dirpathRtremoveRRtNoneRterrorRtgetcwdtchdirt main_infoRR$R#tIOErrortwarningtanyRRtdebug(t only_cleanRtfpatht javac_pathtcur_dirt java_fnamet return_codes((s>/home/birch/BIRCH/local/install/quast-4.6.3/quast_libs/gage.pyt compile_gageAs8     $x   # csbtjj|dƒ‰tjjˆƒs7tjˆƒntjjtdƒ‰tjƒtjdƒdddddd d d d d dddddddddddddddg}t j j t j j t j j t j jt j jt j jt j jt j jt j jt j jt j jt j jt j jt j jt j jt j jt j jt j jt j jt j jt j jt j j t j j!t j j"t j j#g}tjjˆdƒ‰tjj$ˆƒsÒtj%ˆƒnt&tƒ söt'tƒ rt(ƒ rtj)dƒdSt*t+|ƒt,j-ƒ}t.ƒrAddl/m0}m1‰nddl2m0}m1‰t,j3sš|d|ƒ‡‡‡‡‡fd†t4|ƒDƒƒ}n:gt4|ƒD]'\}} t5|| ˆˆˆˆƒ^q§}d|krtj)d tjjˆd!ƒƒdSx$t4|ƒD]\}} t6j7| ƒ} t j8| ƒ} tjjˆd"| d#ƒ} t9| d$ƒ} d}x­| D]¥}|||krt||j:dƒrÏ| j;|||j<||d%ƒd&j=ƒƒn'| j;|||j<d%ƒd&j=ƒƒ|d&7}|t+|ƒkrPqqtqtW| j>ƒqWt j?|ƒt,j@sQtAjBˆƒntjd'ƒdS((NR8sgetCorrectnessStats.shsRunning GAGE...s Total unitstMintMaxtN50s Genome Sizes Assembly Sizes Chaff basessMissing Reference BasessMissing Assembly BasessMissing Assembly ContigssDuplicated Reference BasessCompressed Reference BasessBad TrimsAvg IdytSNPss Indels < 5bps Indels >= 5t Inversionst Relocationt Translocationt BasesInFastattmpsTGAGE module was not installed properly, so it is disabled and you cannot use --gage.iÿÿÿÿ(tParalleltdelayedtn_jobsc3s6|],\}}ˆtƒ||ˆˆˆˆƒVqdS(N(R7(R:R*R+(RXR,R-t ref_fpatht tmp_dirpath(s>/home/birch/BIRCH/local/install/quast-4.6.3/quast_libs/gage.pys ’sisXError occurred while GAGE was processing assemblies. See GAGE error logs for details: %ss gage_*.stderrRs.stdouttrt:isDone.(CRRRtisdirtmkdirR<Rtprint_timestampRBRtFieldstGAGE_NUMCONTIGStGAGE_MINCONTIGtGAGE_MAXCONTIGtGAGE_N50tGAGE_GENOMESIZEtGAGE_ASSEMBLY_SIZEtGAGE_CHAFFBASEStGAGE_MISSINGREFBASEStGAGE_MISSINGASMBLYBASEStGAGE_MISSINGASMBLYCONTIGStGAGE_DUPREFBASEStGAGE_COMPRESSEDREFBASESt GAGE_BADTRIMt GAGE_AVGIDYt GAGE_SNPStGAGE_SHORTINDELStGAGE_LONGINDELStGAGE_INVERSIONStGAGE_RELOCATIONtGAGE_TRANSLOCATIONtGAGE_NUMCORCONTIGStGAGE_CORASMBLYSIZEtGAGE_MINCORCONTIGtGAGE_MAXCORCOTINGt GAGE_CORN50texiststmakedirsRRRMR?tmintlenRt max_threadsRtjoblibRWRXtjoblib3tmemory_efficientt enumerateR7RRtgetR#t startswitht add_fieldtsplittstripR)t save_gageRFtshutiltrmtree(RZtcontigs_fpathstoutput_dirpathtmetricstmetrics_in_reportingRYRWRLR*R+R1treportR2t logfile_outt cur_metric_idtline((RXR,R-RZR[s>/home/birch/BIRCH/local/install/quast-4.6.3/quast_libs/gage.pytdogsv    $   $7    2'    (tloggingRRŠt quast_libsRRRRtquast_libs.ca_utils.miscRRRtos.pathRRtquast_libs.logR tLOGGER_DEFAULT_NAMERR RR7Rt LIBS_LOCATIONR<RRMR”(((s>/home/birch/BIRCH/local/install/quast-4.6.3/quast_libs/gage.pyts   "  ! &