[_Sc@@sVddlmZmZddlZddlZddlmZddlmZddl m Z ddl m Z m Z mZddlTddlmZdd lmZmZd dlTd d lmZmZd d lmZmZmZmZmZmZm Z m!Z!m"Z"m#Z#m$Z$m%Z%m&Z&d Z'ddZ(dZ)dZ*dZ+dS(i(tabsolute_importtprint_functionN(t OptionGroup(tTemplate(tseq_io(tSeqtSeqListtnucleic_alphabet(t*(tDeOptionParser(t iteritemstStringIOi(t ColorSchemet ColorGroup( t LogoOptionstLogoDatat LogoFormatt parse_priort descriptiontrelease_descriptiont formatterstdefault_formattert std_alphabetst std_unitst std_sizeststd_color_schemest read_seq_datacC@s1t}|jtjd\}}|r?|jd|n|jret|jtjdnypt |}t ||}|j }|||}tj ddkr|j jj|n|j j|WnUtk r }td|dtjtjdn tk r,}tjdnXdS( sWebLogo command line interface isUnparsable arguments: %s iisError:tfileiN(t_build_option_parsert parse_argstsystargvterrortservethttpd_serve_forevertporttexitt_build_logodatat_build_logoformatt formattert version_infotfouttbuffertwritet ValueErrortprinttstderrtKeyboardInterrupt(tparsertoptstargstdatatformatR'tlogoterr((s9/home/psgendb/BIRCHDEV/pkg/weblogo-3.4/weblogolib/_cli.pytmainAs(     ic C@s3tjddkr4ddlj}ddlj}nddl}ddl}d|jfdY}tjdd}|dtj j tj dj d7}|tj d(R?R<((s9/home/psgendb/BIRCHDEV/pkg/weblogo-3.4/weblogolib/_cli.pyt is_python|s(t__name__t __module__R@RA(((s9/home/psgendb/BIRCHDEV/pkg/weblogo-3.4/weblogolib/_cli.pyR8rs t PYTHONPATHR9t:thtdocss.WebLogo server running at http://localhost:%d/(RR(t http.servertservertBaseHTTPServert CGIHTTPServertCGIHTTPRequestHandlertostgetenvR<tabspathtsplittenvirontresource_filenameRBt__file__tchdirt HTTPServerR-t serve_foreverR/R$( R#RHt cgiserverR8t pythonpathRFt HandlerClasst ServerClassthttpd((s9/home/psgendb/BIRCHDEV/pkg/weblogo-3.4/weblogolib/_cli.pyR"fs&  +    c C@sEt}|j}|dkr3ttjj}ny2ddlm}|j |d|j }t }WnRt k r}|j dkr|nt||j jd|j d|j}nX|rL|jrt dn|js|jr|jn|js |jr|jnt|j|j |j}tj|j ||}n|js^|jrtg|D]} | j^qh|j }n|js|jrtj|j st dn|j } t|_ tg|D]} t| |j j^q|j }| |_ nt|j|j |j}tj||}|S(Ni(tMotiftalphabetttransfactignore_lower_cases@error: option --ignore-lower-case incompatible with matrix inputs+non-nucleic sequence cannot be complemented(R:tfintNoneR Rtstdintreadtcorebio.matrixR[t read_transfacR\R>R,t input_parserRR^treversetrevcompt complementRt compositiontweightRt from_countsRRt alphabeticRt from_seqs( toptionst motif_flagR_R[tmotift motif_errtseqstpriorR3tstaaa((s9/home/psgendb/BIRCHDEV/pkg/weblogo-3.4/weblogolib/_cli.pyR%sH         .  7 c C@smi}ddddddddd d d d d dddddddddddddddg}x|D]}|j|||dd?|jd@dAddBddddCddDddEddF|jd@dGddHddddIddddJdd?|jdKdLddMddddNddOddPddQ|jdRdSddTddddNddUddQ|jdVdWddXddddNddYddQ|jd@dZdd[dd\dtdd]|jd@d^dd_dd\dtdd`|jd@daddbdd\dtddc|jd@ddddedd\dtddf|jdgdhddidddd#d$tddjd|jddk|jdldmddnddddNddod|jddp|jdqdrddsddddCddtd|jddu|jd@dvddwddddCddxd|jddu|jdydzdddd{d|jdd|dd}|jd~dddddddCd|jddddu|jd@dddddddCddddddu|jddddddddIdddd|jdddddd{ddd|jdd|dd|jddddddddCddddu|jdddddd{d|j dd|dd|jddddddddCd|j!ddddu|jd@dddddddId|j"dddd?|jd@ddddddd{d|j#dd|dd|jd@ddddddd{d|j#dd|ddtt$j} | j%|jdddddddd#d$t$ddddddd&j | |jdddddddddddgdd|jd@dddddddd|j&dd|jddddiddddId|jdd|jdd|jd@dddddddId|j'dd|j'dd|jd@ddddddd{dtdd|dd|jd@dddddddIdddddd|jd@ddddddd{dt(dd|dd|jd@ddddd{d|j)dd|dd|jd@ddddd\dtdd|jd@dddddddNdddddd|S(Ntusages7%prog [options] < sequence_data.fa > sequence_logo.epsRtversiontadd_verbose_optionssInput/Output OptionssLogo Data OptionstTransformationss.Optional transformations of the sequence data.sLogo Format Optionss9These options control the format and display of the logo.s Color OptionssFColors can be specified using CSS2 syntax. e.g. 'red', '#FF0000', etc.sAdvanced Format OptionssAThese options provide fine control over the display of the logo. sWebLogo Servers+Run a standalone webserver on a local port.s-fs--fintdestR_tactiontstorettypetfile_intdefaultthelps$Sequence input file (default: stdin)tmetavartFILENAMER]s-Ds --datatypeRetdicttchoicessRType of multiple sequence alignment or position weight matrix file: (%s, transfac)s, itFORMATs-os--foutR)tfile_outsOutput file (default: stdout)s-Fs--formatR'sIFormat of output: eps (default), png, png_print, pdf, jpeg, svg, logodatas-As--sequence-typeR\s5The type of sequence data: 'protein', 'rna' or 'dna'.tTYPEs-as --alphabetsThe set of symbols to count, e.g. 'AGTC'. All characters not in the alphabet are ignored. If neither the alphabet nor sequence-type are specified then weblogo will examine the input data and make an educated guess. See also --sequence-type, --ignore-lower-cases-Us--unitsRtchoicesA unit of entropy ('bits' (default), 'nats', 'digits'), or a unit of free energy ('kT', 'kJ/mol', 'kcal/mol'), or 'probability' for probabilitiestNUMBERR9s --compositionRitstringtautos{The expected composition of the sequences: 'auto' (default), 'equiprobable', 'none' (do not perform any compositional adjustment), a CG percentage, a species name (e.g. 'E. coli', 'H. sapiens'), or an explicit distribution (e.g. "{'A':10, 'C':40, 'G':40, 'T':10}"). The automatic option uses a typical distribution for proteins and equiprobable distribution for everything else. sCOMP.s--weightRjtfloats=The weight of prior data. Default depends on alphabet lengths-is --first-indexRtintis5Index of first position in sequence data (default: 1)tINDEXs-ls--lowerRs"Lower bound of sequence to displays-us--upperRs"Upper bound of sequence to displays--ignore-lower-caseR^t store_truesLDisregard lower case letters and only count upper case letters in sequences.s --reverseRfsreverse sequencess --complementRhscomplement nucleic sequencess --revcompRgs$reverse complement nucleic sequencess-ss--sizeRtLOGOSIZEs=Specify a standard logo size (small, medium (default), large)s-ns--stacks-per-lineRvs@Maximum number of logo stacks per logo line. (default: %default)tCOUNTs-ts--titleRwsLogo title text.tTEXTs--labelRsA figure label, e.g. '2a's-Xs --show-xaxistbooleansYES/NOs:Display sequence numbers along x-axis? (default: %default)s-xs--xlabelR{s X-axis labels --annotateRszA comma separated list of custom stack annotations, e.g. '1,3,4,5,6,7'. Annotation list must be same length as sequences.s-Ss--yaxisRsLHeight of yaxis in units. (Default: Maximum value with uninformative prior.)tUNITs-Ys --show-yaxisRzs7Display entropy scale along y-axis? (default: %default)s-ys--ylabelRxs5Y-axis label (default depends on plot type and units)s-Es --show-endss3Label the ends of the sequence? (default: %default)s-Ps --fineprintR}s)The fine print (default: weblogo version)s --ticmarksRs-Distance between ticmarks (default: %default)s --errorbarsR~s'Display error bars? (default: %default)s--reverse-stacksRs<Draw stacks with largest letters on top? (default: %default)s-cs--color-schemeRtSCHEMEs$Specify a standard color scheme (%s)s-Cs--colorRRsCOLOR SYMBOLS DESCRIPTION tnargsisRSpecify symbol colors, e.g. --color black AG 'Purine' --color red TC 'Pyrimidine' s--default-colorRtCOLORs(Symbol color if not otherwise specified.s-Ws --stack-widths#Width of a logo stack (default: %s)tPOINTSs--aspect-ratioRs,Ratio of stack height to width (default: %s)s--boxRs(Draw boxes around symbols? (default: no)s --resolutionRi`sBitmap resolution in dots per inch (DPI). (Default: 96 DPI, except png_print, 600 DPI) Low resolution bitmaps (DPI<300) are antialiased.tDPIs --scale-widthRsScale the visible stack width by the fraction of symbols in the column? (I.e. columns with many gaps of unknowns are narrow.) (Default: yes)s--debugs=Output additional diagnostic information. (Default: %default)s--serveR!s>Start a standalone WebLogo server for creating sequence logos.s--portR#is.Listen to this local port. (Default: %default)tPORT(*RR RRR:Rtadd_option_groupt add_optionR`RRt format_namestjointformatstnamesRtstdoutRRRtlistRtkeysRRRRvRwRRyR{RzR|R}RR~RtsortRRR>R( tdefaultsR0tio_grptdata_grpt trans_grpt format_grpt color_grpt advanced_grpt server_grpt fin_choicestftcolor_scheme_choices((s9/home/psgendb/BIRCHDEV/pkg/weblogo-3.4/weblogolib/_cli.pyRsb                    -                     (,t __future__RRRLRtoptparseRRRtcorebioRt corebio.seqRRRt corebio.utilstcorebio.utils.deoptparseR t corebio._py3kR R Rt colorschemeR R R9RRRRRRRRRRRRRR7R"R%R&R(((s9/home/psgendb/BIRCHDEV/pkg/weblogo-3.4/weblogolib/_cli.pyt+s"    X % 1 3 N