make: *** No targets. Stop. make -k T0106-blast.rdb T0106-double-blast.rdb make[1]: Entering directory `/auto/projects/compbio/experiments/casp4/t106' /projects/compbio/experiments/casp4/scripts/single-blast -q T0106.seq -interdb /projects/compbio/data/nrp/nr -db /projects/compbio/data/pdb/all-protein > T0106-blast.rdb starting T0106 with E=10 /projects/compbio/experiments/casp4/scripts/double-blast -q T0106.seq -interdb /projects/compbio/data/nrp/nr -db /projects/compbio/data/pdb/all-protein > T0106-double-blast.rdb starting T0106 with E=0.000050 T0106 T0106 blasting: T0106 gi|120529|sp|P18537|FRIZ_DROME_53:175 blasting: T0106 gi|7983|emb|CAA38459.1|_53:175 blasting: T0106 gi|103213|pir||S15709_53:175 blasting: T0106 gi|310115|gb|AAA41172.1|_44:156 blasting: T0106 gi|310117|gb|AAA41173.1|_111:222 blasting: T0106 gi|348411|pir||A45054_111:222 blasting: T0106 gi|618430|gb|AAC52179.1|_374:480 blasting: T0106 gi|841240|gb|AAA91635.1|_374:480 blasting: T0106 gi|882059|gb|AAB38383.1|_59:171 blasting: T0106 gi|1151180|gb|AAC52429.1|_27:134 blasting: T0106 gi|1151250|gb|AAC46993.1|_26:131 blasting: T0106 gi|1167906|gb|AAC52903.1|_374:480 blasting: T0106 gi|1518051|gb|AAC47273.1|_64:170 blasting: T0106 gi|1589780|gb|AAC47301.1|_28:141 blasting: T0106 gi|1620539|gb|AAC48662.1|_33:160 blasting: T0106 gi|1688095|gb|AAC51217.1|_33:160 blasting: T0106 gi|1688099|gb|AAC60113.1|_29:157 blasting: T0106 gi|1747302|dbj|BAA09053.1|_40:153 blasting: T0106 gi|1854030|gb|AAC60114.1|_29:157 blasting: T0106 gi|1917007|gb|AAB51298.1|_33:160 blasting: T0106 gi|1917009|gb|AAB51299.1|_1:98 blasting: T0106 gi|2226372|gb|AAB61752.1|_48:155 blasting: T0106 gi|7500013|pir||T31690_1:105 blasting: T0106 gi|2293564|gb|AAB65431.1|_22:144 blasting: T0106 gi|2337937|gb|AAB67062.1|_56:162 blasting: T0106 gi|2415415|gb|AAB70792.1|_40:151 blasting: T0106 gi|2415417|gb|AAB70793.1|_62:168 blasting: T0106 gi|2415421|gb|AAB70795.1|_40:153 blasting: T0106 gi|2440219|emb|CAA04977.1|_26:133 blasting: T0106 gi|7494607|pir||T42210_37:146 blasting: T0106 gi|2655274|gb|AAB87968.1|_70:186 blasting: T0106 gi|2655276|gb|AAB87969.1|_47:159 blasting: T0106 gi|2921088|gb|AAC04668.1|_44:163 blasting: T0106 gi|2921090|gb|AAC04669.1|_44:163 blasting: T0106 gi|3089551|gb|AAC15063.1|_31:137 blasting: T0106 gi|3323438|gb|AAC31121.1|_27:132 blasting: T0106 gi|3434981|dbj|BAA32424.1|_49:169 blasting: T0106 gi|7511863|pir||T13484_25:111 blasting: T0106 gi|3869266|gb|AAC77361.1|_29:134 blasting: T0106 gi|3983449|gb|AAC83940.1|_31:136 blasting: T0106 gi|4105677|gb|AAD02492.1|_58:165 blasting: T0106 gi|4164471|gb|AAD05435.1|_30:135 blasting: T0106 gi|4335927|gb|AAD17520.1|_31:136 blasting: T0106 gi|4469405|gb|AAD21247.1|_37:149 blasting: T0106 gi|4545254|gb|AAD22462.1|AF116853_1_29:156 blasting: T0106 gi|4583384|gb|AAD25051.1|AF117757_1_55:162 blasting: T0106 gi|4583388|gb|AAD25053.1|AF117759_1_54:161 blasting: T0106 gi|4585277|gb|AAD25357.1|AF117387_1_45:150 blasting: T0106 gi|4585279|gb|AAD25358.1|AF117388_1_35:139 blasting: T0106 gi|4585281|gb|AAD25359.1|AF117389_1_30:135 blasting: T0106 gi|4585283|gb|AAD25360.1|AF117390_1_31:136 blasting: T0106 gi|4689161|gb|AAD27789.1|AF088850_1_38:158 blasting: T0106 gi|5042380|gb|AAB87503.2|_2:115 blasting: T0106 gi|5262975|emb|CAB45875.1|_37:149 blasting: T0106 gi|5305407|gb|AAD41636.1|AF072872_1_117:229 blasting: T0106 gi|5524663|gb|AAD44331.1|_36:148 blasting: T0106 gi|5524665|gb|AAD44332.1|_34:138 blasting: T0106 gi|5815434|gb|AAD52671.1|_37:149 blasting: T0106 gi|5821259|dbj|BAA83742.1|_38:142 blasting: T0106 gi|5931528|dbj|BAA84677.1|_40:166 blasting: T0106 gi|7511690|pir||T37325_25:144 blasting: T0106 gi|6007775|gb|AAF01036.1|AF183910_1_46:150 blasting: T0106 gi|6273776|gb|AAF06359.1|AF139165_1_38:150 blasting: T0106 gi|6539733|gb|AAF15998.1|_1:106 blasting: T0106 gi|6572469|emb|CAB63117.1|_30:135 blasting: T0106 gi|6578963|gb|AAF18154.1|AF209200_1_27:133 blasting: T0106 gi|6636325|gb|AAF20143.1|AF207665_1_1:113 blasting: T0106 gi|6649948|gb|AAF21644.1|AF039411_1_35:139 blasting: T0106 gi|6649950|gb|AAF21645.1|AF039412_1_30:135 blasting: T0106 gi|6746598|gb|AAF27642.1|AF218056_1_37:150 blasting: T0106 gi|6746600|gb|AAF27643.1|AF218057_1_25:152 blasting: T0106 gi|6863203|gb|AAC12873.2|_111:223 blasting: T0106 gi|7106481|dbj|BAA92184.1|_126:231 blasting: T0106 gi|7110524|gb|AAF36979.1|AF231711_1_51:163 blasting: T0106 gi|7290081|gb|AAF45547.1|_40:166 blasting: T0106 gi|7290800|gb|AAF46245.1|_46:151 blasting: T0106 gi|7293818|gb|AAF49184.1|_64:170 blasting: T0106 gi|7294400|gb|AAF49746.1|_53:175 blasting: T0106 gi|7304194|gb|AAF59230.1|_761:862 blasting: T0106 gi|7331805|gb|AAF60493.1|_26:131 blasting: T0106 gi|7340122|gb|AAF61095.1|AF224315_1_6:118 blasting: T0106 gi|7340124|gb|AAF61096.1|AF224316_1_33:137 blasting: T0106 gi|7340126|gb|AAF61097.1|AF224317_1_47:159 blasting: T0106 gi|7381171|gb|AAF61420.1|_42:149 blasting: T0106 gi|7381173|gb|AAF61421.1|AF136184_1_26:136 blasting: T0106 gi|7446033|pir||B56101_374:480 blasting: T0106 gi|7451362|pir||S71786_64:170 blasting: T0106 gi|7512441|pir||JE0175_57:164 blasting: T0106 gi|7512442|pir||JE0174_40:153 blasting: T0106 gi|7527456|gb|AAF63152.1|AF039215_1_37:149 blasting: T0106 gi|7532831|gb|AAF63250.1|AF195242_1_40:166 blasting: T0106 gi|7672423|gb|AAF66480.1|AF140346_1_22:144 blasting: T0106 gi|7672425|gb|AAF66481.1|_22:144 blasting: T0106 gi|7839499|gb|AAF70299.1|AF255339_1_43:121 blasting: T0106 gi|8248134|gb|AAF74056.1|AF206321_1_44:165 blasting: make[1]: Leaving directory `/auto/projects/compbio/experiments/casp4/t106' make -k T0106.t2k.a2m.gz \ T0106.t2k.pa \ T0106.t2k.pa.html \ T0106.t2k.2d \ T0106.t2k.tree \ T0106.t2k_sorted.pa \ T0106.t2k.tree.ps \ T0106.t2k.tree-unroot.ps make[1]: Entering directory `/auto/projects/compbio/experiments/casp4/t106' /projects/compbio/experiments/models.97/scripts2k/target2k -out T0106.t2k \ -seed T0106.seq \ -all @@@@ chgrp protein /projects/compbio/tmp/target2k-eclipse-16089 @@@@ cp /projects/compbio/experiments/casp4/t106/T0106.seq m0.a2m @@@@ checkseq foo -alphabet protein -db m0.a2m > init.check SAM: checkseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:44 # About to count sequences in /projects/compbio/data/nrp/nr /projects/compbio/experiments/models.97/scripts2k/target2k: E-value thresholds used are prefilter: 0.01, actual: 0.0001 prefilter: 1, actual: 0.0002 prefilter: 10, actual: 0.001 prefilter: 400, actual: 0.005 /projects/compbio/experiments/models.97/scripts2k/target2k: aweight_bits are 0.8, 0.7, 0.6, 0.5 /projects/compbio/experiments/models.97/scripts2k/target2k: db_size= 514491 reverse_diff= 4 /projects/compbio/experiments/models.97/scripts2k/target2k: blast_max_report= 20000 /projects/compbio/experiments/models.97/scripts2k/target2k: db=/projects/compbio/data/nrp/nr /projects/compbio/experiments/models.97/scripts2k/target2k: aweight_method=1 aweight_exponent=10 /projects/compbio/experiments/models.97/scripts2k/target2k: force_seed=1 constraints=1 /projects/compbio/experiments/models.97/scripts2k/target2k: jump_in=0.2 jump_out=1 fimstrength=1 fimtrans=-1 @@@@ blast-prefilter -out prefilter -seed m0.a2m -prefilter_thresholds 0.01,1,10,400 -db /projects/compbio/data/nrp/nr -blast_max_report 20000 @@@@ chgrp protein /projects/compbio/tmp/blast-prefilter-eclipse-6710 /projects/compbio/experiments/models.97/scripts2k/blast-prefilter: E-value thresholds used are For /projects/compbio/tmp/target2k-eclipse-16089/prefilter_1.fa, 0.01 For /projects/compbio/tmp/target2k-eclipse-16089/prefilter_2.fa, 1 For /projects/compbio/tmp/target2k-eclipse-16089/prefilter_3.fa, 10 For /projects/compbio/tmp/target2k-eclipse-16089/prefilter_4.fa, 400 /projects/compbio/experiments/models.97/scripts2k/blast-prefilter: db=/projects/compbio/data/nrp/nr @@@@ uniqueseq init -alignfile /projects/compbio/tmp/target2k-eclipse-16089/m0.a2m -a protein -a2mdots 0 SAM: uniqueseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:50 Reading alignment file /projects/compbio/tmp/target2k-eclipse-16089/m0.a2m (1 sequences, 128 columns) as A2M alignment. Writing sequence output to init.a2m. No sequences have been droppd. @@@@ /projects/compbio/bin/wu-blastp /projects/compbio/data/nrp/nr 1.fasta E=400 V=20000 B=0 -gi > 1.fasta-blast.out Searched 514491 sequences # parsing wu-blastp output @@@@ rm -f 1.fasta < /dev/null @@@@ rm -f 1.fasta-blast.out < /dev/null # Extracting hits in FASTA format from /projects/compbio/data/nrp/nr @@@@ rm -rf /projects/compbio/tmp/blast-prefilter-eclipse-6710 < /dev/null @@@@ modelfromalign unused -alignfile m0.a2m \ -constraints_from_align 1 -constraints_out m0.cst SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file m0.a2m (1 sequences, 128 columns) as A2M alignment. # About to count sequences in m0.a2m # About to count sequences in prefilter_1.fa @@@@ modelfromalign tmp_1-a -alignfile m0.a2m \ -insert /projects/compbio/lib/gap1.5.regularizer \ -aweight_bits 0.8\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 -constraints m0.cst -constraints_out tmp_1-a.cst Reading parameter file /projects/compbio/lib/gap1.5.regularizer /projects/compbio/lib/gap1.5.regularizer(1): Reading REGULARIZER: gap1.5. Weak regularizer, gaps average 1.5, matches average 10 SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file m0.a2m (1 sequences, 128 columns) as A2M alignment. @@@@ hmmscore tmp_1-a -i tmp_1-a.mod -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -dpstyle 0 \ -subtract_null 4 -sort 2 \ -db_size 514491 \ -simple_threshold -26.3612737131737 -Emax 0.0001 \ -adpstyle 1 \ -align_short 0 -mdEmax 0.0001\ -select_score 4 \ -select_seq 4 \ -select_mdalign 4 -select_md 4 -Motifcutoff 0.15625 -db prefilter_1.fa -constraints tmp_1-a.cst -constraints_out tmp_1-a-md.cst Reading parameter file tmp_1-a.mod tmp_1-a.mod(23): Reading MODEL -- Model from alignment file m0.a2m SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -95.254944 Database has 124 sequences with 67183 residues. @@@@ sortseq tmp_1-a-sort -alignfile tmp_1-a.mult \ -NLLfile tmp_1-a.mstat -a2mdots 0 SAM: sortseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:38 Reading alignment file tmp_1-a.mult (122 sequences, 128 columns) as A2M alignment. Reading scores from file tmp_1-a.mstat Writing sequence output to tmp_1-a-sort.a2m. @@@@ cat m0.a2m tmp_1-a-sort.a2m > tmp_1-a.train.seq @@@@ cat m0.cst tmp_1-a-md.cst > tmp_1-a.train.cst # About to count sequences in tmp_1-a.train.seq # num_seqs=123 nll_thresh=-22.3612737131737 frac_id=0.888259788823955 @@@@ uniqueseq thinned_tmp -alignfile tmp_1-a.train.seq -a protein \ -percent_id 0.888259788823955 -a2mdots 0 SAM: uniqueseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:50 Reading alignment file tmp_1-a.train.seq (123 sequences, 128 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 34 (of 123) duplicate sequences with differing IDs Dropping 47 (of 123) sequences with > 88.8% id 42 sequences left after dropping 81 of 123 sequences. @@@@ modelfromalign junk_for_weights \ -alignfile tmp_1-a.train.seq \ -aweight_bits 0.8\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file tmp_1-a.train.seq (42 sequences, 128 columns) as A2M alignment. @@@@ buildmodel tmp_1-b -insert tmp_1-a.mod\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -modellength 0 \ -internal_weight 0 \ -anneal_length 2 \ -many_files 1 \ -Nseq 20000\ -nsurgery 0 \ -fracinsert 0.2 \ -cutinsert 0.25 -cutmatch 0.2 -initial_noise 0.0001 -anneal_noise 0.000001 \ -internal_weight 0 \ -train tmp_1-a.train.seq \ -nmodels 1\ -insert /projects/compbio/lib/gap1.5.regularizer -sequence_weights tmp_1-a.w -constraints tmp_1-a.train.cst Reading parameter file tmp_1-a.mod tmp_1-a.mod(23): Reading MODEL -- Model from alignment file m0.a2m Command line option(48): internal_weight has already been set. Using 0. Reading parameter file /projects/compbio/lib/gap1.5.regularizer /projects/compbio/lib/gap1.5.regularizer(1): Reading REGULARIZER: gap1.5. Weak regularizer, gaps average 1.5, matches average 10 SAM: buildmodel v3.1b (February 24, 1999) compiled 04/18/00_11:43:59 -376.90 -56.35 -129.20 45.03 10 0 129 @@@@ hmmscore tmp_1-b -i tmp_1-b.mod -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -dpstyle 0 \ -subtract_null 4 -sort 2 \ -db_size 514491 \ -db m0.a2m -db tmp_1-a.sel \ -select_score 4 \ -select_seq 4 \ -simple_threshold -26.3612737131737 -Emax 0.0001 \ -select_mdalign 4 -select_md 4 \ -adpstyle 1 \ -Motifcutoff 0.15625 \ -align_short 5 -mdEmax 0.0001 -constraints tmp_1-a.train.cst -constraints_out /projects/compbio/experiments/casp4/t106/T0106.t2k_1.cst Reading parameter file tmp_1-b.mod tmp_1-b.mod(20): Reading MODEL -- Final model for run tmp_1-b SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -139.799896 Database has 124 sequences with 66517 residues. @@@@ cp -f tmp_1-b.mult /projects/compbio/experiments/casp4/t106/T0106.t2k_1.a2m # About to count sequences in prefilter_2.fa @@@@ modelfromalign tmp_2-a -alignfile /projects/compbio/experiments/casp4/t106/T0106.t2k_1.a2m \ -insert /projects/compbio/lib/stiff-gap5.regularizer \ -aweight_bits 0.7\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 -constraints /projects/compbio/experiments/casp4/t106/T0106.t2k_1.cst -constraints_out tmp_2-a.cst Reading parameter file /projects/compbio/lib/stiff-gap5.regularizer /projects/compbio/lib/stiff-gap5.regularizer(1): Reading REGULARIZER: gap5. stiff regularizer, gaps averge 5, matches average 10 SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file /projects/compbio/experiments/casp4/t106/T0106.t2k_1.a2m (126 sequences, 128 columns) as A2M alignment. @@@@ hmmscore tmp_2-a -i tmp_2-a.mod -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -dpstyle 0 \ -subtract_null 4 -sort 2 \ -db_size 514491 \ -simple_threshold -25.6681265324194 -Emax 0.0002 \ -adpstyle 1 \ -align_short 0 -mdEmax 0.0002\ -select_score 4 \ -select_seq 4 \ -select_mdalign 4 -select_md 4 -Motifcutoff 0.15625 -db prefilter_2.fa -constraints tmp_2-a.cst -constraints_out tmp_2-a-md.cst Reading parameter file tmp_2-a.mod tmp_2-a.mod(23): Reading MODEL -- Model from alignment file /projects/compbio/experiments/casp4/t106/T0106.t2k_1.a2m SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -129.729782 Database has 126 sequences with 68128 residues. @@@@ sortseq tmp_2-a-sort -alignfile tmp_2-a.mult \ -NLLfile tmp_2-a.mstat -a2mdots 0 SAM: sortseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:38 Reading alignment file tmp_2-a.mult (128 sequences, 128 columns) as A2M alignment. Reading scores from file tmp_2-a.mstat Writing sequence output to tmp_2-a-sort.a2m. @@@@ cat m0.a2m tmp_2-a-sort.a2m > tmp_2-a.train.seq @@@@ cat m0.cst tmp_2-a-md.cst > tmp_2-a.train.cst # About to count sequences in tmp_2-a.train.seq # num_seqs=129 nll_thresh=-21.6681265324194 frac_id=0.874291813185283 @@@@ uniqueseq thinned_tmp -alignfile tmp_2-a.train.seq -a protein \ -percent_id 0.874291813185283 -a2mdots 0 SAM: uniqueseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:50 Reading alignment file tmp_2-a.train.seq (129 sequences, 128 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 33 (of 129) duplicate sequences with differing IDs Dropping 53 (of 129) sequences with > 87.4% id 43 sequences left after dropping 86 of 129 sequences. @@@@ modelfromalign junk_for_weights \ -alignfile tmp_2-a.train.seq \ -aweight_bits 0.7\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file tmp_2-a.train.seq (43 sequences, 128 columns) as A2M alignment. @@@@ buildmodel tmp_2-b -insert tmp_2-a.mod\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -modellength 0 \ -internal_weight 0 \ -anneal_length 2 \ -many_files 1 \ -Nseq 20000\ -nsurgery 0 \ -fracinsert 0.2 \ -cutinsert 0.25 -cutmatch 0.2 -initial_noise 0.0001 -anneal_noise 0.000001 \ -internal_weight 0 \ -train tmp_2-a.train.seq \ -nmodels 1\ -insert /projects/compbio/lib/stiff-gap5.regularizer -sequence_weights tmp_2-a.w -constraints tmp_2-a.train.cst Reading parameter file tmp_2-a.mod tmp_2-a.mod(23): Reading MODEL -- Model from alignment file /projects/compbio/experiments/casp4/t106/T0106.t2k_1.a2m Command line option(48): internal_weight has already been set. Using 0. Reading parameter file /projects/compbio/lib/stiff-gap5.regularizer /projects/compbio/lib/stiff-gap5.regularizer(1): Reading REGULARIZER: gap5. stiff regularizer, gaps averge 5, matches average 10 SAM: buildmodel v3.1b (February 24, 1999) compiled 04/18/00_11:43:59 -375.77 -64.42 -115.41 45.76 5 0 129 @@@@ hmmscore tmp_2-b -i tmp_2-b.mod -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -dpstyle 0 \ -subtract_null 4 -sort 2 \ -db_size 514491 \ -db m0.a2m -db tmp_2-a.sel \ -select_score 4 \ -select_seq 4 \ -simple_threshold -25.6681265324194 -Emax 0.0002 \ -select_mdalign 4 -select_md 4 \ -adpstyle 1 \ -Motifcutoff 0.15625 \ -align_short 5 -mdEmax 0.0002 -constraints tmp_2-a.train.cst -constraints_out /projects/compbio/experiments/casp4/t106/T0106.t2k_2.cst Reading parameter file tmp_2-b.mod tmp_2-b.mod(20): Reading MODEL -- Final model for run tmp_2-b SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -123.868576 Database has 127 sequences with 68256 residues. @@@@ cp -f tmp_2-b.mult /projects/compbio/experiments/casp4/t106/T0106.t2k_2.a2m # About to count sequences in prefilter_3.fa @@@@ modelfromalign tmp_3-a -alignfile /projects/compbio/experiments/casp4/t106/T0106.t2k_2.a2m \ -insert /projects/compbio/lib/stiff-gap5.regularizer \ -aweight_bits 0.6\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 -constraints /projects/compbio/experiments/casp4/t106/T0106.t2k_2.cst -constraints_out tmp_3-a.cst Reading parameter file /projects/compbio/lib/stiff-gap5.regularizer /projects/compbio/lib/stiff-gap5.regularizer(1): Reading REGULARIZER: gap5. stiff regularizer, gaps averge 5, matches average 10 SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file /projects/compbio/experiments/casp4/t106/T0106.t2k_2.a2m (129 sequences, 128 columns) as A2M alignment. @@@@ hmmscore tmp_3-a -i tmp_3-a.mod -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -dpstyle 0 \ -subtract_null 4 -sort 2 \ -db_size 514491 \ -simple_threshold -24.0586886184303 -Emax 0.001 \ -adpstyle 1 \ -align_short 0 -mdEmax 0.001\ -select_score 4 \ -select_seq 4 \ -select_mdalign 4 -select_md 4 -Motifcutoff 0.15625 -db prefilter_3.fa -constraints tmp_3-a.cst -constraints_out tmp_3-a-md.cst Reading parameter file tmp_3-a.mod tmp_3-a.mod(23): Reading MODEL -- Model from alignment file /projects/compbio/experiments/casp4/t106/T0106.t2k_2.a2m SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -120.640198 Database has 128 sequences with 68318 residues. @@@@ sortseq tmp_3-a-sort -alignfile tmp_3-a.mult \ -NLLfile tmp_3-a.mstat -a2mdots 0 SAM: sortseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:38 Reading alignment file tmp_3-a.mult (128 sequences, 128 columns) as A2M alignment. Reading scores from file tmp_3-a.mstat Writing sequence output to tmp_3-a-sort.a2m. @@@@ cat m0.a2m tmp_3-a-sort.a2m > tmp_3-a.train.seq @@@@ cat m0.cst tmp_3-a-md.cst > tmp_3-a.train.cst # About to count sequences in tmp_3-a.train.seq # num_seqs=129 nll_thresh=-20.0586886184303 frac_id=0.860500527442034 @@@@ uniqueseq thinned_tmp -alignfile tmp_3-a.train.seq -a protein \ -percent_id 0.860500527442034 -a2mdots 0 SAM: uniqueseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:50 Reading alignment file tmp_3-a.train.seq (129 sequences, 128 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 30 (of 129) duplicate sequences with differing IDs Dropping 55 (of 129) sequences with > 86.1% id 44 sequences left after dropping 85 of 129 sequences. @@@@ modelfromalign junk_for_weights \ -alignfile tmp_3-a.train.seq \ -aweight_bits 0.6\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file tmp_3-a.train.seq (44 sequences, 128 columns) as A2M alignment. @@@@ buildmodel tmp_3-b -insert tmp_3-a.mod\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -modellength 0 \ -internal_weight 0 \ -anneal_length 2 \ -many_files 1 \ -Nseq 20000\ -nsurgery 0 \ -fracinsert 0.2 \ -cutinsert 0.25 -cutmatch 0.2 -initial_noise 0.0001 -anneal_noise 0.000001 \ -internal_weight 0 \ -train tmp_3-a.train.seq \ -nmodels 1\ -insert /projects/compbio/lib/stiff-gap5.regularizer -sequence_weights tmp_3-a.w -constraints tmp_3-a.train.cst Reading parameter file tmp_3-a.mod tmp_3-a.mod(23): Reading MODEL -- Model from alignment file /projects/compbio/experiments/casp4/t106/T0106.t2k_2.a2m Command line option(48): internal_weight has already been set. Using 0. Reading parameter file /projects/compbio/lib/stiff-gap5.regularizer /projects/compbio/lib/stiff-gap5.regularizer(1): Reading REGULARIZER: gap5. stiff regularizer, gaps averge 5, matches average 10 SAM: buildmodel v3.1b (February 24, 1999) compiled 04/18/00_11:43:59 -376.03 -59.30 -110.21 45.94 6 0 129 @@@@ hmmscore tmp_3-b -i tmp_3-b.mod -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -dpstyle 0 \ -subtract_null 4 -sort 2 \ -db_size 514491 \ -db m0.a2m -db tmp_3-a.sel \ -select_score 4 \ -select_seq 4 \ -simple_threshold -24.0586886184303 -Emax 0.001 \ -select_mdalign 4 -select_md 4 \ -adpstyle 1 \ -Motifcutoff 0.15625 \ -align_short 5 -mdEmax 0.001 -constraints tmp_3-a.train.cst -constraints_out /projects/compbio/experiments/casp4/t106/T0106.t2k_3.cst Reading parameter file tmp_3-b.mod tmp_3-b.mod(20): Reading MODEL -- Final model for run tmp_3-b SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -118.750031 Database has 127 sequences with 68256 residues. @@@@ cp -f tmp_3-b.mult /projects/compbio/experiments/casp4/t106/T0106.t2k_3.a2m # About to count sequences in prefilter_4.fa @@@@ modelfromalign tmp_4-a -alignfile /projects/compbio/experiments/casp4/t106/T0106.t2k_3.a2m \ -insert /projects/compbio/lib/fssp-trained.regularizer \ -aweight_bits 0.5\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 -constraints /projects/compbio/experiments/casp4/t106/T0106.t2k_3.cst -constraints_out tmp_4-a.cst Reading parameter file /projects/compbio/lib/fssp-trained.regularizer /projects/compbio/lib/fssp-trained.regularizer(1): Reading REGULARIZER: transitions from FSSP and pseudocounts for diverse weighted HSSP SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file /projects/compbio/experiments/casp4/t106/T0106.t2k_3.a2m (129 sequences, 128 columns) as A2M alignment. @@@@ hmmscore tmp_4-a -i tmp_4-a.mod -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -dpstyle 0 \ -subtract_null 4 -sort 2 \ -db_size 514491 \ -simple_threshold -22.4492506982216 -Emax 0.005 \ -adpstyle 1 \ -align_short 0 -mdEmax 0.005\ -select_score 4 \ -select_seq 4 \ -select_mdalign 4 -select_md 4 -Motifcutoff 0.15625 -db prefilter_4.fa -constraints tmp_4-a.cst -constraints_out tmp_4-a-md.cst Reading parameter file tmp_4-a.mod tmp_4-a.mod(23): Reading MODEL -- Model from alignment file /projects/compbio/experiments/casp4/t106/T0106.t2k_3.a2m SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -51.503258 Database has 322 sequences with 95484 residues. @@@@ sortseq tmp_4-a-sort -alignfile tmp_4-a.mult \ -NLLfile tmp_4-a.mstat -a2mdots 0 SAM: sortseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:38 Reading alignment file tmp_4-a.mult (128 sequences, 128 columns) as A2M alignment. Reading scores from file tmp_4-a.mstat Writing sequence output to tmp_4-a-sort.a2m. @@@@ cat m0.a2m tmp_4-a-sort.a2m > tmp_4-a.train.seq @@@@ cat m0.cst tmp_4-a-md.cst > tmp_4-a.train.cst # About to count sequences in tmp_4-a.train.seq # num_seqs=129 nll_thresh=-18.4492506982216 frac_id=0.845542799346682 @@@@ uniqueseq thinned_tmp -alignfile tmp_4-a.train.seq -a protein \ -percent_id 0.845542799346682 -a2mdots 0 SAM: uniqueseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:50 Reading alignment file tmp_4-a.train.seq (129 sequences, 128 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 30 (of 129) duplicate sequences with differing IDs Dropping 57 (of 129) sequences with > 84.6% id 42 sequences left after dropping 87 of 129 sequences. @@@@ modelfromalign junk_for_weights \ -alignfile tmp_4-a.train.seq \ -aweight_bits 0.5\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file tmp_4-a.train.seq (42 sequences, 128 columns) as A2M alignment. @@@@ buildmodel tmp_4-b -insert tmp_4-a.mod\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -modellength 0 \ -internal_weight 0 \ -anneal_length 2 \ -many_files 1 \ -Nseq 20000\ -nsurgery 0 \ -fracinsert 0.2 \ -cutinsert 0.25 -cutmatch 0.2 -initial_noise 0.0001 -anneal_noise 0.000001 \ -internal_weight 0 \ -train tmp_4-a.train.seq \ -nmodels 1\ -insert /projects/compbio/lib/fssp-trained.regularizer -sequence_weights tmp_4-a.w -constraints tmp_4-a.train.cst Reading parameter file tmp_4-a.mod tmp_4-a.mod(23): Reading MODEL -- Model from alignment file /projects/compbio/experiments/casp4/t106/T0106.t2k_3.a2m Command line option(48): internal_weight has already been set. Using 0. Reading parameter file /projects/compbio/lib/fssp-trained.regularizer /projects/compbio/lib/fssp-trained.regularizer(1): Reading REGULARIZER: transitions from FSSP and pseudocounts for diverse weighted HSSP SAM: buildmodel v3.1b (February 24, 1999) compiled 04/18/00_11:43:59 -376.88 -54.51 -105.51 47.56 5 0 129 @@@@ hmmscore tmp_4-b -i tmp_4-b.mod -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -dpstyle 0 \ -subtract_null 4 -sort 2 \ -db_size 514491 \ -db m0.a2m -db tmp_4-a.sel \ -select_score 4 \ -select_seq 4 \ -simple_threshold -22.4492506982216 -Emax 0.005 \ -select_mdalign 4 -select_md 4 \ -adpstyle 1 \ -Motifcutoff 0.15625 \ -align_short 5 -mdEmax 0.005 -constraints tmp_4-a.train.cst -constraints_out /projects/compbio/experiments/casp4/t106/T0106.t2k_4.cst Reading parameter file tmp_4-b.mod tmp_4-b.mod(20): Reading MODEL -- Final model for run tmp_4-b SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -112.144600 Database has 127 sequences with 68256 residues. @@@@ cp -f tmp_4-b.mult /projects/compbio/experiments/casp4/t106/T0106.t2k_4.a2m @@@@ modelfromalign tmp_5-a -alignfile /projects/compbio/experiments/casp4/t106/T0106.t2k_4.a2m \ -insert /projects/compbio/lib/fssp-trained.regularizer \ -aweight_bits 0.5\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 -constraints /projects/compbio/experiments/casp4/t106/T0106.t2k_4.cst -constraints_out tmp_5-a.cst Reading parameter file /projects/compbio/lib/fssp-trained.regularizer /projects/compbio/lib/fssp-trained.regularizer(1): Reading REGULARIZER: transitions from FSSP and pseudocounts for diverse weighted HSSP SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file /projects/compbio/experiments/casp4/t106/T0106.t2k_4.a2m (129 sequences, 128 columns) as A2M alignment. @@@@ hmmscore tmp_5-a -i tmp_5-a.mod -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -dpstyle 0 \ -subtract_null 4 -sort 2 \ -db_size 514491 \ -simple_threshold -22.4492506982216 -Emax 0.005 \ -adpstyle 5 \ -align_short 0 -mdEmax 0.005\ -select_score 4 \ -select_seq 4 \ -select_mdalign 4 -select_md 4 -Motifcutoff 0.15625 -db prefilter_4.fa -constraints tmp_5-a.cst -constraints_out tmp_5-a-md.cst Reading parameter file tmp_5-a.mod tmp_5-a.mod(23): Reading MODEL -- Model from alignment file /projects/compbio/experiments/casp4/t106/T0106.t2k_4.a2m SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -51.483276 Database has 322 sequences with 95484 residues. @@@@ sortseq tmp_5-a-sort -alignfile tmp_5-a.mult \ -NLLfile tmp_5-a.mstat -a2mdots 0 SAM: sortseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:38 Reading alignment file tmp_5-a.mult (127 sequences, 128 columns) as A2M alignment. Reading scores from file tmp_5-a.mstat Writing sequence output to tmp_5-a-sort.a2m. @@@@ cat m0.a2m tmp_5-a-sort.a2m > tmp_5-a.train.seq @@@@ cat m0.cst tmp_5-a-md.cst > tmp_5-a.train.cst @@@@ cp tmp_5-a.train.seq /projects/compbio/experiments/casp4/t106/T0106.t2k.a2m @@@@ cp tmp_5-a.train.cst /projects/compbio/experiments/casp4/t106/T0106.t2k.cst @@@@ rm -rf /projects/compbio/tmp/target2k-eclipse-16089 < /dev/null gzip -f T0106.t2k.a2m T0106.t2k_*.a2m prettyalign T0106.t2k.a2m.gz -m8 -i -n -L3333 > T0106.t2k.pa SAM: prettyalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:15 gunzip -c T0106.t2k.a2m.gz > tmp.a2m a2m2html -a2m_in tmp.a2m -retrieve > T0106.t2k.pa.html SAM: /projects/compbio/bin/alpha/prettyalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:15 rm tmp.a2m echo ReadNeuralNet /projects/compbio/usr/karplus/predict-2nd/testing/networks/overrep-2500-IDaa13-5-10-7-10-5-10-11-ehl2-seeded9-stride-trained.net > tmp.script echo ReadA2M T0106.t2k.a2m.gz >> tmp.script echo PrintPrediction T0106.t2k.2d 3670-4530-6947 >> tmp.script echo PrintPredictionFasta T0106.t2k.2d.seq >> tmp.script echo PrintRDB T0106.t2k.2d.rdb >> tmp.script predict-2nd < tmp.script # command:# Neural network set to overrep-2500-IDaa13-5-10-7-10-5-10-11-ehl2-seeded9-stride-trained.net # command:# Reading A2M format from T0106.t2k.a2m.gz # Using SequenceWeight EntropyWeight(1.4, 10) # Using regularizer /projects/compbio/lib/recode3.20comp for sequence weight regularizer. # T0106.t2k.a2m with 128 sequences, total weight= 35.5883 avg weight= 0.278034 46 iterations # AdjustWeights couldn't save exactly 1.4 bits/position, saving 1.39913 bits. Regularizing alignment for T0106.t2k.a2m.gz # After reading T0106.t2k.a2m.gz, have 128 columns in 1 chains # command:# Initializing Gain for FinalLayer # Initializing Gain for ordering_layer # Initializing Gain for SecondLayer # Initializing Gain for FirstLayer # Network initialization done # Printing prediction to T0106.t2k.2d # command:# Network initialization done # Printing prediction in FASTA format to T0106.t2k.2d.seq # command:# Network initialization done # Printing prediction to T0106.t2k.2d.rdb # command:rm tmp.script gunzip -f T0106.t2k.a2m.gz /projects/compbio/usr/karplus/src/phytree/phytree -o -i -r flat T0106.t2k T0106.t2k.a2m /projects/compbio/lib/recode2.20comp Prior library /projects/compbio/lib/recode2.20comp read. Reading alignment file T0106.t2k.a2m (128 sequences) as A2M alignment. Alignment T0106.t2k.a2m read. Adding sequences:0 10 20 30 40 50 60 70 80 90 100 110 120 One dot for every 10 merges: ............. Merges done. CPU = user 0: 0:43 system 0: 0: 0 gzip -f T0106.t2k.a2m gzip -f T0106.t2k_sorted.a2m rm T0106.t2k.phytrace T0106.t2k.tree_weight prettyalign T0106.t2k_sorted.a2m.gz -m8 -i -n -L3333 > T0106.t2k_sorted.pa SAM: prettyalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:15 /projects/compbio/usr/karplus/src/phytree/dg T0106.t2k.tree /projects/compbio/usr/karplus/src/phytree/dtree T0106.t2k.tree make[1]: Leaving directory `/auto/projects/compbio/experiments/casp4/t106' /projects/compbio/experiments/models.97/scripts99/score-target -alignment T0106.t2k.a2m.gz -build w0.5 \ -db /projects/compbio/data/pdb/all-protein \ -Emax 40.0 \ -viterbi 0 -sw_score 2 \ -output T0106.t2k-w0.5-pdb.rdb PATH=.:/projects/compbio/experiments/models.97/scripts99:/projects/compbio/bin:/projects/compbio/bin/alpha:/projects/compbio/bin/hmmscripts:/projects/compbio/bin/scripts:/sbin:/usr/sbin:/bin:/usr/bin:/usr/local/gnu/bin:/usr/local/share/request/bin:/usr/local/bin @@@@ mkdir /projects/compbio/tmp/score-target-15681 @@@@ mkdir /projects/compbio/tmp/score-target-15681/for-scorescript @@@@ /projects/compbio/experiments/models.97/scripts99/score-mods-vs-seqs.pl -root_tmp /projects/compbio/tmp/score-target-15681/for-scorescript -a2msfile /projects/compbio/tmp/score-target-15681/score-target.a2m \ -buildscript /projects/compbio/bin/scripts/w0.5 -seqfile /projects/compbio/data/pdb/all-protein -Emax 40.0 -viterbi 0 \ -sw_score 2 -output /projects/compbio/tmp/score-target-15681/score-target.op @@@@ chgrp protein /projects/compbio/tmp/score-target-15681/for-scorescript/split-tmp @@@@ chgrp protein /projects/compbio/tmp/score-target-15681/for-scorescript/split-tmp @@@@ /projects/compbio/bin/scripts/w0.5 /projects/compbio/experiments/casp4/t106/T0106.t2k.a2m.gz /projects/compbio/tmp/score-target-15681/for-scorescript/split-tmp/eclipse-2692/T0106.t2k.mod Reading /projects/compbio/bin/scripts/sam-t99.conf @@@@ chgrp protein /projects/compbio/tmp/tmp-build-weighted-model-eclipse-22261 @@@@ uniqueseq /projects/compbio/tmp/tmp-build-weighted-model-eclipse-22261/tmp -alignfile /projects/compbio/experiments/casp4/t106/T0106.t2k.a2m.gz \ -a protein -percent_id 0.8 SAM: uniqueseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:50 Reading alignment file /projects/compbio/experiments/casp4/t106/T0106.t2k.a2m.gz (128 sequences, 128 columns) as A2M alignment. Writing sequence output to /projects/compbio/tmp/tmp-build-weighted-model-eclipse-22261/tmp.a2m. Dropping 32 (of 128) duplicate sequences with differing IDs Dropping 61 (of 128) sequences with > 80.0% id 35 sequences left after dropping 93 of 128 sequences. @@@@ modelfromalign /projects/compbio/tmp/tmp-build-weighted-model-eclipse-22261/tmp -alignfile /projects/compbio/tmp/tmp-build-weighted-model-eclipse-22261/tmp.a2m \ -insert /projects/compbio/lib/fssp-trained.regularizer \ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -prior_library /projects/compbio/lib/recode3.20comp \ -binary_output 1 -a protein -aweight_method 1 -aweight_bits 0.5 -aweight_exponent 10 Reading parameter file /projects/compbio/lib/fssp-trained.regularizer /projects/compbio/lib/fssp-trained.regularizer(1): Reading REGULARIZER: transitions from FSSP and pseudocounts for diverse weighted HSSP SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file /projects/compbio/tmp/tmp-build-weighted-model-eclipse-22261/tmp.a2m (35 sequences, 128 columns) as A2M alignment. @@@@ cp -f /projects/compbio/tmp/tmp-build-weighted-model-eclipse-22261/tmp.mod /projects/compbio/tmp/score-target-15681/for-scorescript/split-tmp/eclipse-2692/T0106.t2k.mod @@@@ rm -rf /projects/compbio/tmp/tmp-build-weighted-model-eclipse-22261