ELF4 4 (444حح00 00((( Qtd/lib/ld-linux.so.2GNU%,()*#% $+ &   ! '"O>}h<%$o#@Y~$a; ;3vEw"k9H>=4}B2\l 46f (> ' H"libm.so.6logexp_Jv_RegisterClasses__gmon_start__libc.so.6putcharstrcpy__strtod_internalstdoutungetcgetenv__strtol_internalgetpidfscanffgetspcloseputsfeofmallocpopenstrrchrfprintffseekstrncmpstrncpyrealloc_IO_getcstrtoksscanfftellstrcmpsprintffclosestderrfputc__ctype_b_locfwriteexitfopen__ctype_toupper_loc_IO_stdin_used__libc_start_mainfreeGLIBC_2.0GLIBC_2.3GLIBC_2.1 ii 6ii ii ii 0*; ;000111 111 1 1 1 $1 (1,1014181<1@1D1H1L1P1T1X1\1`1d1h1 l1!p1#t1$x1%|1&1(1)1*1+Uv50%0%0h%0h%0h%1h%1h %1h(% 1h0%1h8p%1h@`%1hHP%1hP@% 1hX0%$1h` %(1hh%,1hp%01hx%41h%81h%<1h%@1h%D1h%H1h%L1h%P1hp%T1h`%X1hP%\1h@%`1h0%d1h %h1h%l1h%p1h%t1h%x1h%|1h%1h%1h %1h(%1h01^PTRhhQVh8USQ[&tX[U=;tv1ҡ1u;ÉU0tt h0ЃÐUu uXKcoa$ZU UHg+PD{UWVS }] sC==g+t]RWh05 ; uh5 ; uhX5 ;C<{utUe[^_Ív1PWh5 ;N uh5 ;8 uhX5 ;"v jbQWhp5 ; uh5 ; uh5 ;JUVSE] ut&VSPhEe[^5 ;j jhL>x=(;U@=(; ;DT T G=(;ƅe=(;,;J;Tƅ=(;53v@ 1c4(;;؋E@=(;-A@% ;QAQ p,;U@;E9(;B;@ƅ.SjItQ53h5 ;5 ;j(jhD hP P53h5 ;Sj,;;ƅ% jq;SE@9WjV,P53hT5 ;$P53h|ۍvUWVS 3<<-<6 h3<,;tD<<_5 ;j.jh8 5 ;j!jhh  e[^_á(;~1ېC ;t4‹;4Sh  ;(;|΍e[^_É,;<<g(;~11ۉFP;4Vh;tX4X;;4h ;5(;|6 h c$(;1ۉC ;t4‹;4Sh  ;(;|(;11FR;4Vh#;tX4X;;4hA ;5(;|n h Y[53h h tXZ53hb(;Ex4Q;4uh#;t040;;4hAE (;Ex}4vP;4uh;t040;;4hlE 1EvW;4uh#+;t4;;4hAE(; 9|1EP;4uh;t4;;4hE(; 9|5 ;j-jh '5 ;jjh^jUWVSLU8Eh h ju:Eh h ju#Eă u̅6}~\MMUB8-EuPh5 ;5 ;jjh aE Eȍe[^_ËU @M 82E EEEEU 4>uHMMEE}E@EM 8EEU 4>tE EEEEEvE̅19}EIMDE׋UEЀ}EMMME܋pE܀>u}E5 ;jjhb *}O1PEU 4h 5 ;=EUEvE̋8uE]̅9EU9UtgU9UUUMU8-]̅~@E؋]ċ}E7Q3VubEE9EuuME4h5 ;j5 ;j*jhL E Eȍe[^_E]̅~E]1 G;}t5{uU<:uQ3h 5 ;EG;}uRh h u h h uEȍe[^_ËMM1Ҿ1EM9IEĉ BE9UuUD:E}t5 ;jjh0  5 ;jjhJ  5 ;jjh 5 ;j&jh v]1v3VG;}u5 ;jjh0 5 ;jjhJ } 5 ;j'jhx f5 ;j$jh RjEЅF5 ;j(jh< !5 ;j)jhh  }O 1WEM 4ho 5 ; EEE8BM<t- Ruuuw WE E-E1̋]1P3hg 5 ;F9u<] ]1P3hg 5 ;zF9u.PVh 5 ;X5 ;j(jh E UWVS}] E0h| 5 ;t ;:GEƉÉv:t=xu@t 8PRh 5 ;:uƍv:teu&WRh 5 ;vߋ:t8uGtπ8tPRh 5 ;@ߋ:uʉ jrPRh 5 ;$PSh5 ;UUEu EÉRE4h5 ;5 ;j*jhL e1҃ UUEu EÉRE4h5 ;i5 ;j*jhL 1҃ UUEu EÉRE4h5 ;5 ;j*jhL 1҃ UUEu EÉRE4h5 ;5 ;j*jhL u1҃ UUMu EÉRE4h5 ;y5 ;j*jhL %1҃ UWVSu]E;u |$5 ;j4jh 1e[^_ÍvEjht Sj>jhSjjh# SjFjhSj-jh< Sj,jhlSj;jh Sj#jhSj6jh SjYjh4xYXSj ~XZSj tSjjhU Sjjh<CSjjh4Sj 9SjIjh SjJjhSj=jh(Sj SjIjhh SjKjhSjBjhSj SjGjhD SjGjhSj4jhSj SjIjh iSj nSj-jhXOSj TSjKjh5 Sj*jh#YXSj )Sj=jh Sj(jh@SjFjhl Sj jhUXZSj SjMjh SjKjhSjHjhP SjjhYXSj SjFjhr Sj!jh`XZSj fYXSj \Sjjh= Sjjhc+SjjhSj !Sj jh(Sj Sjjh} SjjhSj,jhL Sjjh+Sjjh Sj.jh|Sjjh+ SjEjhsSj#jhdSj iSjEjhJ Sj#jh8XZSj >SjjhSj $SjAjh` SjjhYXSj SjHjh SjAjhSjFjh4 SjFjh|SjIjh SjFjhSjFjhXw SjHjheSjEjhV SjJjh4 DSjEjh 5 SjCjh #Sj*jh !Sj Sjvtǃ( #f f tt UWVS ]+]hp+h+jCPƃx1D@9}E ~.1ۋE 1ҋ+EƋ}׋} B9U uC9uԋE0e[^_ÍvUSE] tә؅u[UWVSEDU @}M [-%s] [-%s <%s>] %d%cis not an integer. integer at its end. is not a positive integer. point number at its end. PAGERwmorers more than once on the command line. does not match any of the legal options. An option in the "exclusive" array does not appear in the "options" array. A required option is mutually exclusive with a non-required option. and "-%s" options are mutually exclusive. One of the following mutually exclusive options is required on the command line: The required option "-%s" is missing The command line is missing an integer at its end. "%s" (from item %d on the command line) The command line is missing a positive The command line is missing a non-negative is not a non-negative integer. The command line is missing a floating is not a floating point number. The command line is missing a positive floating is not a positive floating point number. The command line is missing a non-negative floating is not a non-negative floating point number. The command line is missing a string of characters at its end ########################### # README FOR MATALIGN-V2A # Author: Ting Wang and Gary Stormo Department of Genetics Campus Box 8232 St. Louis, MO 63110 stormo@ural.wustl.edu twang@ural.wustl.edu MatrixAligner-v2b # BASIC OPTIONS # -f1 -f2 [-h ] # GENERAL INFORMATION # matrices. # COMMAND LINE OPTIONS # 1) Required options -f1 filename -f2 filename 2) Alphabet options information. letter letter normalization letter:complement A:T C:G A:T 0.3 C:G 0.2 4) Algorithm options # Example # ***** PID: 4450 ***** Algorithm options: Compare both orientations. Global alignment. Matrix 1 type: count matrix Matrix 2 type: count matrix Matrices Alignment: Comparison scores: ALLR = 8.3894 Distance = 1.9606 E_value = 3.01e-07 P_value = 3.01e-07 NEW MATRIX: number of sequences = 40 width = 7 Consensus: GCGGGGc Copyright 2002--2005 Ting Wang and Gary Stormo May be copied for noncommercial purposes. Washington University in St. Louis -t1 -t2 [-n ] [-c0 ] [-c1 ] [-g ] [-a ] [-A ] [-CS ] Matrix Aligner is a program to compare two positional specific matrices. The precursor of this program is "CompareTwo". Matalign-v2b made several improvements over the previous establishment of the program. Two input PSSMs should be formatted according to consensus output. They can be either count matrices, or frequency matrices. If a frequency matrix is used, a user-defined "total count" is assumed (default = 10). The scoring function between two positions of the two matrices is ALLR statistic. The alignment algorithm can be either "local alignment" or "global alignment". Both options do not allow gaps. The comparison between the two matrices results in the following output: 1) An alignment between two input matrices. 2) An ALLR score of the comparison. In general, the higher the ALLR score, the more similar are the two matrices. 3) A distance score. Distance between A and B is defined as: ALLR(A,A) + ALLR(B,B) - 2xALLR(A,B). In general, the smaller the distance, the more similar are the two 4) P-value and E-value of the observed ALLR score. These significance values are calculated based on Karlin-Altschul statistics. The meaning of the p-value is, given two random PSSMs, the probability of observing an equal or higher ALLR score. 5) The aligned parts of the two matrices will be merged into one new matrix and printed as output. 0) -h: print these directions. "filename" contains the one matrix. "filename" contains the second matrix. -t1 integer: the type of matrix 1. 0 means count matrix, 1 means frequency matrix. -t2 integer: the type of matrix 1. 0 means count matrix, 1 means -a filename: file containing the alphabet and normalization Each line contains a letter (a symbol in the alphabet) followed by an optional normalization number (default: 1.0). The normalization is based on the relative prior probabilities of the letters. For nucleic acids, this might be be the genomic frequency of the bases; however, if the "-d" option is not used, the frequencies observed in your own sequence data are used. In nucleic acid alphabets, a letter and its complement appear on the same line, separated by a colon (a letter can be its own complement, e.g. when using a dimer alphabet). Complementary letters may use the same normalization number. Only the standard 26 letters are permissible; however, when the "-CS" option is used, the alphabet is case sensitive so that a total of 52 different characters are possible. POSSIBLE LINE FORMATS WITHOUT COMPLEMENTARY LETTERS: POSSIBLE LINE FORMATS WITH COMPLEMENTARY LETTERS: letter:complement normalization letter:complement normalization:complement's_normalization Example alphabet file 1: Example alphabet file 2: 3) Options for handling the complement of the matrices --- -c0: ignore the complement -c1: compare both orientation (the default option) -n: if a frequency matrix is compared, a total count is assumed. "n" is the user-defined total count. Default is 10. -g: global alignment. Default is local alignment. (bifrost)[3:11pm]src_v2b 207>>more matrix1 A | 1 0 0 0 0 0 0 C | 0 14 2 0 0 0 9 G | 15 2 14 16 16 16 2 T | 0 0 0 0 0 0 5 (bifrost)[3:23pm]src_v2b 208>>more matrix2 A | 0 0 1 0 1 0 0 C | 3 22 0 0 0 1 15 G | 18 1 23 24 23 23 1 T | 3 1 0 0 0 0 8 (bifrost)[3:23pm]src_v2b 209>>./matalign-v2b -f1 matrix1 -f2 matrix2 -t1 0 -t2 0 -g COMMAND LINE: ./matalign-v2b -f1 matrix1 -f2 matrix2 -t1 0 -t2 0 -g MATRIX: A | 1 0 0 0 0 0 0 C | 0 14 2 0 0 0 9 G | 15 2 14 16 16 16 2 T | 0 0 0 0 0 0 5 Consensus: G C G G G G c MATCH: | | | | | | | MATRIX: A | 0 0 1 0 1 0 0 C | 3 22 0 0 0 1 15 G | 18 1 23 24 23 23 1 T | 3 1 0 0 0 0 8 crude information = 9.7450 (bits) A | 1 0 1 0 1 0 0 C | 3 36 2 0 0 1 24 G | 33 3 37 40 39 39 3 T | 3 1 0 0 0 0 13 alloc_column()columnfrequencymake_mirror_profile()mirrordcodeError in opening file %s read_int_matrix()read_matrix()int_matrix[]int_matrixdouble_matrix[]double_matrixinit_profiles()hsptrace_back()trace_back_g()hsp->matchDP_functionDPDP_function()DP[]score_function()STscoreSTscore[]free_scores()merge_HSP()new_profileNEW MATRIX:number of sequences = %d width = %d %c | %3d%3d |Consensus: %6sMATRIX: %5d %5cMatrices Alignment:MATCH: Comparison scores:ALLR = %.4f Full Distance = %.4f Aligned Distance = %.4f E_value = %.4g P_value = %.4g B?9B.?crude information = %.4f (bits) ffffff?ffffff?find_score_spectrum()spectrumkarlin_lambda_K_H()p0P0karlin_altschul()smsm[]Error: The lowest score in the substitution matrix must be negative.Error: The highest score in the substitution matrix must be positive.Error: A positive score is impossible with the given substitution matrix and residue composition.Error: Negative probabilities are disallowed.Error: Probabilities sum to %.5f and will be normalized to 1. Error: Invalid (non-negative) expected score: %#0.21gError: Value for K may be too large due to insufficient iteration.St$?9m4?h㈵>6  HHT  08ćoLoo0*:JZjzʉډ *:JZjzʊڊ *:JZjz 0d %) %)%)%))<',$;*@, ;CY,;T\,;(,11 ع1, ع;P _<{\3%)I# lSjI$m kRs4Xl1Ds4lVD1$<|#<}p!/p/!o44FF8z{ mXlzl lV{Sk1D}p&jRD1|o'qpptts<Yqpptts;Xpporrri]pporrsh[trTvWrXqtrvTrWXqtrWrTvXqtrrWvTXqrsXXXX;vssqqqqvl<6;6i68h68YX][ׁs?(\µ?Af?Af?<,Ԛ?ZӼ?U؟?Ǻ?ea?_vO?St$? ^)?<,Ԛ?QI?`"?ATCGMKRYWSatcgnGCC: (GNU) 4.0.0 20050525 (Red Hat 4.0.0-9)GCC: (GNU) 4.0.0 20050525 (Red Hat 4.0.0-9)GCC: (GNU) 4.0.0 20050519 (Red Hat 4.0.0-8)GCC: (GNU) 4.0.0 20050519 (Red Hat 4.0.0-8)GCC: (GNU) 4.0.0 20050519 (Red Hat 4.0.0-8)GCC: (GNU) 4.0.0 20050519 (Red Hat 4.0.0-8)GCC: (GNU) 4.0.0 20050519 (Red Hat 4.0.0-8)GCC: (GNU) 4.0.0 20050519 (Red Hat 4.0.0-8)GCC: (GNU) 4.0.0 20050519 (Red Hat 4.0.0-8)GCC: (GNU) 4.0.0 20050519 (Red Hat 4.0.0-8)GCC: (GNU) 4.0.0 20050519 (Red Hat 4.0.0-8)GCC: (GNU) 4.0.0 20050519 (Red Hat 4.0.0-8)GCC: (GNU) 4.0.0 20050519 (Red Hat 4.0.0-8)GCC: (GNU) 4.0.0 20050519 (Red Hat 4.0.0-8)GCC: (GNU) 4.0.0 20050519 (Red Hat 4.0.0-8)GCC: (GNU) 4.0.0 20050519 (Red Hat 4.0.0-8)GCC: (GNU) 4.0.0 20050519 (Red Hat 4.0.0-8)GCC: (GNU) 4.0.0 20050525 (Red Hat 4.0.0-9).symtab.strtab.shstrtab.interp.note.ABI-tag.hash.dynsym.dynstr.gnu.version.gnu.version_r.rel.dyn.rel.plt.init.text.fini.rodata.eh_frame.ctors.dtors.jcr.dynamic.got.got.plt.data.bss.comment#(( 1HHL7 ?TTGoXToLL`c l ć8 up { sHHhhl.-ԭ0000001h ;*2l`K  (HTLć    H h-0000001; 0*080E;T1[܋ q }0 0-0  l t ,k Ԏ &1;Na1,n2y 3D ܓL ( X p t   h < , =J_m@y 0; ;xP 8;<  3\ > ;5 $ 2(;9}I@;P  al t<06s %@F hH;;, $P; #(G 6H D0U, _؏% j8 z1   0 ,  ,; 9 3B +Y< L;R;YX;_o u@ ,S  P ~TV  $(P ȺP `;[ ' /h;4 >h Ko [;m ;-  :xȻ  ;  3z 3p;  ,\ >> F<N NEkW x9 < t" wt;00 O ; (r  ` 82 !>x;@0Q:Wfyl 1 @ 9H >   1 p  00  9 I L V da g  q  } ; 4 ; 0  ; 1 $;  2 0 l ;# t 6 H 4[ fm 1z   > 0 @3   ; TE   8 . P <  G P P عP X l d s "call_gmon_startcrtstuff.c__CTOR_LIST____DTOR_LIST____JCR_LIST__completed.4583p.4582__do_global_dtors_auxframe_dummy__CTOR_END____DTOR_END____FRAME_END____JCR_END____do_global_ctors_auxmain.ctw_alloc_helper.cdetermine_gross_sizedetermine_gross_arraymalloc_checkrealloc_checktw_error.ctw_file.cmatalign_options.cmatalign_command.cST_exclusiveST_optionsalpha.cSTalpha_flagALPHA_read_letterALPHA_read_numberALPHA_read_colonALPHA_allocALPHA_reallocALPHA_truncateALPHA_get_letterALPHA_get_doubleALPHA_get_colonALPHA_get_intALPHA_read_asciiparse-line.cmatalign_direction.cread_matrix.ccompare_profile.cmatalign_allr.cmatalign_stats.cvalE_valuefile2pl_Int_3Ngreatest_common_denominatorAlpha_filepl_Alpha_ac__strtod_internal@@GLIBC_2.0type1get_ALLR_E_valueFile_get_pathA_sizefeof@@GLIBC_2.0Lambdacompare_profilescalculate_distanceungetc@@GLIBC_2.0getpid@@GLIBC_2.0_DYNAMICSM1print_alphastrcmp@@GLIBC_2.0usage_fp_hwPIDAparse_linefprintf@@GLIBC_2.0Kgetenv@@GLIBC_2.0init_profilespclose@@GLIBC_2.1ALLR_functioncheck_options__fini_array_endprint_HSPbug_reportread_int_matrix__dso_handlememory_integrity_check__libc_csu_finialloc_columnpl_Nn_Inte_to_p_helper_1Comp_flagputchar@@GLIBC_2.0find_score_spectrumprint_indentAsciicount_2_frequencyftell@@GLIBC_2.0karlin_altschultype2A_compDist1puts@@GLIBC_2.0_initcompare_2_profilespl_Stringpopen@@GLIBC_2.1pl_Alpha_aftrace_back_gmalloc@@GLIBC_2.0to_nearest_intDP_function_gfscanf@@GLIBC_2.0pl_Int_2pl_Int_4Dist2frequency_2_counttw_freeALLRfind_linepl_Nn_Doublee_to_p_helper_2stdout@@GLIBC_2.0stderr@@GLIBC_2.0File_existsFile_put_pathDcode_fpl_Intadjust_alphabetfseek@@GLIBC_2.0Globalmake_mirror_profile_startfgets@@GLIBC_2.0print_directionsCase_sensitivehsppl_P_Doublekarlin_lambda_K_Hpl_Charpl_Help__strtol_internal@@GLIBC_2.0command_linetw_reallocpl_Alpha_ifcomplete_profilestrncmp@@GLIBC_2.0profile__fini_array_startscore_functionfputc@@GLIBC_2.0__libc_csu_init__bss_startcal_IDP_functionmain__libc_start_main@@GLIBC_2.0H__init_array_endDcodeexp@@GLIBC_2.0realloc@@GLIBC_2.0cal_infodata_startmerge_HSPprintf@@GLIBC_2.0_finifclose@@GLIBC_2.1__ctype_toupper_loc@@GLIBC_2.3strrchr@@GLIBC_2.0tw_malloctext_directionsNum_fpl_P_Int__preinit_array_endpl_ICharexit@@GLIBC_2.0print_MATRIXlog_poisson_unittw_callocread_matrixMsscanf@@GLIBC_2.0_edata_GLOBAL_OFFSET_TABLE_free@@GLIBC_2.0_endComp_statusfile1strncpy@@GLIBC_2.0fopen@@GLIBC_2.1__init_array_start_IO_stdin_usedP_valueread_double_matrixstrtok@@GLIBC_2.0sprintf@@GLIBC_2.0fwrite@@GLIBC_2.0__data_startmap_column_dcode_IO_getc@@GLIBC_2.0_Jv_RegisterClasses__ctype_b_loc@@GLIBC_2.3__preinit_array_startSM2e_value_to_p_valuelog@@GLIBC_2.0Ptw_recallocpl_Doubleprint_profileprint_optionstrace_backpl_NBoolpl_Boolfree_scores__gmon_start__strcpy@@GLIBC_2.0