ELFT44 (444`` ``((( Qtd/lib/ld-linux.so.2GNU%+'()"$ #* %    &!}Vv<%$]#@Y~}$Ok k3dEwY96>+4k02J4$f > ' H"libm.so.6logexp_Jv_RegisterClasses__gmon_start__libc.so.6putcharstrcpystdoutungetcgetenv__strtol_internalgetpidfscanffgetspcloseputsfeofmallocpopenstrrchrfprintffseekstrncmpstrncpyrealloc_IO_getcstrtoksscanfftellstrcmpsprintffclosestderrfputc__ctype_b_locfwriteexitfopen__ctype_toupper_loc_IO_stdin_used__libc_start_mainfreeGLIBC_2.0GLIBC_2.3GLIBC_2.1 ii p6ii zii ii p`)kk```aaa aaa a a a $a (a,a0a4a8aL>x=8kU@=8k {DT T G=8kƅe=8kuHMMEE}E@EM 8EEU 4>tE EEEEEvE̅19}EIMDE׋UEЀ}EMMME܋pE܀>u}E5kjjh}O1PEU 4hl5k)EUEvE̋8uE]̅9EU9UtgU9UUUMU8-]̅~@E؋]ċ}E7Q3VuNEE9EuuME4h5kV5kj*jhE Eȍe[^_E]̅~E]1 G;}t5{uU<:uQ3h5kEG;}uRhHhUu hHhZuEȍe[^_ËMM1Ҿ1EM9IEĉ BE9UuUD:E}t5kjjh5kjjh 5kjjh$5kj&jhDv]1v3VG;}u5kjjh}5kjjhi 5kj'jhR5kj$jh>jEЅF5kj(jh 5kj)jh}O 1WEM 4h5k EEE8BM<t- Ruuuw WE E-E1̋]1P3h5kF9u<] ]1P3h5kfF9u.PVh~5kD5kj(jh|E UWVS}] E0h5kt ;:GEƉÉv:t=xu@t 8PRh5k:uƍv:teu&WRh5kbߋ:t8uGtπ8tPRh5k,ߋ:uʉ j^PRh5k$PSh5kUUEu EÉRE4h5k5kj*jhQ1҃ UUEu EÉRE4h5kU5kj*jh1҃ UUEu EÉRE4h5k5kj*jh1҃ UUEu EÉRE4h5k5kj*jha1҃ UUMu EÉRE4h5ke5kj*jh1҃ UWVSu]E;u |$5kj4jh1e[^_ÍvE SjJjh,SjOjh SjIjh8 Sjjh7Sj SjGjh SjHjhSjGjh  SjGjh` SjIjh  SjFjh SjjhRSj YXSj zSjjhT[ SjjhnISjjhT:Sj ?Sjjh Sj,jh Sj jh3,Sjjh Sjjh Sj jh3 SjjhSjjh۾Sj SjSjpjh6/ Sjpjh7SjZjht7 Sj,jh7Sj,jh8 Sj,jh08ۼSj,jh`8̼ SjOjh8躼SjIjh8諼 Sjjh虼Sjjh芼 Sj>jh,9xSjjhi SjjhWSjjhH Sjjh6Sj-jhl9' Sj$jh9Shjh9 Shjh`:Sjzjh:߻ Shjhx;ʻShjh;踻 Shjh<裻Shjh=葻 Sjjjh=Sjrjh=p Sjrjhh>^Sjrjh>O SjrjhP?=Sjzjh?. Sjzjh@@Sjzjh@ Sjzjh8AShjhA Shjh@BԺShjhBº ShjhPC譺ShjhC蛺 ShjhhD膺SjcjhDw Sj]jhXEeSjjhV SjjhDSjXjhE5 Sjjh#Sjjh SjjhSjjh Sj;jhFSj$jhPFҹ Sj0jhxFSj+jhF豹 Sj+jhF蟹Sj*jhG萹 Sj*jh0G~Sj*jh\Go Sj*jhG]Sj*jhGN Sj+jhG Sj+jh`M,Sj+jhM Sj+jhM Sj*jhM Sj*jhNSj*jhu 1e[^_h9OhKOj1Q蕹E1E1I9v3v# ыDP tLE8uC1I9wЋUE h9OhKOjE@PRع e[^_ÐUE @ECUWVS ]hSOhdOjj E^_hOS9ERuhh@k(tV h@k1эY|/t' hSOhKOjAPR"/D E 1۾ E|Puhh@k観t=C h@kIǍF9|Ń hSOhdOjVu许E 뢃 u׳EUEe[^_QhSOhdOuȶ ShrO5kx$輳UWVSPExpPUhOhOj苷EhOhOjuTEă hOhOju:EȃEЅzhOhOjuEhOhOjuEԃ EЅMȉM܋EPUlj}ƉuEEZMىM؋UBMuԋU؉ M1űU؉EEmmmmM9M8k~5{uEUċEMIMu܋>E@EUډUؐuLuuDu΋TDTC98k/uUMĉJuȉrEЉBMԉJủr EBЍe[^_1hOhOj58k薵MĉhOhOj58kyUȉC 9]uvtǃ( #f f tt UWVS }hRhv j58kƋ 8k~;11҉@9uR${1ҍvB9u؉e[^_ÐUWVSHuwE^u kEY[Pu]XZuu]؋58k]U]1PhRhRuL hRhRu7EM]EM]EES4$UЍe[^_Ð111ɉ} } B9uS$] $$]sUWVS ]+]hRhRjCP蘕ƃ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. PAGERmorers 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-V4a # Author: Ting Wang and Gary Stormo Department of Genetics Campus Box 8232 St. Louis, MO 63110 stormo@ural.wustl.edu twang@ural.wustl.edu NOTES ON THIS UPDATE: a 50~100x gain in speed. would be more consistent. # BASIC OPTIONS # Usage: ./matalign-v4a [-h ] # GENERAL INFORMATION # E_value: # COMMAND LINE OPTIONS # 1) Required options -f1 filename a matrix file name. -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 # GCN4_01.matrix GCN4_C.matrix HSF_01.matrix HSF_02.matrix HSF_03.matrix HSF_04.matrix HSF_05.matrix ***** PID: 19688 ***** Algorithm options: Compare both orientations. Global alignment. Output option: ***** PID: 32161 ***** Copyright 2002--2006 Ting Wang and Gary Stormo May be copied for noncommercial purposes. Washington University in St. Louis MatrixAligner-v4a, or "MegaMatAlign". MAIN DIFFERENCE FROM V2: batch process of many matrix comparisons. This version allows comparison among many matrices, and fix some bugs. Credit should also be given to Alok Saldanha [alok@caltech.edu] who did lots of initial experimenting with this idea. Strategy: user now provides one or two files that contains file names of matrices to compare. The first line of the file always contains the directory information of the matrices (therefore all matrices should be under the same directory). If one list is provided, all pair-wise comparison will be performed; if two lists are provided, all matrices in list file 1 will be compared to all matrices in list file 2. I also did some software engineering to speed up the comparison. In this version, all matrices should have the format of a count vector. I will allow more variations to the format in later modifications. To illustrate the speed gain, I compared all TRANSFAC matrices against themselves (636 matrices, 201930 comparisons). It took MegaMatAlign about 30 seconds to complete the comparison, while if I run a perl script that calls MatAlign that many times, it took almost an hour. So, there is at least I should also point out that another trick for speeding up the process is to calculate a ALLR lookup table upfront. However, since pseudocount treatment is different between the lookup table and real calculation, the ALLR scores are slightly different. I don't think this will cause significant different, but have to watch out. I may turn this off. Then it would run for twice amount the time, but the logrithm calculation Use the "-h" option for detailed directions -f1 [-f2 ] [-c0 ] [-c1 ] [-g ] [-ez ] [-a ] [-i ] [-A ] [-CS ] Matrix Aligner is a program to compare two positional specific matrices. The precursor of this program is "CompareTwo". Matalign-v4a made several improvements over the previous establishment of the program. "MegaMatAlign" enables comparison of many matrices. The input file contains matrix names (and their location, i.e. directory information). If there is one input list file, MegaMatAlign pair-wise compares all matrices listed in this file; if there are two input list files, then MegaMatAlign compares every matrix listed in file 1 to every matrix listed in file 2. Input PSSMs should be formatted according to consensus output. They should be count matrix like this: (bifrost)[11:06am]src_v4a 574>>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 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 program generates a data table as output. The table contains the following columns, tab-delimited: Matrix_1: Name of the first matrix in comparison Consensus: Consensus pattern of the first matrix Matrix_2: Name of the second matrix in comparison Consensus: Consensus pattern of the second matrix ALLR: ALLR score of the comparison. In general, the higher the ALLR score, the more similar are the two matrices. Overlap: The overlap length (bp) of the aligned part Distance: 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 matrices. Aligned_Dist: Similar defined as Distance, but only calculated for the aligned part. Shared_Consensus: The aligned part of the two matrices are merged, and a new matrix generated. This is the consensus of the new matrix. P_value: 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. If "-ez" is defined at command line, the output table will only contain Matrix_1, Matrix_2, ALLR, and Overlap for simplicity. 0) -h: print these directions. "filename" contains information about matrices: the first line of the file contains directory, after which every line contains "filename" contains the second list of matrices. -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) -g: global alignment. Default is local alignment. -ez: only generate ALLR score and Overlap for each comparison. (bifrost)[11:32am]src_v4a 598>>more sample_list /home3/twang/CompareTwo/src_v4a/sample (bifrost)[11:32am]src_v4a 599>>more /home3/twang/CompareTwo/src_v4a/sample/GCN4_01.matrix A | 4 2 7 2 3 5 5 4 6 19 0 0 43 0 2 1 43 1 12 4 12 6 7 4 5 8 6 C | 6 7 6 7 6 4 11 12 9 6 0 0 0 43 0 42 0 11 16 14 10 9 10 8 10 7 7 G | 2 5 2 9 6 10 7 9 11 15 0 43 0 0 0 0 0 5 2 13 4 11 4 11 6 3 4 T | 2 1 3 2 8 6 4 5 11 2 43 0 0 0 41 0 0 26 8 3 8 7 10 5 6 7 6 (bifrost)[11:33am]src_v4a 600>>more /home3/twang/CompareTwo/src_v4a/sample/GCN4_C.matrix A | 13 25 21 0 0 38 0 0 0 3 C | 12 0 0 0 0 0 30 0 38 1 G | 8 12 12 0 38 0 0 0 0 0 T | 4 0 4 38 0 0 8 38 0 1 (bifrost)[11:33am]src_v4a 601>>./matalign-v4a -f1 sample_list -a ~/alphabet -g COMMAND LINE: ./matalign-v4a -f1 sample_list -a /home3/twang/alphabet -g Pair-wise comparison among matrices specified in sample_list. Complete output, including all statistics. Number of pair-wise comparisons: 21 Matrix_1 Consensus Matrix_2 Consensus ALLR Overlap Distance Aligned_Dist Shared_Consensus E_value P_value GCN4_01.matrix nnnnnnnnnRTGACTCAtnSnnnnnnn GCN4_C.matrix naaTGACTCa 9.362 10 8.874 7.115 nnaTGACTCA 2.165e-07 2.165e-07 GCN4_01.matrix nnnnnnnnnRTGACTCAtnSnnnnnnn HSF_01.matrix aGAAn -0.018 1 20.037 0.138 n 36.51 1 GCN4_01.matrix nnnnnnnnnRTGACTCAtnSnnnnnnn HSF_02.matrix aGAAnaGAAnaGAAn -0.018 1 29.985 0.138 n 109.5 1 GCN4_01.matrix nnnnnnnnnRTGACTCAtnSnnnnnnn HSF_03.matrix aGAAnaGAAnnTTCt -0.215 1 30.379 0.743 n 165.4 1 GCN4_01.matrix nnnnnnnnnRTGACTCAtnSnnnnnnn HSF_04.matrix aGAAnnTTCtaGAAn -0.018 1 29.985 0.138 n 109.5 1 GCN4_01.matrix nnnnnnnnnRTGACTCAtnSnnnnnnn HSF_05.matrix nTTCtaGAAnaGAAn -0.018 1 29.985 0.138 n 109.5 1 GCN4_C.matrix naaTGACTCa HSF_01.matrix aGAAn -0.044 1 17.633 0.210 n 14.3 1 GCN4_C.matrix naaTGACTCa HSF_02.matrix aGAAnaGAAnaGAAn -0.044 1 27.581 0.210 n 42.89 1 GCN4_C.matrix naaTGACTCa HSF_03.matrix aGAAnaGAAnnTTCt -0.225 1 27.943 1.079 a 62.66 1 GCN4_C.matrix naaTGACTCa HSF_04.matrix aGAAnnTTCtaGAAn -0.044 1 27.581 0.210 n 42.89 1 GCN4_C.matrix naaTGACTCa HSF_05.matrix nTTCtaGAAnaGAAn -0.044 1 27.581 0.210 n 42.89 1 HSF_01.matrix aGAAn HSF_02.matrix aGAAnaGAAnaGAAn 4.974 5 9.948 -0.000 aGAAn 0.000587 0.0005868 HSF_01.matrix aGAAn HSF_03.matrix aGAAnaGAAnnTTCt 4.974 5 9.948 -0.000 aGAAn 0.000587 0.0005868 HSF_01.matrix aGAAn HSF_04.matrix aGAAnnTTCtaGAAn 4.974 5 9.948 -0.000 aGAAn 0.000587 0.0005868 HSF_01.matrix aGAAn HSF_05.matrix nTTCtaGAAnaGAAn 4.974 5 9.948 -0.000 aGAAn 0.000587 0.0005868 HSF_02.matrix aGAAnaGAAnaGAAn HSF_03.matrix aGAAnaGAAnnTTCt 9.948 10 9.948 -0.000 aGAAnaGAAn 5.288e-08 5.288e-08 HSF_02.matrix aGAAnaGAAnaGAAn HSF_04.matrix aGAAnnTTCtaGAAn 4.974 5 19.896 -0.000 aGAAn 0.001761 0.001759 HSF_02.matrix aGAAnaGAAnaGAAn HSF_05.matrix nTTCtaGAAnaGAAn 9.948 10 9.948 -0.000 aGAAnaGAAn 5.288e-08 5.288e-08 HSF_03.matrix aGAAnaGAAnnTTCt HSF_04.matrix aGAAnnTTCtaGAAn 9.948 10 9.948 -0.000 aGAAnnTTCt 5.288e-08 5.288e-08 HSF_03.matrix aGAAnaGAAnnTTCt HSF_05.matrix nTTCtaGAAnaGAAn 9.948 10 9.948 -0.000 aGAAnaGAAn 5.288e-08 5.288e-08 HSF_04.matrix aGAAnnTTCtaGAAn HSF_05.matrix nTTCtaGAAnaGAAn 9.948 10 9.948 -0.000 nTTCtaGAAn 5.288e-08 5.288e-08 (bifrost)[11:33am]src_v4a 603>>./matalign-v4a -f1 sample_list -f2 sample_list -a ~/alphabet -g -ez COMMAND LINE: ./matalign-v4a -f1 sample_list -f2 sample_list -a /home3/twang/alphabet -g -ez Compare all matrices specified in sample_list to all matrices specified in sample_list. Simple output, only names of matrices, ALLR and Overlap. Number of pair-wise comparisons: 49 Matrix_1 Matrix_2 ALLR Overlap GCN4_01.matrix GCN4_01.matrix 15.027 27 GCN4_01.matrix GCN4_C.matrix 9.362 10 GCN4_01.matrix HSF_01.matrix -0.018 1 GCN4_01.matrix HSF_02.matrix -0.018 1 GCN4_01.matrix HSF_03.matrix -0.215 1 GCN4_01.matrix HSF_04.matrix -0.018 1 GCN4_01.matrix HSF_05.matrix -0.018 1 GCN4_C.matrix GCN4_01.matrix 9.362 10 GCN4_C.matrix GCN4_C.matrix 12.570 10 GCN4_C.matrix HSF_01.matrix -0.044 1 GCN4_C.matrix HSF_02.matrix -0.044 1 GCN4_C.matrix HSF_03.matrix -0.225 1 GCN4_C.matrix HSF_04.matrix -0.044 1 GCN4_C.matrix HSF_05.matrix -0.044 1 HSF_01.matrix GCN4_01.matrix -0.018 1 HSF_01.matrix GCN4_C.matrix -0.044 1 HSF_01.matrix HSF_01.matrix 4.974 5 HSF_01.matrix HSF_02.matrix 4.974 5 HSF_01.matrix HSF_03.matrix 4.974 5 HSF_01.matrix HSF_04.matrix 4.974 5 HSF_01.matrix HSF_05.matrix 4.974 5 HSF_02.matrix GCN4_01.matrix -0.018 1 HSF_02.matrix GCN4_C.matrix -0.044 1 HSF_02.matrix HSF_01.matrix 4.974 5 HSF_02.matrix HSF_02.matrix 14.922 15 HSF_02.matrix HSF_03.matrix 9.948 10 HSF_02.matrix HSF_04.matrix 4.974 5 HSF_02.matrix HSF_05.matrix 9.948 10 HSF_03.matrix GCN4_01.matrix -0.215 1 HSF_03.matrix GCN4_C.matrix -0.225 1 HSF_03.matrix HSF_01.matrix 4.974 5 HSF_03.matrix HSF_02.matrix 9.948 10 HSF_03.matrix HSF_03.matrix 14.922 15 HSF_03.matrix HSF_04.matrix 9.948 10 HSF_03.matrix HSF_05.matrix 9.948 10 HSF_04.matrix GCN4_01.matrix -0.018 1 HSF_04.matrix GCN4_C.matrix -0.044 1 HSF_04.matrix HSF_01.matrix 4.974 5 HSF_04.matrix HSF_02.matrix 4.974 5 HSF_04.matrix HSF_03.matrix 9.948 10 HSF_04.matrix HSF_04.matrix 14.922 15 HSF_04.matrix HSF_05.matrix 9.948 10 HSF_05.matrix GCN4_01.matrix -0.018 1 HSF_05.matrix GCN4_C.matrix -0.044 1 HSF_05.matrix HSF_01.matrix 4.974 5 HSF_05.matrix HSF_02.matrix 9.948 10 HSF_05.matrix HSF_03.matrix 9.948 10 HSF_05.matrix HSF_04.matrix 9.948 10 HSF_05.matrix HSF_05.matrix 14.922 15 (bifrost)[11:34am]src_v4a 604>> read_1st_string()string1read_list_file()first_stringsCannot open file %s. Exit. make_revcomp_profile()mirrorint_matrixf_matrixint_matrix[]f_matrix[]dcodenum_seqError in opening file %s read_int_matrix()read_matrix()read_matrix_files()profiles?Warning: Matrix file %s doesn't exist. DP_functionDPDP_function()DP[]DP_function_g()trace_back()hspscore_function()STscoreSTscore[]merge_HSP()new_profileNewfree_PROFILE()Matrix_1%s Matrix_2ALLROverlapDistanceAligned_DistShared_ConsensusE_valueP_value%.3f %d %.4g free_HSP()free_scores()compare_profiles()Hsp_QNumber of pair-wise comparisons: %d Warning: HSP numbers not consistent -- %d vs %d 9B.?ffffff?ffffff?init_allr_lookup_table()C_LN_CC_LN_C[]C_LN_C[][]C_LN_NC_LN_N[]C_LN_N[][]C_LN_PC_LN_P[]cal_f()ALLR_score()f_1f_2find_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 ̈ pH@  `0o$ooΆ` *:JZjzʉډ *:JZjzʊڊ *:J `/2lQ579lQ;7lQlQ_>O(kRm$kp/\aT2a|X4k0k5790 ;\b0 lQ eI# 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#(( 1HHH7 ?@@GoΆVTo$$`c l 0 üpp{TT zpppOU``````aH k *LO  (H@Ά$ ̈  T p U``````akx `*`8`EkTa[ qЋ }` `U`H 0 8 k  &1;Naa,nbyb L  X ؕp H   <  ,ؠ =J_m@kvH{D{@{,y $kLP P{(  \ b 0 5 &ď 48k;}KX{R c v B <`,kfs %F `{{ $h{##5 E! S  a`r% } a kP  p{@    } N E\  `9 0k w` (O 8jDr J$` V7 [x{z`ja P1 9p > $ | . `) > D G `[  d t a \  { 4 j `  { PG   a k (k" 5 k; 2L `_ {g v  k 4 f a t  > `, c0 Dw 6  I X {Z E f  p $ ~ |+ n \P P l "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.cSTbuffercompare_profile.cmatalign_allr.cC_LN_CC_LN_NC_LN_Pmatalign_stats.cvalfile2pl_Int_3Ngreatest_common_denominatorprint_HSP_infoAlpha_filepl_Alpha_acget_ALLR_E_valueFile_get_pathA_sizefeof@@GLIBC_2.0Lambdacompare_profilescalculate_distanceprint_consensusungetc@@GLIBC_2.0getpid@@GLIBC_2.0_DYNAMICList_numSM1print_alphastrcmp@@GLIBC_2.0usage_fp_hwPIDAparse_linefprintf@@GLIBC_2.0Kgetenv@@GLIBC_2.0pclose@@GLIBC_2.1read_1st_stringALLR_functioncheck_options__fini_array_endbug_reportread_int_matrix__dso_handleinit_allr_lookup_tableMat_names_2memory_integrity_checkProfiles_1__libc_csu_finipl_Nn_Inte_to_p_helper_1Comp_flagputchar@@GLIBC_2.0find_score_spectrumget_hsp_E_valueAsciiProfiles_2ftell@@GLIBC_2.0karlin_altschulA_compputs@@GLIBC_2.0_initcompare_2_profilespl_Stringpopen@@GLIBC_2.1pl_Alpha_aftrace_back_gmalloc@@GLIBC_2.0DP_function_gfscanf@@GLIBC_2.0pl_Int_2pl_Int_4Table_sizepath1tw_freefind_linepl_Nn_Doublee_to_p_helper_2stdout@@GLIBC_2.0stderr@@GLIBC_2.0File_existsFile_put_pathDcode_fpl_Intadjust_alphabetread_list_fileHsp_QMat_num_1fseek@@GLIBC_2.0Global_startfgets@@GLIBC_2.0read_matrix_filesprint_directionsCase_sensitivemake_revcomp_profilepl_P_Doublekarlin_lambda_K_Hpl_Charread_matricespl_Help__strtol_internal@@GLIBC_2.0command_lineALLR_scoretw_reallocEZpl_Alpha_ifstrncmp@@GLIBC_2.0__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.0data_startmerge_HSPprintf@@GLIBC_2.0_finifclose@@GLIBC_2.1__ctype_toupper_loc@@GLIBC_2.3strrchr@@GLIBC_2.0tw_malloctext_directionspl_P_Int__preinit_array_endpl_ICharexit@@GLIBC_2.0log_poisson_unittw_callocread_matrixMsscanf@@GLIBC_2.0_edata_GLOBAL_OFFSET_TABLE_free@@GLIBC_2.0_endcalculate_ALLR_doublefree_PROFILEComp_statusMat_names_1file1strncpy@@GLIBC_2.0path2fopen@@GLIBC_2.1__init_array_startHsp_num_IO_stdin_usedstrtok@@GLIBC_2.0Mat_num_2sprintf@@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_startSM2cal_fe_value_to_p_valuelog@@GLIBC_2.0Ptw_recallocpl_Doubleprint_optionsfree_HSPtrace_backpl_NBoolpl_Boolfree_scores__gmon_start__strcpy@@GLIBC_2.0