make -k T0123-blast.rdb T0123-double-blast.rdb make[1]: Entering directory `/export/projects/compbio/experiments/casp4/t123' /projects/compbio/experiments/casp4/scripts/single-blast -q T0123.seq -interdb /projects/compbio/data/nrp/nr -db /projects/compbio/data/pdb/all-protein > T0123-blast.rdb starting T0123 with E=10 /projects/compbio/experiments/casp4/scripts/double-blast -q T0123.seq -interdb /projects/compbio/data/nrp/nr -db /projects/compbio/data/pdb/all-protein > T0123-double-blast.rdb starting T0123 with E=0.000050 T0123 T0123 blasting: T0123 gi|520|emb|CAA32835.1|_19:178 blasting: T0123 gi|125912|sp|P02756|LACB_CAPHI_21:180 blasting: T0123 gi|1146|emb|CAA33282.1|_17:128 blasting: T0123 gi|125918|sp|P02757|LACB_SHEEP_21:180 blasting: T0123 gi|111366|pir||S05440_29:153 blasting: T0123 gi|137823|sp|P20289|VEG1_RAT_24:172 blasting: T0123 gi|266472|sp|Q01584|LIPO_BUFMA_23:161 blasting: T0123 gi|544253|sp|P35578|ESP4_LACVV_24:167 blasting: T0123 gi|127533|sp|P02761|MUP_RAT_21:153 blasting: T0123 gi|125911|sp|P02755|LACB_BUBAR_3:162 blasting: T0123 gi|72082|pir||LGHO_7:162 blasting: T0123 gi|125913|sp|P13613|LACB_EQUAS_7:162 blasting: T0123 gi|72084|pir||LGHO2_2:164 blasting: T0123 gi|125916|sp|P11944|LACB_MACGI_9:155 blasting: T0123 gi|89098|pir||A05144_1:77 blasting: T0123 gi|89098|pir||A05144_31:131 blasting: T0123 gi|89159|pir||A05309_1:99 blasting: T0123 gi|89159|pir||A05309_92:129 blasting: T0123 gi|89160|pir||A25804_1:159 blasting: T0123 gi|89908|pir||S00132_3:162 blasting: T0123 gi|127532|sp|P04939|MUPM_MOUSE_24:156 blasting: T0123 gi|127529|sp|P11591|MUP5_MOUSE_29:152 blasting: T0123 gi|130701|sp|P09466|PAEP_HUMAN_19:179 blasting: T0123 gi|125905|sp|P21664|LACA_FELCA_5:163 blasting: T0123 gi|125904|sp|P19647|LACA_EQUAS_2:161 blasting: T0123 gi|125910|sp|P02754|LACB_BOVIN_19:178 blasting: T0123 gi|162748|gb|AAA30411.1|_2:151 blasting: T0123 gi|165839|gb|AAA31510.1|_21:180 blasting: T0123 gi|182093|gb|AAA35802.1|_28:157 blasting: T0123 gi|2119647|pir||I52504_21:153 blasting: T0123 gi|2119648|pir||I65319_21:153 blasting: T0123 gi|204261|gb|AAA41198.1|_17:149 blasting: T0123 gi|2119649|pir||I52212_7:131 blasting: T0123 gi|223780|prf||0911271A_3:162 blasting: T0123 gi|229460|prf||732164A_3:162 blasting: T0123 gi|3122325|sp|Q29614|LACB_MACEU_27:174 blasting: T0123 gi|346460|pir||A45542_19:178 blasting: T0123 gi|462472|sp|P33685|LACA_CANFA_2:161 blasting: T0123 gi|462475|sp|P33688|LACC_FELCA_7:162 blasting: T0123 gi|462473|sp|P33687|LACB_FELCA_7:162 blasting: T0123 gi|448344|prf||1916447B_5:163 blasting: T0123 gi|462474|sp|P33686|LACC_CANFA_7:161 blasting: T0123 gi|448393|prf||1917182A_20:163 blasting: T0123 gi|731102|sp|P41244|VEG2_RAT_24:172 blasting: T0123 gi|530873|gb|AAA30412.1|_1:57 blasting: T0123 gi|1079295|pir||S52354_23:163 blasting: T0123 gi|1082934|pir||S33878_7:161 blasting: T0123 gi|2497698|sp|Q29487|PGHD_FELCA_38:169 blasting: T0123 gi|6226255|sp|Q29147|TRIC_TRIVU_36:159 blasting: T0123 gi|1196815|gb|AAA88508.1|_29:153 blasting: T0123 gi|1213589|dbj|BAA12075.1|_23:163 blasting: T0123 gi|3913984|sp|Q29146|LACB_TRIVU_27:174 blasting: T0123 gi|1279615|emb|CAA30059.1|_21:141 blasting: T0123 gi|1346407|sp|P04119|LACB_PIG_19:178 blasting: T0123 gi|1403310|emb|CAA67099.1|_29:158 blasting: T0123 gi|2851437|sp|P02758|LACB_HORSE_25:180 blasting: T0123 gi|2851599|sp|P07380|LACA_HORSE_20:179 blasting: T0123 gi|3121758|sp|Q95182|ALL1_HORSE_29:153 blasting: T0123 gi|2119651|pir||S33875_7:161 blasting: T0123 gi|2194088|pdb|1BEB|A_3:162 blasting: T0123 gi|3121746|sp|O18874|ALL2_CANFA_31:176 blasting: T0123 gi|3483096|gb|AAC33251.1|_25:186 blasting: T0123 gi|3694896|gb|AAC62487.1|_18:162 blasting: T0123 gi|3694898|gb|AAC62488.1|_14:159 blasting: T0123 gi|3703061|gb|AAC63024.1|_1:46 blasting: T0123 gi|3808097|emb|CAA06532.1|_21:180 blasting: T0123 gi|4322130|gb|AAD15971.1|_3:154 blasting: T0123 gi|4322131|gb|AAD15972.1|_3:154 blasting: T0123 gi|4322132|gb|AAD15973.1|_2:156 blasting: T0123 gi|4322133|gb|AAD15974.1|_2:156 blasting: T0123 gi|4322134|gb|AAD15975.1|_3:155 blasting: T0123 gi|6729725|pdb|1BSQ|A_3:162 blasting: T0123 gi|6980895|pdb|1BSO|A_3:162 blasting: T0123 gi|4884164|emb|CAB43305.1|_27:187 blasting: T0123 gi|7766897|pdb|1CJ5|A_3:162 blasting: T0123 gi|7766863|pdb|1EW3|A_7:131 blasting: T0123 gi|8250627|emb|CAB93679.1|_25:153 blasting: T0123 gi|8307686|dbj|BAA96480.1|_29:153 blasting: T0123 gi|8307688|dbj|BAA96481.1|_30:153 blasting: T0123 gi|8307692|dbj|BAA96483.1|_29:153 blasting: T0123 gi|8307694|dbj|BAA96484.1|_29:153 blasting: T0123 gi|8307698|dbj|BAA96486.1|_21:154 blasting: make[1]: Leaving directory `/export/projects/compbio/experiments/casp4/t123' make -k T0123.t2k.a2m.gz \ T0123.t2k.pa \ T0123.t2k.pa.html \ T0123.t2k.2d \ T0123.t2k.tree \ T0123.t2k_sorted.pa \ T0123.t2k.tree.ps \ T0123.t2k.tree-unroot.ps make[1]: Entering directory `/export/projects/compbio/experiments/casp4/t123' /projects/compbio/experiments/models.97/scripts2k/target2k -out T0123.t2k \ -seed T0123.seq \ -all @@@@ chgrp protein /projects/compbio/tmp/target2k-beta-5581 @@@@ cp /export/projects/compbio/experiments/casp4/t123/T0123.seq m0.a2m @@@@ checkseq foo -alphabet protein -db m0.a2m > init.check SAM: checkseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:52 # 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= 541007 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-beta-6012 /projects/compbio/experiments/models.97/scripts2k/blast-prefilter: E-value thresholds used are For /export/projects/compbio/tmp/target2k-beta-5581/prefilter_1.fa, 0.01 For /export/projects/compbio/tmp/target2k-beta-5581/prefilter_2.fa, 1 For /export/projects/compbio/tmp/target2k-beta-5581/prefilter_3.fa, 10 For /export/projects/compbio/tmp/target2k-beta-5581/prefilter_4.fa, 400 /projects/compbio/experiments/models.97/scripts2k/blast-prefilter: db=/projects/compbio/data/nrp/nr @@@@ uniqueseq init -alignfile /export/projects/compbio/tmp/target2k-beta-5581/m0.a2m -a protein -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file /export/projects/compbio/tmp/target2k-beta-5581/m0.a2m (1 sequences, 160 columns) as A2M alignment. Writing sequence output to init.a2m. No sequences have been dropped. @@@@ /projects/compbio/bin/wu-blastp /projects/compbio/data/nrp/nr 1.fasta E=400 V=20000 B=0 -gi > 1.fasta-blast.out Searched 541007 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-beta-6012 < /dev/null @@@@ modelfromalign unused -alignfile m0.a2m \ -constraints_from_align 1 -constraints_out m0.cst SAM: modelfromalign v3.2 (July 31, 2000) compiled 07/31/00_14:08:27 Reading alignment file m0.a2m (1 sequences, 160 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.2 (July 31, 2000) compiled 07/31/00_14:08:27 Reading alignment file m0.a2m (1 sequences, 160 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 541007 \ -simple_threshold -26.4115278685381 -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.125 -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(21): Reading MODEL -- Model from alignment file m0.a2m SAM: hmmscore v3.2 (July 31, 2000) compiled 08/04/00_11:31:33 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -83.640457 Database has 105 sequences with 17230 residues. @@@@ sortseq tmp_1-a-sort -alignfile tmp_1-a.mult \ -NLLfile tmp_1-a.mstat -a2mdots 0 SAM: sortseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:48 Reading alignment file tmp_1-a.mult (70 sequences, 160 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=71 nll_thresh=-22.4115278685381 frac_id=0.967332929461138 @@@@ uniqueseq thinned_tmp -alignfile tmp_1-a.train.seq -a protein \ -percent_id 0.967332929461138 -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file tmp_1-a.train.seq (71 sequences, 160 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 8 (of 71) duplicate sequences with differing IDs Dropping 32 (of 71) sequences with > 96.7% id 31 sequences left after dropping 40 of 71 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.2 (July 31, 2000) compiled 07/31/00_14:08:27 Reading alignment file tmp_1-a.train.seq (31 sequences, 160 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(21): 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.2 (July 31, 2000) compiled 07/31/00_14:08:12 -450.19 -18.17 -158.39 71.88 12 0 161 @@@@ 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 541007 \ -db m0.a2m -db tmp_1-a.sel \ -select_score 4 \ -select_seq 4 \ -simple_threshold -26.4115278685381 -Emax 0.0001 \ -select_mdalign 4 -select_md 4 \ -adpstyle 1 \ -Motifcutoff 0.125 \ -align_short 5 -mdEmax 0.0001 -constraints tmp_1-a.train.cst -constraints_out /export/projects/compbio/experiments/casp4/t123/T0123.t2k_1.cst Reading parameter file tmp_1-b.mod tmp_1-b.mod(18): Reading MODEL -- Final model for run tmp_1-b SAM: hmmscore v3.2 (July 31, 2000) compiled 08/04/00_11:31:33 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -159.455368 Database has 86 sequences with 13985 residues. @@@@ cp -f tmp_1-b.mult /export/projects/compbio/experiments/casp4/t123/T0123.t2k_1.a2m # About to count sequences in prefilter_2.fa @@@@ modelfromalign tmp_2-a -alignfile /export/projects/compbio/experiments/casp4/t123/T0123.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 /export/projects/compbio/experiments/casp4/t123/T0123.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.2 (July 31, 2000) compiled 07/31/00_14:08:27 Reading alignment file /export/projects/compbio/experiments/casp4/t123/T0123.t2k_1.a2m (86 sequences, 160 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 541007 \ -simple_threshold -25.7183806877933 -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.125 -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(21): Reading MODEL -- Model from alignment file /export/projects/compbio/experiments/casp4/t123/T0123.t2k_1.a2m SAM: hmmscore v3.2 (July 31, 2000) compiled 08/04/00_11:31:33 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -123.011597 Database has 147 sequences with 23975 residues. @@@@ sortseq tmp_2-a-sort -alignfile tmp_2-a.mult \ -NLLfile tmp_2-a.mstat -a2mdots 0 SAM: sortseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:48 Reading alignment file tmp_2-a.mult (138 sequences, 160 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=139 nll_thresh=-21.7183806877933 frac_id=0.861366390197216 @@@@ uniqueseq thinned_tmp -alignfile tmp_2-a.train.seq -a protein \ -percent_id 0.861366390197216 -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file tmp_2-a.train.seq (139 sequences, 160 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 15 (of 139) duplicate sequences with differing IDs Dropping 73 (of 139) sequences with > 86.1% id 51 sequences left after dropping 88 of 139 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.2 (July 31, 2000) compiled 07/31/00_14:08:27 Reading alignment file tmp_2-a.train.seq (51 sequences, 160 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(21): Reading MODEL -- Model from alignment file /export/projects/compbio/experiments/casp4/t123/T0123.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.2 (July 31, 2000) compiled 07/31/00_14:08:12 -450.37 2.37 -126.43 55.63 6 0 161 @@@@ 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 541007 \ -db m0.a2m -db tmp_2-a.sel \ -select_score 4 \ -select_seq 4 \ -simple_threshold -25.7183806877933 -Emax 0.0002 \ -select_mdalign 4 -select_md 4 \ -adpstyle 1 \ -Motifcutoff 0.125 \ -align_short 5 -mdEmax 0.0002 -constraints tmp_2-a.train.cst -constraints_out /export/projects/compbio/experiments/casp4/t123/T0123.t2k_2.cst Reading parameter file tmp_2-b.mod tmp_2-b.mod(18): Reading MODEL -- Final model for run tmp_2-b SAM: hmmscore v3.2 (July 31, 2000) compiled 08/04/00_11:31:33 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -129.770523 Database has 142 sequences with 23829 residues. @@@@ cp -f tmp_2-b.mult /export/projects/compbio/experiments/casp4/t123/T0123.t2k_2.a2m # About to count sequences in prefilter_3.fa @@@@ modelfromalign tmp_3-a -alignfile /export/projects/compbio/experiments/casp4/t123/T0123.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 /export/projects/compbio/experiments/casp4/t123/T0123.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.2 (July 31, 2000) compiled 07/31/00_14:08:27 Reading alignment file /export/projects/compbio/experiments/casp4/t123/T0123.t2k_2.a2m (139 sequences, 160 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 541007 \ -simple_threshold -24.1089427738805 -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.125 -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(21): Reading MODEL -- Model from alignment file /export/projects/compbio/experiments/casp4/t123/T0123.t2k_2.a2m SAM: hmmscore v3.2 (July 31, 2000) compiled 08/04/00_11:31:33 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -111.399551 Database has 172 sequences with 31269 residues. @@@@ sortseq tmp_3-a-sort -alignfile tmp_3-a.mult \ -NLLfile tmp_3-a.mstat -a2mdots 0 SAM: sortseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:48 Reading alignment file tmp_3-a.mult (155 sequences, 160 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=156 nll_thresh=-20.1089427738805 frac_id=0.827187018485051 @@@@ uniqueseq thinned_tmp -alignfile tmp_3-a.train.seq -a protein \ -percent_id 0.827187018485051 -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file tmp_3-a.train.seq (156 sequences, 160 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 20 (of 156) duplicate sequences with differing IDs Dropping 79 (of 156) sequences with > 82.7% id 57 sequences left after dropping 99 of 156 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.2 (July 31, 2000) compiled 07/31/00_14:08:27 Reading alignment file tmp_3-a.train.seq (57 sequences, 160 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(21): Reading MODEL -- Model from alignment file /export/projects/compbio/experiments/casp4/t123/T0123.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.2 (July 31, 2000) compiled 07/31/00_14:08:12 -450.60 4.94 -113.61 53.26 5 0 161 @@@@ 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 541007 \ -db m0.a2m -db tmp_3-a.sel \ -select_score 4 \ -select_seq 4 \ -simple_threshold -24.1089427738805 -Emax 0.001 \ -select_mdalign 4 -select_md 4 \ -adpstyle 1 \ -Motifcutoff 0.125 \ -align_short 5 -mdEmax 0.001 -constraints tmp_3-a.train.cst -constraints_out /export/projects/compbio/experiments/casp4/t123/T0123.t2k_3.cst Reading parameter file tmp_3-b.mod tmp_3-b.mod(18): Reading MODEL -- Final model for run tmp_3-b SAM: hmmscore v3.2 (July 31, 2000) compiled 08/04/00_11:31:33 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -117.236549 Database has 160 sequences with 27030 residues. @@@@ cp -f tmp_3-b.mult /export/projects/compbio/experiments/casp4/t123/T0123.t2k_3.a2m # About to count sequences in prefilter_4.fa @@@@ modelfromalign tmp_4-a -alignfile /export/projects/compbio/experiments/casp4/t123/T0123.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 /export/projects/compbio/experiments/casp4/t123/T0123.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.2 (July 31, 2000) compiled 07/31/00_14:08:27 Reading alignment file /export/projects/compbio/experiments/casp4/t123/T0123.t2k_3.a2m (156 sequences, 160 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 541007 \ -simple_threshold -22.4995048540528 -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.125 -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(21): Reading MODEL -- Model from alignment file /export/projects/compbio/experiments/casp4/t123/T0123.t2k_3.a2m SAM: hmmscore v3.2 (July 31, 2000) compiled 08/04/00_11:31:33 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -51.172310 Database has 423 sequences with 95006 residues. @@@@ sortseq tmp_4-a-sort -alignfile tmp_4-a.mult \ -NLLfile tmp_4-a.mstat -a2mdots 0 SAM: sortseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:48 Reading alignment file tmp_4-a.mult (186 sequences, 160 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=187 nll_thresh=-18.4995048540528 frac_id=0.783156285060481 @@@@ uniqueseq thinned_tmp -alignfile tmp_4-a.train.seq -a protein \ -percent_id 0.783156285060481 -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file tmp_4-a.train.seq (187 sequences, 160 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 27 (of 187) duplicate sequences with differing IDs Dropping 97 (of 187) sequences with > 78.3% id 63 sequences left after dropping 124 of 187 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.2 (July 31, 2000) compiled 07/31/00_14:08:27 Reading alignment file tmp_4-a.train.seq (63 sequences, 160 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(21): Reading MODEL -- Model from alignment file /export/projects/compbio/experiments/casp4/t123/T0123.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.2 (July 31, 2000) compiled 07/31/00_14:08:12 -451.47 3.76 -104.50 51.81 4 0 161 @@@@ 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 541007 \ -db m0.a2m -db tmp_4-a.sel \ -select_score 4 \ -select_seq 4 \ -simple_threshold -22.4995048540528 -Emax 0.005 \ -select_mdalign 4 -select_md 4 \ -adpstyle 1 \ -Motifcutoff 0.125 \ -align_short 5 -mdEmax 0.005 -constraints tmp_4-a.train.cst -constraints_out /export/projects/compbio/experiments/casp4/t123/T0123.t2k_4.cst Reading parameter file tmp_4-b.mod tmp_4-b.mod(18): Reading MODEL -- Final model for run tmp_4-b SAM: hmmscore v3.2 (July 31, 2000) compiled 08/04/00_11:31:33 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -105.258461 Database has 193 sequences with 32874 residues. @@@@ cp -f tmp_4-b.mult /export/projects/compbio/experiments/casp4/t123/T0123.t2k_4.a2m @@@@ modelfromalign tmp_5-a -alignfile /export/projects/compbio/experiments/casp4/t123/T0123.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 /export/projects/compbio/experiments/casp4/t123/T0123.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.2 (July 31, 2000) compiled 07/31/00_14:08:27 Reading alignment file /export/projects/compbio/experiments/casp4/t123/T0123.t2k_4.a2m (190 sequences, 160 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 541007 \ -simple_threshold -22.4995048540528 -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.125 -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(21): Reading MODEL -- Model from alignment file /export/projects/compbio/experiments/casp4/t123/T0123.t2k_4.a2m SAM: hmmscore v3.2 (July 31, 2000) compiled 08/04/00_11:31:33 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -53.904678 Database has 423 sequences with 95006 residues. @@@@ sortseq tmp_5-a-sort -alignfile tmp_5-a.mult \ -NLLfile tmp_5-a.mstat -a2mdots 0 SAM: sortseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:48 Reading alignment file tmp_5-a.mult (182 sequences, 160 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 /export/projects/compbio/experiments/casp4/t123/T0123.t2k.a2m @@@@ cp tmp_5-a.train.cst /export/projects/compbio/experiments/casp4/t123/T0123.t2k.cst @@@@ rm -rf /projects/compbio/tmp/target2k-beta-5581 < /dev/null gzip -f T0123.t2k.a2m T0123.t2k_*.a2m prettyalign T0123.t2k.a2m.gz -m8 -i -n -L3333 > T0123.t2k.pa SAM: prettyalign v3.2 (July 31, 2000) compiled 07/31/00_16:52:38 gunzip -c T0123.t2k.a2m.gz > tmp.a2m a2m2html -a2m_in tmp.a2m -retrieve > T0123.t2k.pa.html SAM: /projects/compbio/bin/alpha/prettyalign v3.2 (July 31, 2000) compiled 07/31/00_16:52:38 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 T0123.t2k.a2m.gz >> tmp.script echo PrintPrediction T0123.t2k.2d 3670-4530-6947 >> tmp.script echo PrintPredictionFasta T0123.t2k.2d.seq >> tmp.script echo PrintRDB T0123.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 T0123.t2k.a2m.gz # Using SequenceWeight EntropyWeight(1.4, 10) # Using regularizer /projects/compbio/lib/recode3.20comp for sequence weight regularizer. # T0123.t2k.a2m with 183 sequences, total weight= 108.694 avg weight= 0.593957 clipped 46 iterations # AdjustWeights couldn't save exactly 1.4 bits/position, saving 1.2649 bits. Regularizing alignment for T0123.t2k.a2m.gz # After reading T0123.t2k.a2m.gz, have 160 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 T0123.t2k.2d # command:# Network initialization done # Printing prediction in FASTA format to T0123.t2k.2d.seq # command:# Network initialization done # Printing prediction to T0123.t2k.2d.rdb # command:rm tmp.script gunzip -f T0123.t2k.a2m.gz /projects/compbio/usr/karplus/src/phytree/phytree -o -i -r flat T0123.t2k T0123.t2k.a2m /projects/compbio/lib/recode2.20comp Prior library /projects/compbio/lib/recode2.20comp read. Reading alignment file T0123.t2k.a2m (183 sequences) as A2M alignment. Alignment T0123.t2k.a2m read. Adding sequences:0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 One dot for every 10 merges: .................. Merges done. CPU = user 0: 3: 0 system 0: 0: 0 gzip -f T0123.t2k.a2m gzip -f T0123.t2k_sorted.a2m rm T0123.t2k.phytrace T0123.t2k.tree_weight prettyalign T0123.t2k_sorted.a2m.gz -m8 -i -n -L3333 > T0123.t2k_sorted.pa SAM: prettyalign v3.2 (July 31, 2000) compiled 07/31/00_16:52:38 /projects/compbio/usr/karplus/src/phytree/dg T0123.t2k.tree /projects/compbio/usr/karplus/src/phytree/dtree T0123.t2k.tree make[1]: Leaving directory `/export/projects/compbio/experiments/casp4/t123' echo "WeightingRegularizer /projects/compbio/lib/recode3.20comp" > tmp.script echo "SequenceWeight EntropyWeight 0.5 10.0" >> tmp.script echo "ReadA2M T0123.t2k.a2m.gz" >> tmp.script echo "SavingsRegularizer /projects/compbio/lib/recode3.20comp" >> tmp.script echo "PrintSavings T0123.t2k.saves" >> tmp.script echo "exit" >> tmp.script estimate-dist < tmp.script # command:Weighting regularizer set to recode3.20comp # command:Sequence weighting on input set to EntropyWeight(0.5, 10) # command:Reading A2M format from T0123.t2k.a2m.gz % using regularizer recode3.20comp % T0123.t2k.a2m.gz with 183 sequences, total weight= 7.46556 avg weight= 0.0407954 After reading T0123.t2k.a2m.gz, have 160 columns, totalling 1115.82 Using SequenceWeight EntropyWeight(0.5, 10) Using regularizer recode3.20compfor sequence weight regularizer. # command:Savings regularizer set to recode3.20comp # command:% Computing savings using regularizer recode3.20comp Printed 160 columns to T0123.t2k.saves # command:rm tmp.script if test -e T0123.t2k-logo.labels ;\ then makelogo -b 60 -t T0123.t2k -o T0123.t2k-logo.eps -l T0123.t2k-logo.labels T0123.t2k.saves ;\ else makelogo -b 60 -t T0123.t2k -o T0123.t2k-logo.eps -lf T0123.t2k.2d.seq T0123.t2k.saves ;\ fi /projects/compbio/experiments/models.97/scripts99/score-target -alignment T0123.t2k.a2m.gz -build w0.5 \ -db /projects/compbio/data/pdb/all-protein \ -Emax 40.0 \ -viterbi 0 -sw_score 2 \ -output T0123.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-18818 @@@@ mkdir /projects/compbio/tmp/score-target-18818/for-scorescript @@@@ /projects/compbio/experiments/models.97/scripts99/score-mods-vs-seqs.pl -root_tmp /projects/compbio/tmp/score-target-18818/for-scorescript -a2msfile /projects/compbio/tmp/score-target-18818/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-18818/score-target.op @@@@ chgrp protein /projects/compbio/tmp/score-target-18818/for-scorescript/split-tmp @@@@ chgrp protein /projects/compbio/tmp/score-target-18818/for-scorescript/split-tmp @@@@ /projects/compbio/bin/scripts/w0.5 /export/projects/compbio/experiments/casp4/t123/T0123.t2k.a2m.gz /projects/compbio/tmp/score-target-18818/for-scorescript/split-tmp/beta-16969/T0123.t2k.mod Reading /projects/compbio/bin/scripts/sam-t99.conf @@@@ chgrp protein /projects/compbio/tmp/tmp-build-weighted-model-beta-6514 @@@@ uniqueseq /projects/compbio/tmp/tmp-build-weighted-model-beta-6514/tmp -alignfile /export/projects/compbio/experiments/casp4/t123/T0123.t2k.a2m.gz \ -a protein -percent_id 0.8 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file /export/projects/compbio/experiments/casp4/t123/T0123.t2k.a2m.gz (183 sequences, 160 columns) as A2M alignment. Writing sequence output to /projects/compbio/tmp/tmp-build-weighted-model-beta-6514/tmp.a2m. Dropping 22 (of 183) duplicate sequences with differing IDs Dropping 98 (of 183) sequences with > 80.0% id 63 sequences left after dropping 120 of 183 sequences. @@@@ modelfromalign /projects/compbio/tmp/tmp-build-weighted-model-beta-6514/tmp -alignfile /projects/compbio/tmp/tmp-build-weighted-model-beta-6514/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.2 (July 31, 2000) compiled 07/31/00_14:08:27 Reading alignment file /projects/compbio/tmp/tmp-build-weighted-model-beta-6514/tmp.a2m (63 sequences, 160 columns) as A2M alignment. @@@@ cp -f /projects/compbio/tmp/tmp-build-weighted-model-beta-6514/tmp.mod /projects/compbio/tmp/score-target-18818/for-scorescript/split-tmp/beta-16969/T0123.t2k.mod @@@@ rm -rf /projects/compbio/tmp/tmp-build-weighted-model-beta-6514