make -k T0116-500-end-blast.rdb T0116-500-end-double-blast.rdb make[1]: Entering directory `/export/projects/compbio/experiments/casp4/t116/t116-500-end' /projects/compbio/experiments/casp4/scripts/single-blast -q T0116-500-end.seq -interdb /projects/compbio/data/nrp/nr -db /projects/compbio/data/pdb/all-protein > T0116-500-end-blast.rdb starting T0116-500-end with E=10 /projects/compbio/experiments/casp4/scripts/double-blast -q T0116-500-end.seq -interdb /projects/compbio/data/nrp/nr -db /projects/compbio/data/pdb/all-protein > T0116-500-end-double-blast.rdb starting T0116-500-end with E=0.000050 T0116-500-end T0116-500-end blasting: T0116-500-end gi|135075|sp|P26359|SWI4_SCHPO_671:967 blasting: T0116-500-end gi|79102|pir||A28668_542:832 blasting: T0116-500-end gi|123080|sp|P10564|HEXA_STRPN_513:781 blasting: T0116-500-end gi|123080|sp|P10564|HEXA_STRPN_816:841 blasting: T0116-500-end gi|118951|sp|P20585|MSH3_HUMAN_807:1011 blasting: T0116-500-end gi|118951|sp|P20585|MSH3_HUMAN_1053:1086 blasting: T0116-500-end gi|110917|pir||A32495_776:925 blasting: T0116-500-end gi|127555|sp|P27345|MUTS_AZOVI_526:805 blasting: T0116-500-end gi|127555|sp|P27345|MUTS_AZOVI_827:852 blasting: T0116-500-end gi|127556|sp|P23909|MUTS_ECOLI_540:832 blasting: T0116-500-end gi|127556|sp|P23909|MUTS_ECOLI_826:852 blasting: T0116-500-end gi|172000|gb|AAA34801.1|_692:958 blasting: T0116-500-end gi|172002|gb|AAA34802.1|_616:873 blasting: T0116-500-end gi|400971|sp|P13705|MSH3_MOUSE_776:965 blasting: T0116-500-end gi|400971|sp|P13705|MSH3_MOUSE_1007:1040 blasting: T0116-500-end gi|387849|gb|AAB60711.1|_776:965 blasting: T0116-500-end gi|387849|gb|AAB60711.1|_1007:1040 blasting: T0116-500-end gi|417330|sp|P10339|MUTS_SALTY_543:833 blasting: T0116-500-end gi|448278|prf||1916402A_746:1043 blasting: T0116-500-end gi|730066|sp|P40965|MSH4_YEAST_552:823 blasting: T0116-500-end gi|1174416|sp|P43248|MSH2_DROME_593:839 blasting: T0116-500-end gi|726086|gb|AAA75027.1|_590:845 blasting: T0116-500-end gi|7448104|pir||I64827_590:787 blasting: T0116-500-end gi|1002520|gb|AAB19235.1|_508:778 blasting: T0116-500-end gi|1171081|sp|P44834|MUTS_HAEIN_543:795 blasting: T0116-500-end gi|1171081|sp|P44834|MUTS_HAEIN_826:859 blasting: T0116-500-end gi|1076927|pir||S22569_579:851 blasting: T0116-500-end gi|1079288|pir||S53609_597:845 blasting: T0116-500-end gi|2135744|pir||I37550_590:849 blasting: T0116-500-end gi|1082386|pir||A56742_1021:1247 blasting: T0116-500-end gi|1083296|pir||JC4019_811:1000 blasting: T0116-500-end gi|1083296|pir||JC4019_1042:1075 blasting: T0116-500-end gi|7506963|pir||T16784_1:78 blasting: T0116-500-end gi|1709122|sp|P54275|MSH2_RAT_590:849 blasting: T0116-500-end gi|2497995|sp|Q56215|MUTS_THEAQ_501:811 blasting: T0116-500-end gi|1582703|prf||2119252A_638:762 blasting: T0116-500-end gi|1582703|prf||2119252A_791:841 blasting: T0116-500-end gi|1582703|prf||2119252A_543:616 blasting: T0116-500-end gi|1588283|prf||2208298A_943:1179 blasting: T0116-500-end gi|1592569|gb|AAB97931.1|_540:832 blasting: T0116-500-end gi|1592569|gb|AAB97931.1|_826:852 blasting: T0116-500-end gi|3024198|sp|P70755|MUTS_AQUPY_520:831 blasting: T0116-500-end gi|1619909|gb|AAB16999.1|_503:768 blasting: T0116-500-end gi|3914083|sp|P73625|MUS2_SYNY3_255:428 blasting: T0116-500-end gi|3914083|sp|P73625|MUS2_SYNY3_441:540 blasting: T0116-500-end gi|7448108|pir||S74860_597:869 blasting: T0116-500-end gi|1709123|sp|P52701|MSH6_HUMAN_1089:1315 blasting: T0116-500-end gi|3914084|sp|P94545|MUS2_BACSU_257:544 blasting: T0116-500-end gi|2506881|sp|P54276|MSH6_MOUSE_1087:1313 blasting: T0116-500-end gi|1871501|dbj|BAA09880.1|_509:819 blasting: T0116-500-end gi|1871501|dbj|BAA09880.1|_347:370 blasting: T0116-500-end gi|2104549|gb|AAB57798.1|_1088:1301 blasting: T0116-500-end gi|2104549|gb|AAB57798.1|_1036:1108 blasting: T0116-500-end gi|2145507|pir||S62790_508:818 blasting: T0116-500-end gi|2145507|pir||S62790_346:369 blasting: T0116-500-end gi|2147739|pir||S58881_618:857 blasting: T0116-500-end gi|3914057|sp|O25338|MUS2_HELPY_277:506 blasting: T0116-500-end gi|7448107|pir||T08098_1:116 blasting: T0116-500-end gi|7522430|pir||T38256_508:767 blasting: T0116-500-end gi|2497996|sp|Q56239|MUTS_THETH_509:819 blasting: T0116-500-end gi|2497996|sp|Q56239|MUTS_THETH_347:370 blasting: T0116-500-end gi|3914056|sp|O24617|MSH2_ARATH_573:851 blasting: T0116-500-end gi|2522394|gb|AAB82654.1|_1:86 blasting: T0116-500-end gi|3914053|sp|O13396|MSH2_NEUCR_570:794 blasting: T0116-500-end gi|3914053|sp|O13396|MSH2_NEUCR_823:853 blasting: T0116-500-end gi|7482290|pir||C69102_376:569 blasting: T0116-500-end gi|7404403|sp|P49849|MUTS_BACSU_514:784 blasting: T0116-500-end gi|3914072|sp|O51125|MUS2_BORBU_281:549 blasting: T0116-500-end gi|3914075|sp|O51737|MUTS_BORBU_518:806 blasting: T0116-500-end gi|3914075|sp|O51737|MUTS_BORBU_836:853 blasting: T0116-500-end gi|7512474|pir||JC5839_1089:1315 blasting: T0116-500-end gi|3914062|sp|O65607|MSH3_ARATH_755:893 blasting: T0116-500-end gi|3914062|sp|O65607|MSH3_ARATH_890:1011 blasting: T0116-500-end gi|3914079|sp|O66652|MUTS_AQUAE_520:840 blasting: T0116-500-end gi|3914080|sp|O67287|MUS2_AQUAE_242:545 blasting: T0116-500-end gi|3108220|gb|AAC62533.1|_515:772 blasting: T0116-500-end gi|3914061|sp|O63852|MSHM_SARGL_618:857 blasting: T0116-500-end gi|3319876|emb|CAA07684.1|_742:1016 blasting: T0116-500-end gi|6093453|sp|O83348|MUTS_TREPA_532:813 blasting: T0116-500-end gi|7387913|sp|O84797|MUTS_CHLTR_551:799 blasting: T0116-500-end gi|7492209|pir||T41262_955:1254 blasting: T0116-500-end gi|7490470|pir||T39976_135:365 blasting: T0116-500-end gi|7387927|sp|Q9ZEH5|MUS2_STAAU_202:450 blasting: T0116-500-end gi|3831701|gb|AAC70065.1|_551:828 blasting: T0116-500-end gi|7493777|pir||T18246_512:743 blasting: T0116-500-end gi|6647605|sp|Q9ZDM9|MUTS_RICPR_555:870 blasting: T0116-500-end gi|6647605|sp|Q9ZDM9|MUTS_RICPR_865:890 blasting: T0116-500-end gi|3914066|sp|Q19272|MSH5_CAEEL_652:929 blasting: T0116-500-end gi|3886083|gb|AAC78226.1|_548:820 blasting: T0116-500-end gi|6226648|sp|O04716|MSH6_ARATH_1050:1263 blasting: T0116-500-end gi|6226648|sp|O04716|MSH6_ARATH_998:1070 blasting: T0116-500-end gi|3986757|gb|AAC84154.1|_514:771 blasting: T0116-500-end gi|4034921|emb|CAA10477.1|_1:74 blasting: T0116-500-end gi|4034925|emb|CAA10478.1|_1:74 blasting: T0116-500-end gi|4034929|emb|CAA10479.1|_1:74 blasting: T0116-500-end gi|4099512|gb|AAD00647.1|_1:117 blasting: T0116-500-end gi|4100661|gb|AAD00921.1|_1:130 blasting: T0116-500-end gi|7387928|sp|Q9ZIX6|MUTS_THECA_507:817 blasting: T0116-500-end gi|4139230|gb|AAD04176.1|_573:851 blasting: T0116-500-end gi|4151103|emb|CAA07342.1|_683:913 blasting: T0116-500-end gi|7387931|sp|Q9ZLL4|MUS2_HELPJ_277:548 blasting: T0116-500-end gi|4165574|emb|CAA10614.1|_1:90 blasting: T0116-500-end gi|4336728|gb|AAD17920.1|_620:864 blasting: T0116-500-end gi|7387924|sp|Q9Z6W5|MUTS_CHLPN_557:808 blasting: T0116-500-end gi|4587688|gb|AAD25857.1|_1:98 blasting: T0116-500-end gi|4587694|gb|AAD25860.1|_1:98 blasting: T0116-500-end gi|4587698|gb|AAD25862.1|_1:98 blasting: T0116-500-end gi|4587700|gb|AAD25863.1|_1:98 blasting: T0116-500-end gi|4587702|gb|AAD25864.1|_1:98 blasting: T0116-500-end gi|4587708|gb|AAD25867.1|_1:98 blasting: T0116-500-end gi|4587710|gb|AAD25868.1|_1:98 blasting: T0116-500-end gi|4775576|emb|CAB42554.1|_587:855 blasting: T0116-500-end gi|4775578|emb|CAB42555.1|_892:1128 blasting: T0116-500-end gi|4775580|emb|CAB42556.1|_585:821 blasting: T0116-500-end gi|4887554|emb|CAB43497.1|_540:832 blasting: T0116-500-end gi|4887554|emb|CAB43497.1|_826:852 blasting: T0116-500-end gi|4929286|gb|AAD33934.1|AF144608_3_3:78 blasting: T0116-500-end gi|4929286|gb|AAD33934.1|AF144608_3_108:135 blasting: T0116-500-end gi|7387923|sp|Q9X105|MUS2_THEMA_254:523 blasting: T0116-500-end gi|6226653|sp|P74926|MUTS_THEMA_503:768 blasting: T0116-500-end gi|5305132|emb|CAB46193.1|_1:152 blasting: T0116-500-end gi|5305132|emb|CAB46193.1|_194:220 blasting: T0116-500-end gi|5725250|emb|CAB52406.1|_532:789 blasting: T0116-500-end gi|5814102|gb|AAD52100.1|_46:248 blasting: T0116-500-end gi|5815430|gb|AAD52669.1|AF178755_1_503:774 blasting: T0116-500-end gi|7512546|pir||T34526_39:241 blasting: T0116-500-end gi|6166561|sp|O13921|MSH1_SCHPO_669:928 blasting: T0116-500-end gi|6224917|gb|AAF06013.1|AF193018_1_849:1045 blasting: T0116-500-end gi|6227005|gb|AAF06041.1|AC009360_6_385:575 blasting: T0116-500-end gi|6227005|gb|AAF06041.1|AC009360_6_680:704 blasting: T0116-500-end gi|7473468|pir||A75332_257:519 blasting: T0116-500-end gi|6473463|dbj|BAA87137.1|_47:170 blasting: T0116-500-end gi|6518487|dbj|BAA87902.1|_1:100 blasting: T0116-500-end gi|7012942|gb|AAF35250.1|AF227632_1_351:563 blasting: T0116-500-end gi|7012942|gb|AAF35250.1|AF227632_1_301:322 blasting: T0116-500-end gi|7108493|gb|AAF36440.1|AF126491_1_526:794 blasting: T0116-500-end gi|7108493|gb|AAF36440.1|AF126491_1_827:852 blasting: T0116-500-end gi|7188744|gb|AAF37851.1|AF227730_1_9:298 blasting: T0116-500-end gi|7188744|gb|AAF37851.1|AF227730_1_304:345 blasting: T0116-500-end gi|7189834|gb|AAF38705.1|_557:808 blasting: T0116-500-end gi|7190210|gb|AAF39048.1|_542:799 blasting: T0116-500-end gi|7227419|gb|AAF42468.1|_520:801 blasting: T0116-500-end gi|7227419|gb|AAF42468.1|_829:859 blasting: T0116-500-end gi|7229490|gb|AAF42850.1|_526:795 blasting: T0116-500-end gi|7229490|gb|AAF42850.1|_828:853 blasting: T0116-500-end gi|7243293|dbj|BAA92687.1|_247:494 blasting: T0116-500-end gi|7243293|dbj|BAA92687.1|_536:562 blasting: T0116-500-end gi|7269404|emb|CAB81364.1|_755:893 blasting: T0116-500-end gi|7269404|emb|CAB81364.1|_890:1011 blasting: T0116-500-end gi|7287834|gb|AAF44872.1|AE003406_77_597:766 blasting: T0116-500-end gi|7287834|gb|AAF44872.1|AE003406_77_757:867 blasting: T0116-500-end gi|7294307|gb|AAF49656.1|_874:982 blasting: T0116-500-end gi|7298155|gb|AAF53392.1|_534:780 blasting: T0116-500-end gi|7328295|emb|CAB82462.1|_64:337 blasting: T0116-500-end gi|7328295|emb|CAB82462.1|_401:419 blasting: T0116-500-end gi|7331984|gb|AAF60672.1|_1539:1747 blasting: T0116-500-end gi|7379009|emb|CAB83555.1|_520:801 blasting: T0116-500-end gi|7379009|emb|CAB83555.1|_829:859 blasting: T0116-500-end gi|7404405|sp|P73769|MUTS_SYNY3_563:835 blasting: T0116-500-end gi|8052216|emb|CAB92351.1|_1:130 blasting: T0116-500-end gi|8473172|sp|P56883|MUTS_RHIME_574:863 blasting: T0116-500-end gi|8473172|sp|P56883|MUTS_RHIME_869:910 blasting: T0116-500-end gi|8515837|gb|AAF76200.1|_551:828 blasting: T0116-500-end gi|8809580|dbj|BAA97131.1|_423:634 blasting: T0116-500-end gi|9106782|gb|AAF84525.1|AE003995_7_579:860 blasting: T0116-500-end gi|9106782|gb|AAF84525.1|AE003995_7_871:902 blasting: T0116-500-end gi|9294039|dbj|BAB01996.1|_849:1045 blasting: T0116-500-end gi|9294568|dbj|BAB02831.1|_393:635 blasting: T0116-500-end gi|9294568|dbj|BAB02831.1|_682:699 blasting: T0116-500-end gi|9654962|gb|AAF93703.1|_547:803 blasting: T0116-500-end gi|9654962|gb|AAF93703.1|_833:861 blasting: T0116-500-end gi|9715729|emb|CAC01599.1|_111:281 blasting: T0116-500-end gi|9715729|emb|CAC01599.1|_322:421 blasting: T0116-500-end gi|9864530|gb|AAG00261.1|_681:889 blasting: T0116-500-end gi|9864530|gb|AAG00261.1|_627:702 blasting: make[1]: Leaving directory `/export/projects/compbio/experiments/casp4/t116/t116-500-end' make -k T0116-500-end.t2k.a2m.gz \ T0116-500-end.t2k.pa \ T0116-500-end.t2k.pa.html \ T0116-500-end.t2k.2d \ T0116-500-end.t2k.tree \ T0116-500-end.t2k_sorted.pa \ T0116-500-end.t2k.tree.ps \ T0116-500-end.t2k.tree-unroot.ps make[1]: Entering directory `/export/projects/compbio/experiments/casp4/t116/t116-500-end' /projects/compbio/experiments/models.97/scripts2k/target2k -out T0116-500-end.t2k \ -seed T0116-500-end.seq \ -all @@@@ chgrp protein /projects/compbio/tmp/target2k-beta-20006 @@@@ cp /export/projects/compbio/experiments/casp4/t116/t116-500-end/T0116-500-end.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-19791 /projects/compbio/experiments/models.97/scripts2k/blast-prefilter: E-value thresholds used are For /export/projects/compbio/tmp/target2k-beta-20006/prefilter_1.fa, 0.01 For /export/projects/compbio/tmp/target2k-beta-20006/prefilter_2.fa, 1 For /export/projects/compbio/tmp/target2k-beta-20006/prefilter_3.fa, 10 For /export/projects/compbio/tmp/target2k-beta-20006/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-20006/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-20006/m0.a2m (1 sequences, 311 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-19791 < /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, 311 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, 311 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.0643086816720257 -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: -142.185242 Database has 151 sequences with 120277 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 (148 sequences, 311 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=149 nll_thresh=-22.4115278685381 frac_id=0.854551641727273 @@@@ uniqueseq thinned_tmp -alignfile tmp_1-a.train.seq -a protein \ -percent_id 0.854551641727273 -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file tmp_1-a.train.seq (149 sequences, 311 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 27 (of 149) duplicate sequences with differing IDs Dropping 53 (of 149) sequences with > 85.5% id 69 sequences left after dropping 80 of 149 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 (69 sequences, 311 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 -855.20 33.00 -235.31 110.38 17 0 312 @@@@ 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.0643086816720257 \ -align_short 5 -mdEmax 0.0001 -constraints tmp_1-a.train.cst -constraints_out /export/projects/compbio/experiments/casp4/t116/t116-500-end/T0116-500-end.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: -260.338715 Database has 151 sequences with 120503 residues. @@@@ cp -f tmp_1-b.mult /export/projects/compbio/experiments/casp4/t116/t116-500-end/T0116-500-end.t2k_1.a2m # About to count sequences in prefilter_2.fa @@@@ modelfromalign tmp_2-a -alignfile /export/projects/compbio/experiments/casp4/t116/t116-500-end/T0116-500-end.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/t116/t116-500-end/T0116-500-end.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/t116/t116-500-end/T0116-500-end.t2k_1.a2m (151 sequences, 311 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.0643086816720257 -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/t116/t116-500-end/T0116-500-end.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: -231.619110 Database has 154 sequences with 122366 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 (151 sequences, 311 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=152 nll_thresh=-21.7183806877933 frac_id=0.845376512108611 @@@@ uniqueseq thinned_tmp -alignfile tmp_2-a.train.seq -a protein \ -percent_id 0.845376512108611 -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file tmp_2-a.train.seq (152 sequences, 311 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 24 (of 152) duplicate sequences with differing IDs Dropping 58 (of 152) sequences with > 84.5% id 70 sequences left after dropping 82 of 152 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 (70 sequences, 311 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/t116/t116-500-end/T0116-500-end.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 -854.79 -28.80 -230.33 94.03 3 0 312 @@@@ 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.0643086816720257 \ -align_short 5 -mdEmax 0.0002 -constraints tmp_2-a.train.cst -constraints_out /export/projects/compbio/experiments/casp4/t116/t116-500-end/T0116-500-end.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: -230.809341 Database has 153 sequences with 121158 residues. @@@@ cp -f tmp_2-b.mult /export/projects/compbio/experiments/casp4/t116/t116-500-end/T0116-500-end.t2k_2.a2m # About to count sequences in prefilter_3.fa @@@@ modelfromalign tmp_3-a -alignfile /export/projects/compbio/experiments/casp4/t116/t116-500-end/T0116-500-end.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/t116/t116-500-end/T0116-500-end.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/t116/t116-500-end/T0116-500-end.t2k_2.a2m (152 sequences, 311 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.0643086816720257 -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/t116/t116-500-end/T0116-500-end.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: -197.546463 Database has 173 sequences with 131264 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 (152 sequences, 311 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=153 nll_thresh=-20.1089427738805 frac_id=0.830586782911846 @@@@ uniqueseq thinned_tmp -alignfile tmp_3-a.train.seq -a protein \ -percent_id 0.830586782911846 -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file tmp_3-a.train.seq (153 sequences, 311 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 23 (of 153) duplicate sequences with differing IDs Dropping 61 (of 153) sequences with > 83.1% id 69 sequences left after dropping 84 of 153 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 (69 sequences, 311 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/t116/t116-500-end/T0116-500-end.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 -855.23 -28.73 -218.42 90.64 5 0 312 @@@@ 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.0643086816720257 \ -align_short 5 -mdEmax 0.001 -constraints tmp_3-a.train.cst -constraints_out /export/projects/compbio/experiments/casp4/t116/t116-500-end/T0116-500-end.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: -213.486679 Database has 156 sequences with 122655 residues. @@@@ cp -f tmp_3-b.mult /export/projects/compbio/experiments/casp4/t116/t116-500-end/T0116-500-end.t2k_3.a2m # About to count sequences in prefilter_4.fa @@@@ modelfromalign tmp_4-a -alignfile /export/projects/compbio/experiments/casp4/t116/t116-500-end/T0116-500-end.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/t116/t116-500-end/T0116-500-end.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/t116/t116-500-end/T0116-500-end.t2k_3.a2m (153 sequences, 311 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.0643086816720257 -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/t116/t116-500-end/T0116-500-end.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: -63.753021 Database has 574 sequences with 292957 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 (152 sequences, 311 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=153 nll_thresh=-18.4995048540528 frac_id=0.816114276630591 @@@@ uniqueseq thinned_tmp -alignfile tmp_4-a.train.seq -a protein \ -percent_id 0.816114276630591 -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file tmp_4-a.train.seq (153 sequences, 311 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 23 (of 153) duplicate sequences with differing IDs Dropping 62 (of 153) sequences with > 81.6% id 68 sequences left after dropping 85 of 153 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 (68 sequences, 311 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/t116/t116-500-end/T0116-500-end.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 -856.39 -33.41 -210.25 90.35 3 0 312 @@@@ 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.0643086816720257 \ -align_short 5 -mdEmax 0.005 -constraints tmp_4-a.train.cst -constraints_out /export/projects/compbio/experiments/casp4/t116/t116-500-end/T0116-500-end.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: -166.832108 Database has 193 sequences with 137405 residues. @@@@ cp -f tmp_4-b.mult /export/projects/compbio/experiments/casp4/t116/t116-500-end/T0116-500-end.t2k_4.a2m @@@@ modelfromalign tmp_5-a -alignfile /export/projects/compbio/experiments/casp4/t116/t116-500-end/T0116-500-end.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/t116/t116-500-end/T0116-500-end.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/t116/t116-500-end/T0116-500-end.t2k_4.a2m (153 sequences, 311 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.0643086816720257 -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/t116/t116-500-end/T0116-500-end.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: -64.279953 Database has 574 sequences with 292957 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 (150 sequences, 311 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/t116/t116-500-end/T0116-500-end.t2k.a2m @@@@ cp tmp_5-a.train.cst /export/projects/compbio/experiments/casp4/t116/t116-500-end/T0116-500-end.t2k.cst @@@@ rm -rf /projects/compbio/tmp/target2k-beta-20006 < /dev/null gzip -f T0116-500-end.t2k.a2m T0116-500-end.t2k_*.a2m prettyalign T0116-500-end.t2k.a2m.gz -m8 -i -n -L3333 > T0116-500-end.t2k.pa SAM: prettyalign v3.2 (July 31, 2000) compiled 07/31/00_16:52:38 gunzip -c T0116-500-end.t2k.a2m.gz > tmp.a2m a2m2html -a2m_in tmp.a2m -retrieve > T0116-500-end.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 T0116-500-end.t2k.a2m.gz >> tmp.script echo PrintPrediction T0116-500-end.t2k.2d 3670-4530-6947 >> tmp.script echo PrintPredictionFasta T0116-500-end.t2k.2d.seq >> tmp.script echo PrintRDB T0116-500-end.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 T0116-500-end.t2k.a2m.gz # Using SequenceWeight EntropyWeight(1.4, 10) # Using regularizer /projects/compbio/lib/recode3.20comp for sequence weight regularizer. # T0116-500-end.t2k.a2m with 151 sequences, total weight= 53.9247 avg weight= 0.357118 46 iterations # AdjustWeights couldn't save exactly 1.4 bits/position, saving 1.39847 bits. Regularizing alignment for T0116-500-end.t2k.a2m.gz # After reading T0116-500-end.t2k.a2m.gz, have 311 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 T0116-500-end.t2k.2d # command:# Network initialization done # Printing prediction in FASTA format to T0116-500-end.t2k.2d.seq # command:# Network initialization done # Printing prediction to T0116-500-end.t2k.2d.rdb # command:rm tmp.script gunzip -f T0116-500-end.t2k.a2m.gz /projects/compbio/usr/karplus/src/phytree/phytree -o -i -r flat T0116-500-end.t2k T0116-500-end.t2k.a2m /projects/compbio/lib/recode2.20comp Prior library /projects/compbio/lib/recode2.20comp read. Reading alignment file T0116-500-end.t2k.a2m (151 sequences) as A2M alignment. Alignment T0116-500-end.t2k.a2m read. Adding sequences:0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 One dot for every 10 merges: ............... Merges done. CPU = user 0: 4: 4 system 0: 0: 0 gzip -f T0116-500-end.t2k.a2m gzip -f T0116-500-end.t2k_sorted.a2m rm T0116-500-end.t2k.phytrace T0116-500-end.t2k.tree_weight prettyalign T0116-500-end.t2k_sorted.a2m.gz -m8 -i -n -L3333 > T0116-500-end.t2k_sorted.pa SAM: prettyalign v3.2 (July 31, 2000) compiled 07/31/00_16:52:38 /projects/compbio/usr/karplus/src/phytree/dg T0116-500-end.t2k.tree /projects/compbio/usr/karplus/src/phytree/dtree T0116-500-end.t2k.tree make[1]: Leaving directory `/export/projects/compbio/experiments/casp4/t116/t116-500-end' echo "WeightingRegularizer /projects/compbio/lib/recode3.20comp" > tmp.script echo "SequenceWeight EntropyWeight 0.5 10.0" >> tmp.script echo "ReadA2M T0116-500-end.t2k.a2m.gz" >> tmp.script echo "SavingsRegularizer /projects/compbio/lib/recode3.20comp" >> tmp.script echo "PrintSavings T0116-500-end.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 T0116-500-end.t2k.a2m.gz % using regularizer recode3.20comp % T0116-500-end.t2k.a2m.gz with 151 sequences, total weight= 3.92118 avg weight= 0.0259681 After reading T0116-500-end.t2k.a2m.gz, have 311 columns, totalling 1072.71 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 311 columns to T0116-500-end.t2k.saves # command:rm tmp.script if test -e T0116-500-end.t2k-logo.labels ;\ then makelogo -b 60 -t T0116-500-end.t2k -o T0116-500-end.t2k-logo.eps -l T0116-500-end.t2k-logo.labels T0116-500-end.t2k.saves ;\ else makelogo -b 60 -t T0116-500-end.t2k -o T0116-500-end.t2k-logo.eps -lf T0116-500-end.t2k.2d.seq T0116-500-end.t2k.saves ;\ fi /bin/sh: makelogo: execute permission denied make: *** [T0116-500-end.t2k-logo.eps] Error 1 rm T0116-500-end.t2k.saves make -k T0116-500-end-blast.rdb T0116-500-end-double-blast.rdb make[1]: Entering directory `/export/projects/compbio/experiments/casp4/t116/t116-500-end' make[1]: `T0116-500-end-blast.rdb' is up to date. make[1]: `T0116-500-end-double-blast.rdb' is up to date. make[1]: Leaving directory `/export/projects/compbio/experiments/casp4/t116/t116-500-end' make -k T0116-500-end.t2k.a2m.gz \ T0116-500-end.t2k.pa \ T0116-500-end.t2k.pa.html \ T0116-500-end.t2k.2d \ T0116-500-end.t2k.tree \ T0116-500-end.t2k_sorted.pa \ T0116-500-end.t2k.tree.ps \ T0116-500-end.t2k.tree-unroot.ps make[1]: Entering directory `/export/projects/compbio/experiments/casp4/t116/t116-500-end' make[1]: `T0116-500-end.t2k.a2m.gz' is up to date. make[1]: `T0116-500-end.t2k.pa' is up to date. make[1]: `T0116-500-end.t2k.pa.html' is up to date. make[1]: `T0116-500-end.t2k.2d' is up to date. make[1]: `T0116-500-end.t2k.tree' is up to date. make[1]: `T0116-500-end.t2k_sorted.pa' is up to date. make[1]: `T0116-500-end.t2k.tree.ps' is up to date. make[1]: `T0116-500-end.t2k.tree-unroot.ps' is up to date. make[1]: Leaving directory `/export/projects/compbio/experiments/casp4/t116/t116-500-end' /projects/compbio/experiments/models.97/scripts99/score-target -alignment T0116-500-end.t2k.a2m.gz -build w0.5 \ -db /projects/compbio/data/pdb/all-protein \ -Emax 40.0 \ -viterbi 0 -sw_score 2 \ -output T0116-500-end.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-4676 @@@@ mkdir /projects/compbio/tmp/score-target-4676/for-scorescript @@@@ /projects/compbio/experiments/models.97/scripts99/score-mods-vs-seqs.pl -root_tmp /projects/compbio/tmp/score-target-4676/for-scorescript -a2msfile /projects/compbio/tmp/score-target-4676/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-4676/score-target.op @@@@ chgrp protein /projects/compbio/tmp/score-target-4676/for-scorescript/split-tmp @@@@ chgrp protein /projects/compbio/tmp/score-target-4676/for-scorescript/split-tmp @@@@ /projects/compbio/bin/scripts/w0.5 /export/projects/compbio/experiments/casp4/t116/t116-500-end/T0116-500-end.t2k.a2m.gz /projects/compbio/tmp/score-target-4676/for-scorescript/split-tmp/beta-5076/T0116-500-end.t2k.mod Reading /projects/compbio/bin/scripts/sam-t99.conf @@@@ chgrp protein /projects/compbio/tmp/tmp-build-weighted-model-beta-3781 @@@@ uniqueseq /projects/compbio/tmp/tmp-build-weighted-model-beta-3781/tmp -alignfile /export/projects/compbio/experiments/casp4/t116/t116-500-end/T0116-500-end.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/t116/t116-500-end/T0116-500-end.t2k.a2m.gz (151 sequences, 311 columns) as A2M alignment. Writing sequence output to /projects/compbio/tmp/tmp-build-weighted-model-beta-3781/tmp.a2m. Dropping 22 (of 151) duplicate sequences with differing IDs Dropping 62 (of 151) sequences with > 80.0% id 67 sequences left after dropping 84 of 151 sequences. @@@@ modelfromalign /projects/compbio/tmp/tmp-build-weighted-model-beta-3781/tmp -alignfile /projects/compbio/tmp/tmp-build-weighted-model-beta-3781/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-3781/tmp.a2m (67 sequences, 311 columns) as A2M alignment. @@@@ cp -f /projects/compbio/tmp/tmp-build-weighted-model-beta-3781/tmp.mod /projects/compbio/tmp/score-target-4676/for-scorescript/split-tmp/beta-5076/T0116-500-end.t2k.mod @@@@ rm -rf /projects/compbio/tmp/tmp-build-weighted-model-beta-3781 80.0% id 67 sequences left after dropping 84 of 151 sequences. @@@@ modelfromalign /projects/compbio/tmp/tmp-build-weighted-model-beta-1332/tmp -alignfile /projects/compbio/tmp/tmp-build-weighted-model-beta-1332/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-1332/tmp.a2m (67 sequences, 311 columns) as A2M alignment. @@@@ cp -f /projects/compbio/tmp/tmp-build-weighted-model-beta-1332/tmp.mod T0116-500-end.t2k-w0.5.mod @@@@ rm -rf /projects/compbio/tmp/tmp-build-weighted-model-beta-1332