make -k T0122-blast.rdb T0122-double-blast.rdb make[1]: Entering directory `/export/projects/compbio/experiments/casp4/t122' /projects/compbio/experiments/casp4/scripts/single-blast -q T0122.seq -interdb /projects/compbio/data/nrp/nr -db /projects/compbio/data/pdb/all-protein > T0122-blast.rdb starting T0122 with E=10 /projects/compbio/experiments/casp4/scripts/double-blast -q T0122.seq -interdb /projects/compbio/data/nrp/nr -db /projects/compbio/data/pdb/all-protein > T0122-double-blast.rdb starting T0122 with E=0.000050 T0122 T0122 blasting: T0122 gi|4676|emb|CAA24636.1|_19:167 blasting: T0122 gi|401211|sp|P31204|TRPA_ANTSP_16:236 blasting: T0122 gi|136255|sp|P06562|TRPA_BRELA_17:240 blasting: T0122 gi|136257|sp|P00928|TRPA_ECOLI_15:252 blasting: T0122 gi|136259|sp|P00930|TRPA_KLEAE_15:241 blasting: T0122 gi|136264|sp|P11081|TRPA_PSEPU_18:242 blasting: T0122 gi|79917|pir||S06782_8:86 blasting: T0122 gi|136265|sp|P00929|TRPA_SALTY_15:243 blasting: T0122 gi|136267|sp|P22095|TRPA_VIBPA_17:259 blasting: T0122 gi|68331|pir||TSAEAA_15:241 blasting: T0122 gi|68332|pir||TSPSAA_18:250 blasting: T0122 gi|136253|sp|P19867|TRPA_BACST_12:235 blasting: T0122 gi|136254|sp|P07601|TRPA_BACSU_14:235 blasting: T0122 gi|136260|sp|P17166|TRPA_LACCA_9:253 blasting: T0122 gi|2822116|sp|P16578|TRP_COPCI_25:262 blasting: T0122 gi|94733|pir||B35407_17:248 blasting: T0122 gi|136261|sp|P26920|TRPA_METTM_14:254 blasting: T0122 gi|136258|sp|P18284|TRPA_HALVO_7:240 blasting: T0122 gi|136372|sp|P13228|TRP_NEUCR_18:255 blasting: T0122 gi|136256|sp|P12291|TRPA_CAUCR_22:256 blasting: T0122 gi|136262|sp|P14637|TRPA_METVO_26:267 blasting: T0122 gi|136263|sp|P07344|TRPA_PSEAE_18:251 blasting: T0122 gi|136266|sp|P16608|TRPA_THETH_17:248 blasting: T0122 gi|144100|gb|AAA83990.1|_1:80 blasting: T0122 gi|267167|sp|Q01997|TRPA_LACLA_11:251 blasting: T0122 gi|464911|sp|P34816|TRPA_PSESY_18:237 blasting: T0122 gi|1174781|sp|P42389|TRPA_BUCAP_15:255 blasting: T0122 gi|409540|gb|AAC60451.1|_12:179 blasting: T0122 gi|409778|dbj|BAA04617.1|_3:233 blasting: T0122 gi|421534|pir||PN0500_12:179 blasting: T0122 gi|1174783|sp|P42390|TRPA_MAIZE_103:323 blasting: T0122 gi|541528|pir||S41995_3:233 blasting: T0122 gi|2129755|pir||S59519_65:303 blasting: T0122 gi|7674382|sp|Q44604|TRPA_BUCSC_14:253 blasting: T0122 gi|775154|gb|AAA65158.1|_15:252 blasting: T0122 gi|775168|gb|AAB60041.1|_15:252 blasting: T0122 gi|775174|gb|AAA65170.1|_15:252 blasting: T0122 gi|775181|gb|AAB60046.1|_15:252 blasting: T0122 gi|775188|gb|AAB60052.1|_15:252 blasting: T0122 gi|775193|gb|AAB60056.1|_15:252 blasting: T0122 gi|775198|gb|AAA65176.1|_15:252 blasting: T0122 gi|1717757|sp|P50382|TRPA_SULSO_5:239 blasting: T0122 gi|2120372|pir||S59050_5:238 blasting: T0122 gi|1174782|sp|P43759|TRPA_HAEIN_16:242 blasting: T0122 gi|1717756|sp|P51382|TRPA_PORPU_23:246 blasting: T0122 gi|2129331|pir||E64429_22:267 blasting: T0122 gi|2501410|sp|P77960|TRPA_SYNY3_18:242 blasting: T0122 gi|2098385|pdb|2TSY|A_15:237 blasting: T0122 gi|6226271|sp|O06130|TRPA_MYCTU_26:248 blasting: T0122 gi|2501408|sp|P56141|TRPA_HELPY_17:258 blasting: T0122 gi|3915157|sp|O13831|TRP_SCHPO_20:266 blasting: T0122 gi|3915160|sp|O22018|TRPA_CYAME_2:226 blasting: T0122 gi|3915890|sp|P34793|TRPA_CYACA_3:233 blasting: T0122 gi|7436988|pir||A69089_15:264 blasting: T0122 gi|3334388|sp|O28673|TRPA_ARCFU_3:243 blasting: T0122 gi|6226268|sp|O67502|TRPA_AQUAE_15:258 blasting: T0122 gi|6226270|sp|O68906|TRPA_MYCIT_21:260 blasting: T0122 gi|6226272|sp|O68816|TRPA_STRCO_22:241 blasting: T0122 gi|7674377|sp|O68429|TRPA_BUCDN_15:259 blasting: T0122 gi|3212365|pdb|1A5A|A_15:243 blasting: T0122 gi|7674379|sp|O84173|TRPA_CHLTR_7:214 blasting: T0122 gi|7436987|pir||T01088_27:273 blasting: T0122 gi|4151263|gb|AAD04348.1|_18:129 blasting: T0122 gi|4151265|gb|AAD04349.1|_18:129 blasting: T0122 gi|4151267|gb|AAD04350.1|_18:129 blasting: T0122 gi|4151275|gb|AAD04354.1|_18:131 blasting: T0122 gi|4151277|gb|AAD04355.1|_1:133 blasting: T0122 gi|4151279|gb|AAD04356.1|_1:133 blasting: T0122 gi|4151281|gb|AAD04357.1|_1:133 blasting: T0122 gi|4151283|gb|AAD04358.1|_1:133 blasting: T0122 gi|7674400|sp|Q9ZJV0|TRPA_HELPJ_14:258 blasting: T0122 gi|7674390|sp|Q9X4E8|TRPA_RHOSH_18:252 blasting: T0122 gi|6226718|sp|P50908|TRPA_THEMA_5:238 blasting: T0122 gi|7674392|sp|Q9Y8T3|TRPA_AERPE_11:232 blasting: T0122 gi|5231181|gb|AAD41107.1|AF157817_2_18:217 blasting: T0122 gi|5231184|gb|AAD41109.1|AF157818_2_18:217 blasting: T0122 gi|5231187|gb|AAD41111.1|AF157819_2_18:217 blasting: T0122 gi|5231190|gb|AAD41113.1|AF157820_2_18:217 blasting: T0122 gi|5231193|gb|AAD41115.1|AF157821_2_18:217 blasting: T0122 gi|5231196|gb|AAD41117.1|AF157822_2_18:217 blasting: T0122 gi|5231199|gb|AAD41119.1|AF157823_2_18:217 blasting: T0122 gi|5231202|gb|AAD41121.1|AF157824_2_18:217 blasting: T0122 gi|5231205|gb|AAD41123.1|AF157825_2_18:217 blasting: T0122 gi|5231208|gb|AAD41125.1|AF157826_2_18:220 blasting: T0122 gi|7436994|pir||E75162_1:248 blasting: T0122 gi|7674396|sp|Q9YGA9|TRPA_PYRKO_1:246 blasting: T0122 gi|5764615|gb|AAD51339.1|AF173835_3_18:261 blasting: T0122 gi|5834682|emb|CAB55325.1|_19:76 blasting: T0122 gi|7474050|pir||E75455_30:256 blasting: T0122 gi|6466362|gb|AAF12944.1|AF022186_67_17:244 blasting: T0122 gi|6492284|gb|AAF14255.1|AF132318_5_14:254 blasting: T0122 gi|6599047|emb|CAB63648.1|_15:174 blasting: T0122 gi|6599049|emb|CAB63649.1|_15:174 blasting: T0122 gi|6689543|emb|CAA79517.1|_17:239 blasting: T0122 gi|6967824|emb|CAB74186.1|_8:246 blasting: T0122 gi|7225905|gb|AAF41096.1|_12:246 blasting: T0122 gi|7379594|emb|CAB84159.1|_12:246 blasting: T0122 gi|7381257|gb|AAF61458.1|AF139661_3_29:270 blasting: T0122 gi|7676165|sp|O27697|TRPA_METTH_12:261 blasting: T0122 gi|7676173|sp|Q60180|TRPA_METJA_19:264 blasting: T0122 gi|8272441|dbj|BAA96470.1|_1:188 blasting: T0122 gi|9106381|gb|AAF84185.1|AE003969_2_18:242 blasting: T0122 gi|9651535|gb|AAF91181.1|AF207903_2_15:243 blasting: T0122 gi|9655645|gb|AAF94328.1|_15:241 blasting: make[1]: Leaving directory `/export/projects/compbio/experiments/casp4/t122' make -k T0122.t2k.a2m.gz \ T0122.t2k.pa \ T0122.t2k.pa.html \ T0122.t2k.2d \ T0122.t2k.tree \ T0122.t2k_sorted.pa \ T0122.t2k.tree.ps \ T0122.t2k.tree-unroot.ps make[1]: Entering directory `/export/projects/compbio/experiments/casp4/t122' /projects/compbio/experiments/models.97/scripts2k/target2k -out T0122.t2k \ -seed T0122.seq \ -all @@@@ chgrp protein /projects/compbio/tmp/target2k-beta-5861 @@@@ cp /export/projects/compbio/experiments/casp4/t122/T0122.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-5079 /projects/compbio/experiments/models.97/scripts2k/blast-prefilter: E-value thresholds used are For /export/projects/compbio/tmp/target2k-beta-5861/prefilter_1.fa, 0.01 For /export/projects/compbio/tmp/target2k-beta-5861/prefilter_2.fa, 1 For /export/projects/compbio/tmp/target2k-beta-5861/prefilter_3.fa, 10 For /export/projects/compbio/tmp/target2k-beta-5861/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-5861/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-5861/m0.a2m (1 sequences, 248 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-5079 < /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, 248 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, 248 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.0806451612903226 -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: -116.872368 Database has 105 sequences with 27925 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 (104 sequences, 248 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=105 nll_thresh=-22.4115278685381 frac_id=0.915363956825628 @@@@ uniqueseq thinned_tmp -alignfile tmp_1-a.train.seq -a protein \ -percent_id 0.915363956825628 -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file tmp_1-a.train.seq (105 sequences, 248 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 4 (of 105) duplicate sequences with differing IDs Dropping 36 (of 105) sequences with > 91.5% id 65 sequences left after dropping 40 of 105 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 (65 sequences, 248 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 -699.08 108.53 -222.58 80.33 13 0 249 @@@@ 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.0806451612903226 \ -align_short 5 -mdEmax 0.0001 -constraints tmp_1-a.train.cst -constraints_out /export/projects/compbio/experiments/casp4/t122/T0122.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: -223.389709 Database has 106 sequences with 28173 residues. @@@@ cp -f tmp_1-b.mult /export/projects/compbio/experiments/casp4/t122/T0122.t2k_1.a2m # About to count sequences in prefilter_2.fa @@@@ modelfromalign tmp_2-a -alignfile /export/projects/compbio/experiments/casp4/t122/T0122.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/t122/T0122.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/t122/T0122.t2k_1.a2m (106 sequences, 248 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.0806451612903226 -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/t122/T0122.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: -199.042938 Database has 112 sequences with 29433 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 (107 sequences, 248 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=108 nll_thresh=-21.7183806877933 frac_id=0.905554508734834 @@@@ uniqueseq thinned_tmp -alignfile tmp_2-a.train.seq -a protein \ -percent_id 0.905554508734834 -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file tmp_2-a.train.seq (108 sequences, 248 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 2 (of 108) duplicate sequences with differing IDs Dropping 40 (of 108) sequences with > 90.6% id 66 sequences left after dropping 42 of 108 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 (66 sequences, 248 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/t122/T0122.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 -697.12 10.74 -210.11 77.29 9 0 249 @@@@ 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.0806451612903226 \ -align_short 5 -mdEmax 0.0002 -constraints tmp_2-a.train.cst -constraints_out /export/projects/compbio/experiments/casp4/t122/T0122.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: -202.185303 Database has 109 sequences with 28716 residues. @@@@ cp -f tmp_2-b.mult /export/projects/compbio/experiments/casp4/t122/T0122.t2k_2.a2m # About to count sequences in prefilter_3.fa @@@@ modelfromalign tmp_3-a -alignfile /export/projects/compbio/experiments/casp4/t122/T0122.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/t122/T0122.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/t122/T0122.t2k_2.a2m (108 sequences, 248 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.0806451612903226 -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/t122/T0122.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: -160.600464 Database has 136 sequences with 36221 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 (107 sequences, 248 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=108 nll_thresh=-20.1089427738805 frac_id=0.892418248962251 @@@@ uniqueseq thinned_tmp -alignfile tmp_3-a.train.seq -a protein \ -percent_id 0.892418248962251 -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file tmp_3-a.train.seq (108 sequences, 248 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 2 (of 108) duplicate sequences with differing IDs Dropping 40 (of 108) sequences with > 89.2% id 66 sequences left after dropping 42 of 108 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 (66 sequences, 248 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/t122/T0122.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 -697.40 12.76 -199.03 77.35 7 0 249 @@@@ 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.0806451612903226 \ -align_short 5 -mdEmax 0.001 -constraints tmp_3-a.train.cst -constraints_out /export/projects/compbio/experiments/casp4/t122/T0122.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: -185.590057 Database has 113 sequences with 29458 residues. @@@@ cp -f tmp_3-b.mult /export/projects/compbio/experiments/casp4/t122/T0122.t2k_3.a2m # About to count sequences in prefilter_4.fa @@@@ modelfromalign tmp_4-a -alignfile /export/projects/compbio/experiments/casp4/t122/T0122.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/t122/T0122.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/t122/T0122.t2k_3.a2m (109 sequences, 248 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.0806451612903226 -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/t122/T0122.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: -57.845978 Database has 407 sequences with 109770 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 (112 sequences, 248 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=113 nll_thresh=-18.4995048540528 frac_id=0.86971461097418 @@@@ uniqueseq thinned_tmp -alignfile tmp_4-a.train.seq -a protein \ -percent_id 0.86971461097418 -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file tmp_4-a.train.seq (113 sequences, 248 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 2 (of 113) duplicate sequences with differing IDs Dropping 41 (of 113) sequences with > 87.0% id 70 sequences left after dropping 43 of 113 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 (70 sequences, 248 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/t122/T0122.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 -697.70 -0.19 -176.51 85.64 10 0 249 @@@@ 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.0806451612903226 \ -align_short 5 -mdEmax 0.005 -constraints tmp_4-a.train.cst -constraints_out /export/projects/compbio/experiments/casp4/t122/T0122.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: -139.786423 Database has 148 sequences with 38427 residues. @@@@ cp -f tmp_4-b.mult /export/projects/compbio/experiments/casp4/t122/T0122.t2k_4.a2m @@@@ modelfromalign tmp_5-a -alignfile /export/projects/compbio/experiments/casp4/t122/T0122.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/t122/T0122.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/t122/T0122.t2k_4.a2m (126 sequences, 248 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.0806451612903226 -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/t122/T0122.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: -58.782566 Database has 407 sequences with 109770 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 (117 sequences, 248 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/t122/T0122.t2k.a2m @@@@ cp tmp_5-a.train.cst /export/projects/compbio/experiments/casp4/t122/T0122.t2k.cst @@@@ rm -rf /projects/compbio/tmp/target2k-beta-5861 < /dev/null gzip -f T0122.t2k.a2m T0122.t2k_*.a2m prettyalign T0122.t2k.a2m.gz -m8 -i -n -L3333 > T0122.t2k.pa SAM: prettyalign v3.2 (July 31, 2000) compiled 07/31/00_16:52:38 gunzip -c T0122.t2k.a2m.gz > tmp.a2m a2m2html -a2m_in tmp.a2m -retrieve > T0122.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 T0122.t2k.a2m.gz >> tmp.script echo PrintPrediction T0122.t2k.2d 3670-4530-6947 >> tmp.script echo PrintPredictionFasta T0122.t2k.2d.seq >> tmp.script echo PrintRDB T0122.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 T0122.t2k.a2m.gz # Using SequenceWeight EntropyWeight(1.4, 10) # Using regularizer /projects/compbio/lib/recode3.20comp for sequence weight regularizer. # T0122.t2k.a2m with 118 sequences, total weight= 33.032 avg weight= 0.279932 46 iterations Regularizing alignment for T0122.t2k.a2m.gz # After reading T0122.t2k.a2m.gz, have 248 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 T0122.t2k.2d # command:# Network initialization done # Printing prediction in FASTA format to T0122.t2k.2d.seq # command:# Network initialization done # Printing prediction to T0122.t2k.2d.rdb # command:rm tmp.script gunzip -f T0122.t2k.a2m.gz /projects/compbio/usr/karplus/src/phytree/phytree -o -i -r flat T0122.t2k T0122.t2k.a2m /projects/compbio/lib/recode2.20comp Prior library /projects/compbio/lib/recode2.20comp read. Reading alignment file T0122.t2k.a2m (118 sequences) as A2M alignment. Alignment T0122.t2k.a2m read. Adding sequences:0 10 20 30 40 50 60 70 80 90 100 110 One dot for every 10 merges: ............ Merges done. CPU = user 0: 2:17 system 0: 0: 1 gzip -f T0122.t2k.a2m gzip -f T0122.t2k_sorted.a2m rm T0122.t2k.phytrace T0122.t2k.tree_weight prettyalign T0122.t2k_sorted.a2m.gz -m8 -i -n -L3333 > T0122.t2k_sorted.pa SAM: prettyalign v3.2 (July 31, 2000) compiled 07/31/00_16:52:38 /projects/compbio/usr/karplus/src/phytree/dg T0122.t2k.tree /projects/compbio/usr/karplus/src/phytree/dtree T0122.t2k.tree make[1]: Leaving directory `/export/projects/compbio/experiments/casp4/t122' echo "WeightingRegularizer /projects/compbio/lib/recode3.20comp" > tmp.script echo "SequenceWeight EntropyWeight 0.5 10.0" >> tmp.script echo "ReadA2M T0122.t2k.a2m.gz" >> tmp.script echo "SavingsRegularizer /projects/compbio/lib/recode3.20comp" >> tmp.script echo "PrintSavings T0122.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 T0122.t2k.a2m.gz % using regularizer recode3.20comp % T0122.t2k.a2m.gz with 118 sequences, total weight= 3.88155 avg weight= 0.0328945 After reading T0122.t2k.a2m.gz, have 248 columns, totalling 750.092 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 248 columns to T0122.t2k.saves # command:rm tmp.script if test -e T0122.t2k-logo.labels ;\ then makelogo -b 60 -t T0122.t2k -o T0122.t2k-logo.eps -l T0122.t2k-logo.labels T0122.t2k.saves ;\ else makelogo -b 60 -t T0122.t2k -o T0122.t2k-logo.eps -lf T0122.t2k.2d.seq T0122.t2k.saves ;\ fi /projects/compbio/experiments/models.97/scripts99/score-target -alignment T0122.t2k.a2m.gz -build w0.5 \ -db /projects/compbio/data/pdb/all-protein \ -Emax 40.0 \ -viterbi 0 -sw_score 2 \ -output T0122.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-7150 @@@@ mkdir /projects/compbio/tmp/score-target-7150/for-scorescript @@@@ /projects/compbio/experiments/models.97/scripts99/score-mods-vs-seqs.pl -root_tmp /projects/compbio/tmp/score-target-7150/for-scorescript -a2msfile /projects/compbio/tmp/score-target-7150/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-7150/score-target.op @@@@ chgrp protein /projects/compbio/tmp/score-target-7150/for-scorescript/split-tmp @@@@ chgrp protein /projects/compbio/tmp/score-target-7150/for-scorescript/split-tmp @@@@ /projects/compbio/bin/scripts/w0.5 /export/projects/compbio/experiments/casp4/t122/T0122.t2k.a2m.gz /projects/compbio/tmp/score-target-7150/for-scorescript/split-tmp/beta-2432/T0122.t2k.mod Reading /projects/compbio/bin/scripts/sam-t99.conf @@@@ chgrp protein /projects/compbio/tmp/tmp-build-weighted-model-beta-28295 @@@@ uniqueseq /projects/compbio/tmp/tmp-build-weighted-model-beta-28295/tmp -alignfile /export/projects/compbio/experiments/casp4/t122/T0122.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/t122/T0122.t2k.a2m.gz (118 sequences, 248 columns) as A2M alignment. Writing sequence output to /projects/compbio/tmp/tmp-build-weighted-model-beta-28295/tmp.a2m. Dropping 4 (of 118) duplicate sequences with differing IDs Dropping 48 (of 118) sequences with > 80.0% id 66 sequences left after dropping 52 of 118 sequences. @@@@ modelfromalign /projects/compbio/tmp/tmp-build-weighted-model-beta-28295/tmp -alignfile /projects/compbio/tmp/tmp-build-weighted-model-beta-28295/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-28295/tmp.a2m (66 sequences, 248 columns) as A2M alignment. @@@@ cp -f /projects/compbio/tmp/tmp-build-weighted-model-beta-28295/tmp.mod /projects/compbio/tmp/score-target-7150/for-scorescript/split-tmp/beta-2432/T0122.t2k.mod @@@@ rm -rf /projects/compbio/tmp/tmp-build-weighted-model-beta-28295