ó £žbc@€sãddlmZddlZddlZddlmZddlmZmZm Z m Z m Z ddl m Z ddlmZe e jƒZiZddd„ƒYZd „Zd „Zd „Zd „Zd dd„ƒYZdS(iÿÿÿÿ(twith_statementN(t defaultdict(t fastaparsert genes_parsert reportingtqconfigtqutils(t get_logger(t run_paralleltFeatureContainercB€seZdd„ZRS(tcC€s(||_||_g|_i|_dS(N(tkindtfpathst region_listtchr_names_dict(tselfR R ((sH/home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/genome_analyzer.pyt__init__s   (t__name__t __module__R(((sH/home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/genome_analyzer.pyR scC€stS(N(tref_lengths_by_contigs(((sH/home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/genome_analyzer.pytget_ref_aligned_lengthsscC€sDi}x=|D]5}|j|kr5|j||j8ss{Reference names in file with genomic features of type "%s" do not match any chromosome. Check your genomic feature file(s).s‰Some of the reference names in file with genomic features of type "%s" does not match any chromosome. Check your genomic feature file(s).N( tseqnameRtlentpoptloggertnoticetalltvaluestwarning(tfeaturetregionst chr_namestregion_2_chr_nametregionR((sH/home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/genome_analyzer.pytchromosomes_names_dict"s( ;      c6C€sä tj|ƒ}tj|ƒ}tƒ} ttƒ} tjdtj|ƒ|ƒt j j ||dƒ} t j r|| } n | d} t j j| ƒs·tjd| dddƒd2Si} x/|jƒD]!\}}dg|d| |`streverses¥Analysis of genes and/or operons files (provided with -g and -O) requires extensive RAM usage, consider running QUAST without them if memory consumption is critical.t|i i s?Something went wrong and chromosome names in your coords file (sS) differ from the names in the reference. Try to remove the file and restart QUAST.Rtstarttendtcontigtstart_in_contigt end_in_contigs _gaps.txts /dev/nulltws t t gaps_countt_fullt_partialt_genomic_features_s.txts%s %s %s %s %s sID or #tStarttEndtTypetContigt=i2s# is%s %d %d complete %s toperont,cS€s|jS(N(R,(tblock((sH/home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/genome_analyzer.pyR)Üss%s %d %d partial %s sAnalysis is finished.(NN(sID or #R7R8sTypesContig(;Rtlabel_from_fpathtlabel_from_fpath_for_fnametdictRtintRtinfot index_to_strtostpathtjoinRtuse_all_alignmentstisfileterrorRtitemsRt read_fastatsortedt enumeratetTruetappendRtmemory_efficientR topentsplittstript AlignedBlocktrangetkeystsumtspace_efficienttendswithtremovet analyze_gapstwritet min_gap_sizetstrRtFieldstGENEStOPERONSR R tlowertidtnumbertFalseRR-R,tformat_gene_infotmintmaxtmin_gene_overlaptclose(6t contigs_fpathtindextcoords_dirpathtgenome_stats_dirpathtreference_chromosomestns_by_chromosomest containerstassembly_labeltcorr_assembly_labeltresultst ref_lengthstcoords_base_fpatht coords_fpathtgenome_mappingRtchr_lent contig_tuplestsorted_contig_tuplestsorted_contigs_namest contigs_ordertidxtnamet_tfeatures_in_contigstoperons_in_contigstaligned_blocks_by_contig_nametgene_searching_enabledt coordfiletlinets1te1ts2te2t contig_nametiR3t gaps_fpatht gaps_filet cur_gap_sizet containert total_fullt total_partialt found_fpatht found_filet found_listR%t gene_blockst contig_idtcur_feature_is_foundt cur_blockt contig_infoR>tunsorted_features_in_contigstunsorted_operons_in_contigs((sH/home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/genome_analyzer.pytprocess_single_fileBs     !    #### " ($ 0  =   * $$   &1  :-c; C€s tjj|tjƒ}ddlm}|jrItjj|dƒ}ntj ƒtj dƒtjj |ƒs‚tj |ƒnt j|ƒ\} } } tjj|dƒ} t| dƒ} g}x3|jƒD]%\}}|jt|g|ƒƒqÑW|stjddd ƒn|r5|jt|d ƒƒntjd dd ƒx*|D]"}|jsdqOnx/|jD]$}|jtj||jƒ7_qnWt|jƒd krìtjd |jddd ƒ| jd|jdddƒqOtjdtt|jƒƒd|jdƒ| jd|jdtt|jƒƒdƒt|j|jt| j ƒƒƒ|_!qOWdL\}}x¬|D]¤}t#j$|ƒ}d }xa|D]Y}|jd krðt|jƒ}|j%t#j&j't|jƒƒqª|t|jƒ7}qªW|rˆ|}|j%t#j&j(|ƒqˆqˆWi}i}i}i}g}g}g}tj)}t*t|ƒtj+ƒ} gt,|ƒD]'\}!}||!||| | |f^qˆ}"t-t.|"| dt/ƒ\}#}$|t|ƒt|#ƒ7}|t_)|#stj dƒ| j0ƒdSx?| D]7}%gt1t|#ƒƒD]}&|#|&|%^q7t2|%ƒ|0fƒ|jt?|j<t#j&j=ƒƒƒx˜t#j&j:|1|2ft#j&j;|3|4fgD]l\}5}6}7|6dkrD|7dkrD| jd6dOƒq | jd6|6|7fƒ|j%|5d8|6|7fƒq W| jdƒqœW| j0ƒtj@rôdd9lAmB}8|rÏ|8jC||d:||ƒn|rô|8jC||d*||ƒqôntjDr d;d<lEmF}9dd=lGmH}:|r€|9jI||||d>d?ƒ|9jJ||||:||d@d?ƒ|9jK|||dAdBƒn|rã|9jI||||dCd*ƒ|9jJ||||:||dDd*ƒ|9jK|||dEdFƒn|9jK|||dGdHdIdJƒntj dKƒ|S(PNiÿÿÿÿ(tsearch_references_metatrawsRunning Genome analyzer...sgenome_info.txtR1sbNo file with genomic features were provided. Use the --features option if you want to specify it. Rs R<sPNo file with operons were provided. Use the -O option if you want to specify it.isNo genomic features of type "s" were loaded.sGenomic features of type "s " loaded: Rs s Loaded s genomic features of type "t"tfilter_resultss.Genome analyzer failed for all the assemblies.sreference chromosomes: s s (total length: s bp, stotal length without N's: s bp, maximal covered length: s bp) stotal genome size: s sgap min size: spartial gene/operon min size: s8%-25s| %-10s| %-12s| %-10s| %-10s| %-10s| %-10s| %-10s| tassemblytgenomet duplicationtgapstgenestpartialtoperonsR tfractiontratioReR;ixR3R4R5s%-25s| %-10s| %-12s| %-10s|is %-10s| %-10s|t-s %s + %s part(t html_savertfeaturesi(tplotter(tcontigs_aligned_lengthss/features_cumulative_plotsgenomic featuress/features_frcurve_plots/complete_features_histograms# complete genomic featuress/operons_cumulative_plots/operons_frcurve_plots/complete_operons_histograms# complete operonss/genome_fraction_histogramsGenome fraction, %t top_valueidsDone.(NN(sassemblysgenomeR¥R¦sgenesspartialsoperonsspartial(R sfractionsratiosnumberR sgenesR soperons(R¬R¬(LRERFRGRtaligner_output_dirnamet quast_libsRŸtis_quast_first_runRtprint_timestampt main_infotisdirtmkdirRtget_genome_statsRRRKRPR RR R Rtget_genes_from_fileR RR R]RCR_R&tlistRWRRRtgett add_fieldR`t REF_OPERONSt REF_GENESt_num_nf_errorsRht max_threadsRNRRžRORkRVRRiR^RjtzipRtname_from_fpathRXRaRbt get_fieldt MAPPEDGENOMEtDUPLICATION_RATIOtfloatt html_reporttquast_libs.html_saverR­tsave_features_in_contigst draw_plotsR R¯tquast_libs.ca_utils.miscR°tgenes_operons_plottfrc_plott histogram(;t ref_fpathtaligned_contigs_fpathstoutput_dirpatht features_dicttoperons_fpathst detailed_contigs_reports_dirpathRoRnRŸt genome_sizeRpRqt result_fpathtres_fileRrR!t feature_fpathR‘tfpatht ref_genes_numtref_operons_numRltreporttgenomic_featurestfiles_features_in_contigst"files_unsorted_features_in_contigstfiles_operons_in_contigst!files_unsorted_operons_in_contigst genome_mappedtfull_found_genestfull_found_operonst num_nf_errorstn_jobsRmtparallel_run_argsRvtresults_genes_operons_tuplestrefRRRzt aligned_lenRuRœR‚RRƒt assembly_nameR3t genes_fullt genes_partt operons_fullt operons_parttfieldtfulltpartR­R¯R°((sH/home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/genome_analyzer.pytdoñs         "#//+   " 7    5 O    (      6"%      RUcB€s,eZddddddd„Zd„ZRS(cC€s:||_||_||_||_||_||_dS(N(RR,R-R.R/R0(RRR,R-R.R/R0((sH/home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/genome_analyzer.pyR§s      cC€sÒ|j|j}}|j|jkr[|j|j}||krN||7}q[||8}n|j|jkr¯|jt|j|jƒ}||kr¢||}q¯||}n|jdt|ƒdt|ƒS(Nt:R¬(R/R0R,R-RiR.R_(RR%R,R-t region_shiftt region_size((sH/home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/genome_analyzer.pyRg¯s      N(RRRRRg(((sH/home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/genome_analyzer.pyRU¦s(((t __future__RtloggingREt collectionsRR³RRRRRtquast_libs.logRtquast_libs.qutilsRtLOGGER_DEFAULT_NAMERRR RR&RžRôRU(((sH/home/psgendb/BIRCHDEV/install/quast-5.2.0/quast_libs/genome_analyzer.pyts  (  ¯ µ