make -k T0099-blast.rdb T0099-double-blast.rdb make[1]: Entering directory `/auto/projects/compbio/experiments/casp4/t99' /projects/compbio/experiments/casp4/scripts/single-blast -q T0099.seq -interdb /projects/compbio/data/nrp/nr -db /projects/compbio/data/pdb/all-protein > T0099-blast.rdb starting T0099 with E=10 /projects/compbio/experiments/casp4/scripts/double-blast -q T0099.seq -interdb /projects/compbio/data/nrp/nr -db /projects/compbio/data/pdb/all-protein > T0099-double-blast.rdb starting T0099 with E=0.000050 T0099 T0099 blasting: T0099 gi|1174436|sp|P42686|SRK1_SPOLA_59:112 blasting: T0099 gi|1174439|sp|P42690|SRK4_SPOLA_59:112 blasting: T0099 gi|34295|emb|CAA31884.1|_67:119 blasting: T0099 gi|36808|emb|CAA28691.1|_67:119 blasting: T0099 gi|109465|pir||S20808_86:141 blasting: T0099 gi|280989|pir||A43807_70:125 blasting: T0099 gi|2117800|pir||I48845_67:119 blasting: T0099 gi|112234|pir||S24547_70:125 blasting: T0099 gi|57582|emb|CAA44218.1|_61:111 blasting: T0099 gi|66791|pir||TVFV60_86:141 blasting: T0099 gi|61504|emb|CAA24496.1|_370:425 blasting: T0099 gi|125712|sp|P25020|SRC_RSVH1_86:141 blasting: T0099 gi|628312|pir||S26420_86:141 blasting: T0099 gi|125714|sp|P00524|SRC_RSVSR_86:141 blasting: T0099 gi|462444|sp|Q05876|FYN_CHICK_87:142 blasting: T0099 gi|63363|emb|CAA31595.1|_94:149 blasting: T0099 gi|125869|sp|P09324|YES_CHICK_94:149 blasting: T0099 gi|462471|sp|Q02977|YRK_CHICK_86:141 blasting: T0099 gi|125372|sp|P27446|FYN_XIPHE_87:142 blasting: T0099 gi|125872|sp|P27447|YES_XIPHE_97:152 blasting: T0099 gi|125371|sp|P13406|FYN_XENLA_87:142 blasting: T0099 gi|125871|sp|P10936|YES_XENLA_90:145 blasting: T0099 gi|65305|emb|CAA45924.1|_80:135 blasting: T0099 gi|66784|pir||TVHUSR_87:142 blasting: T0099 gi|125717|sp|P17713|STK_HYDAT_64:116 blasting: T0099 gi|125474|sp|P06239|LCK_HUMAN_67:119 blasting: T0099 gi|125706|sp|P15054|SRC_AVIS2_86:141 blasting: T0099 gi|125713|sp|P00526|SRC_RSVP_86:141 blasting: T0099 gi|125709|sp|P14085|SRC_AVIST_86:141 blasting: T0099 gi|125708|sp|P14084|SRC_AVISS_86:141 blasting: T0099 gi|66797|pir||TVFVG9_86:141 blasting: T0099 gi|125133|sp|P00522|ABL_DROME_209:260 blasting: T0099 gi|71608|pir||MWAXIB_1097:1144 blasting: T0099 gi|77260|pir||S15582_86:141 blasting: T0099 gi|90368|pir||B24773_60:110 blasting: T0099 gi|90369|pir||C24773_61:111 blasting: T0099 gi|90370|pir||D24773_85:135 blasting: T0099 gi|125475|sp|P06240|LCK_MOUSE_67:119 blasting: T0099 gi|125703|sp|P13115|SRC1_XENLA_85:140 blasting: T0099 gi|125705|sp|P13116|SRC2_XENLA_85:140 blasting: T0099 gi|110658|pir||A39719_69:119 blasting: T0099 gi|125702|sp|P00528|SRC1_DROME_101:152 blasting: T0099 gi|125707|sp|P00525|SRC_AVISR_86:141 blasting: T0099 gi|134425|sp|P29355|SEM5_CAEEL_2:57 blasting: T0099 gi|1171093|sp|P19706|MYSB_ACACA_1097:1144 blasting: T0099 gi|179385|gb|AAA35593.1|_64:114 blasting: T0099 gi|179750|gb|AAA51896.1|_39:89 blasting: T0099 gi|180639|gb|AAA35697.1|_76:126 blasting: T0099 gi|2117805|pir||I53715_48:98 blasting: T0099 gi|191559|gb|AAA37135.1|_66:116 blasting: T0099 gi|191561|gb|AAA37136.1|_60:110 blasting: T0099 gi|191563|gb|AAA37137.1|_61:111 blasting: T0099 gi|191565|gb|AAA37138.1|_85:135 blasting: T0099 gi|198941|gb|AAA39471.1|_69:119 blasting: T0099 gi|198943|gb|AAA39472.1|_48:98 blasting: T0099 gi|210226|gb|AAA42573.1|_86:141 blasting: T0099 gi|210247|gb|AAA42581.1|_90:138 blasting: T0099 gi|210271|gb|AAA42585.1|_86:141 blasting: T0099 gi|212706|gb|AAA49078.1|_86:141 blasting: T0099 gi|212713|gb|AAA49081.1|_14:64 blasting: T0099 gi|221014|dbj|BAA01500.1|_86:141 blasting: T0099 gi|223452|prf||0806225A_370:425 blasting: T0099 gi|223627|prf||0903255A_332:387 blasting: T0099 gi|224525|prf||1107272A_115:164 blasting: T0099 gi|249851|gb|AAB96845.1|_86:141 blasting: T0099 gi|251793|gb|AAB22579.1|_59:112 blasting: T0099 gi|400155|sp|P31693|SRC_RSVPA_90:138 blasting: T0099 gi|7434383|pir||I67811_69:119 blasting: T0099 gi|7434384|pir||I67812_48:98 blasting: T0099 gi|2507209|sp|Q07014|LYN_RAT_69:119 blasting: T0099 gi|729629|sp|Q07883|GRB2_CHICK_2:53 blasting: T0099 gi|2144421|pir||TVHUHC_63:116 blasting: T0099 gi|321075|pir||A45501_90:145 blasting: T0099 gi|346653|pir||B45184_220:273 blasting: T0099 gi|347052|gb|AAC37877.1|_86:141 blasting: T0099 gi|420220|pir||JU0228_184:233 blasting: T0099 gi|423455|pir||A46243_166:218 blasting: T0099 gi|423549|pir||JN0471_222:275 blasting: T0099 gi|443263|pdb|1SHF|A_4:59 blasting: T0099 gi|2194099|pdb|1AEY|_10:58 blasting: T0099 gi|460966|gb|AAA18225.1|_67:119 blasting: T0099 gi|479367|pir||S33568_87:142 blasting: T0099 gi|498178|gb|AAC37549.1|_2:53 blasting: T0099 gi|1171094|sp|P42522|MYSC_DICDI_1126:1178 blasting: T0099 gi|553223|gb|AAA35679.1|_82:137 blasting: T0099 gi|554598|gb|AAA42570.1|_86:141 blasting: T0099 gi|554599|gb|AAA42571.1|_86:141 blasting: T0099 gi|555875|gb|AAB60448.1|_60:110 blasting: T0099 gi|555876|gb|AAB60449.1|_85:135 blasting: T0099 gi|555877|gb|AAB60450.1|_61:111 blasting: T0099 gi|555878|gb|AAB60451.1|_66:116 blasting: T0099 gi|2117817|pir||I39457_220:273 blasting: T0099 gi|576073|pdb|1CSK|A_13:69 blasting: T0099 gi|999679|pdb|1PRL|C_10:62 blasting: T0099 gi|2507603|sp|P35991|BTK_MOUSE_220:273 blasting: T0099 gi|625219|pir||TVHUSC_89:150 blasting: T0099 gi|625220|pir||TVCHS_86:141 blasting: T0099 gi|640062|pdb|1HSQ|_11:59 blasting: T0099 gi|1072767|pir||S52313_106:161 blasting: T0099 gi|1072768|pir||S52314_106:161 blasting: T0099 gi|732528|gb|AAC50116.1|_47:106 blasting: T0099 gi|2117803|pir||I48608_54:113 blasting: T0099 gi|775208|gb|AAC50287.1|_67:119 blasting: T0099 gi|802051|gb|AAB33113.1|_87:142 blasting: T0099 gi|901820|gb|AAA70194.1|_86:141 blasting: T0099 gi|914204|gb|AAB33265.1|_64:117 blasting: T0099 gi|914957|dbj|BAA08645.1|_2:53 blasting: T0099 gi|914961|dbj|BAA08646.1|_2:53 blasting: T0099 gi|939625|gb|AAC52725.1|_48:107 blasting: T0099 gi|6686295|sp|Q02779|M3KA_HUMAN_23:77 blasting: T0099 gi|999906|pdb|1GBR|A_11:62 blasting: T0099 gi|7511162|pir||T27877_4:55 blasting: T0099 gi|7511162|pir||T27877_217:268 blasting: T0099 gi|1064965|pdb|1CKA|A_2:54 blasting: T0099 gi|1071850|pir||A39939_65:115 blasting: T0099 gi|1170188|sp|P08631|HCK_HUMAN_84:137 blasting: T0099 gi|1170189|sp|P08103|HCK_MOUSE_82:132 blasting: T0099 gi|1170731|sp|P42683|LCK_CHICK_66:116 blasting: T0099 gi|1174631|sp|P24604|TEC_MOUSE_184:233 blasting: T0099 gi|1311313|pdb|1LCK|A_16:68 blasting: T0099 gi|1311410|pdb|1ABO|A_7:57 blasting: T0099 gi|1351818|sp|P00527|YES_AVISY_86:141 blasting: T0099 gi|1362664|pir||B49114_89:144 blasting: T0099 gi|1362665|pir||A49114_92:147 blasting: T0099 gi|1363239|pir||A57152_25:78 blasting: T0099 gi|1373390|gb|AAB05596.1|_192:243 blasting: T0099 gi|1375041|dbj|BAA12862.1|_2:53 blasting: T0099 gi|1942473|pdb|1NLO|C_10:62 blasting: T0099 gi|2498936|sp|Q15811|ITSN_HUMAN_459:509 blasting: T0099 gi|1438935|gb|AAC50593.1|_41:90 blasting: T0099 gi|1438935|gb|AAC50593.1|_187:237 blasting: T0099 gi|3025159|sp|Q28923|YES_CANFA_92:147 blasting: T0099 gi|1536790|dbj|BAA07705.1|_68:120 blasting: T0099 gi|2499805|sp|P55345|ANM2_HUMAN_34:84 blasting: T0099 gi|6175046|sp|P00523|SRC_CHICK_86:141 blasting: T0099 gi|1710263|gb|AAB50221.1|_34:84 blasting: T0099 gi|1816662|gb|AAC60143.1|_197:247 blasting: T0099 gi|1816662|gb|AAC60143.1|_9:57 blasting: T0099 gi|3913784|sp|P87379|GRB2_XENLA_2:53 blasting: T0099 gi|1913783|dbj|BAA18950.1|_110:158 blasting: T0099 gi|1942719|pdb|1EFN|C_3:58 blasting: T0099 gi|2982073|pdb|1AZG|B_6:61 blasting: T0099 gi|1943539|pdb|1FYN|A_7:62 blasting: T0099 gi|2105002|gb|AAB71345.1|_48:98 blasting: T0099 gi|2114076|dbj|BAA20078.1|_45:95 blasting: T0099 gi|2114412|gb|AAC47535.1|_1000:1046 blasting: T0099 gi|2194103|pdb|1AD5|A_6:59 blasting: T0099 gi|5822499|pdb|2SRC|_5:60 blasting: T0099 gi|7513197|pir||T08855_154:205 blasting: T0099 gi|2507208|sp|P25911|LYN_MOUSE_69:119 blasting: T0099 gi|2554681|pdb|2ABL|_11:61 blasting: T0099 gi|7512152|pir||T09194_990:1039 blasting: T0099 gi|7512152|pir||T09194_1143:1193 blasting: T0099 gi|2695869|emb|CAB05103.1|_2:51 blasting: T0099 gi|2723311|dbj|BAA24063.1|_165:214 blasting: T0099 gi|2723313|dbj|BAA24064.1|_348:397 blasting: T0099 gi|2738777|gb|AAB94550.1|_12:71 blasting: T0099 gi|2780916|pdb|2PTK|_6:61 blasting: T0099 gi|2827464|gb|AAB99858.1|_150:199 blasting: T0099 gi|2914378|pdb|1AWO|_8:58 blasting: T0099 gi|2921549|gb|AAD11958.1|_344:394 blasting: T0099 gi|2921551|gb|AAD11959.1|_363:413 blasting: T0099 gi|2967840|gb|AAC05835.1|_13:69 blasting: T0099 gi|2982012|pdb|1A0N|B_8:63 blasting: T0099 gi|7511852|pir||T13053_781:833 blasting: T0099 gi|7511851|pir||T13055_781:833 blasting: T0099 gi|3003004|gb|AAC08990.1|_86:141 blasting: T0099 gi|3114379|pdb|1AWW|_12:65 blasting: T0099 gi|3114543|pdb|1AVZ|C_3:55 blasting: T0099 gi|3318781|pdb|4HCK|_13:66 blasting: T0099 gi|3411274|gb|AAC31209.1|_57:107 blasting: T0099 gi|3560565|gb|AAC35011.1|_26:79 blasting: T0099 gi|7494125|pir||T32734_1166:1212 blasting: T0099 gi|3721908|dbj|BAA33741.1|_69:121 blasting: T0099 gi|3859853|gb|AAC78610.1|_1159:1209 blasting: T0099 gi|3859853|gb|AAC78610.1|_1006:1055 blasting: T0099 gi|3891430|pdb|1BU1|A_4:57 blasting: T0099 gi|3930217|gb|AAC80284.1|_203:252 blasting: T0099 gi|3930217|gb|AAC80284.1|_8:56 blasting: T0099 gi|4033555|emb|CAA10376.1|_164:214 blasting: T0099 gi|4033557|emb|CAA10377.1|_139:189 blasting: T0099 gi|4100623|gb|AAD00899.1|_164:213 blasting: T0099 gi|4100623|gb|AAD00899.1|_310:360 blasting: T0099 gi|4138839|gb|AAD03595.1|_669:717 blasting: T0099 gi|4139932|pdb|1BBZ|C_3:53 blasting: T0099 gi|4139934|pdb|1BBZ|E_3:53 blasting: T0099 gi|4151946|gb|AAD04633.1|_191:241 blasting: T0099 gi|4249651|gb|AAD13752.1|_197:247 blasting: T0099 gi|4249651|gb|AAD13752.1|_8:56 blasting: T0099 gi|4249651|gb|AAD13752.1|_115:160 blasting: T0099 gi|4378885|gb|AAD19746.1|_1152:1202 blasting: T0099 gi|4378885|gb|AAD19746.1|_999:1048 blasting: T0099 gi|4558122|pdb|1BK2|_5:53 blasting: T0099 gi|8134704|sp|Q9WUD9|SRC_RAT_89:141 blasting: T0099 gi|4731289|gb|AAD28428.1|AF112976_1_112:161 blasting: T0099 gi|4808823|gb|AAD29952.1|AF114487_1_1159:1209 blasting: T0099 gi|4808823|gb|AAD29952.1|AF114487_1_1006:1055 blasting: T0099 gi|4808825|gb|AAD29953.1|AF114488_1_1159:1209 blasting: T0099 gi|4808825|gb|AAD29953.1|AF114488_1_1006:1055 blasting: T0099 gi|4835853|gb|AAD30271.1|AF127798_1_1156:1206 blasting: T0099 gi|4835853|gb|AAD30271.1|AF127798_1_1003:1052 blasting: T0099 gi|4838526|gb|AAD31026.1|AF132672_1_1085:1135 blasting: T0099 gi|5107454|pdb|1AZE|A_2:53 blasting: T0099 gi|5107624|pdb|1PWT|_9:57 blasting: T0099 gi|5262742|emb|CAB45687.1|_160:209 blasting: T0099 gi|7492740|pir||T39427_1111:1160 blasting: T0099 gi|5453029|gb|AAD43402.1|_184:233 blasting: T0099 gi|5453032|gb|AAD43405.1|_184:233 blasting: T0099 gi|5542349|pdb|1QCF|A_12:65 blasting: T0099 gi|5578991|emb|CAB51395.1|_282:335 blasting: T0099 gi|5731279|gb|AAD48847.1|AF169620_1_46:96 blasting: T0099 gi|5731281|gb|AAD48848.1|AF169621_1_187:237 blasting: T0099 gi|5731281|gb|AAD48848.1|AF169621_1_34:83 blasting: T0099 gi|5823552|gb|AAD53183.1|AF180522_1_126:176 blasting: T0099 gi|5853151|gb|AAD54337.1|AF177536_1_44:95 blasting: T0099 gi|6006522|emb|CAA36457.2|_82:137 blasting: T0099 gi|6103277|emb|CAB59279.1|_2:53 blasting: T0099 gi|6651025|gb|AAF22134.1|AF127180_1_160:210 blasting: T0099 gi|6651163|gb|AAF22211.1|AF139492_1_435:488 blasting: T0099 gi|6651165|gb|AAF22212.1|AF139493_1_394:447 blasting: T0099 gi|6651167|gb|AAF22213.1|AF139494_1_433:486 blasting: T0099 gi|6651169|gb|AAF22214.1|AF139495_1_392:445 blasting: T0099 gi|6729971|pdb|1B07|A_4:56 blasting: T0099 gi|6730486|pdb|1QLY|A_5:58 blasting: T0099 gi|6969624|gb|AAF33845.1|AF226866_1_122:171 blasting: T0099 gi|6984209|gb|AAF34794.1|AF228313_1_54:106 blasting: T0099 gi|7160701|emb|CAB04427.2|_62:114 blasting: T0099 gi|7271012|emb|CAB77642.1|_315:368 blasting: T0099 gi|7292519|gb|AAF47922.1|_101:152 blasting: T0099 gi|7294076|gb|AAF49431.1|_192:243 blasting: T0099 gi|7294985|gb|AAF50313.1|_549:598 blasting: T0099 gi|7296156|gb|AAF51449.1|_330:381 blasting: T0099 gi|7298751|gb|AAF53962.1|_784:836 blasting: T0099 gi|7302199|gb|AAF57295.1|_68:120 blasting: T0099 gi|7303067|gb|AAF58135.1|_168:216 blasting: T0099 gi|7329078|gb|AAF59904.1|AF182199_1_984:1033 blasting: T0099 gi|7329078|gb|AAF59904.1|AF182199_1_1130:1180 blasting: T0099 gi|7331201|gb|AAF60319.1|AF236119_1_2:51 blasting: T0099 gi|7503796|pir||T22405_62:114 blasting: T0099 gi|7573534|emb|CAB87592.1|_47:106 blasting: T0099 gi|7578848|gb|AAF64151.1|AF232691_1_117:169 blasting: T0099 gi|7688244|emb|CAB89868.1|_300:349 blasting: make[1]: Leaving directory `/auto/projects/compbio/experiments/casp4/t99' make -k T0099.t2k.a2m.gz \ T0099.t2k.pa \ T0099.t2k.pa.html \ T0099.t2k.2d \ T0099.t2k.tree \ T0099.t2k_sorted.pa \ T0099.t2k.tree.ps \ T0099.t2k.tree-unroot.ps make[1]: Entering directory `/auto/projects/compbio/experiments/casp4/t99' /projects/compbio/experiments/models.97/scripts2k/target2k -out T0099.t2k \ -seed T0099.seq \ -all @@@@ chgrp protein /projects/compbio/tmp/target2k-titanic-29378 @@@@ cp /projects/compbio/experiments/casp4/t99/T0099.seq m0.a2m @@@@ checkseq foo -alphabet protein -db m0.a2m > init.check SAM: checkseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:44 # About to count sequences in /projects/compbio/data/nrp/nr /projects/compbio/experiments/models.97/scripts2k/target2k: E-value thresholds used are prefilter: 0.01, actual: 0.0001 prefilter: 1, actual: 0.0002 prefilter: 10, actual: 0.001 prefilter: 400, actual: 0.005 /projects/compbio/experiments/models.97/scripts2k/target2k: aweight_bits are 0.8, 0.7, 0.6, 0.5 /projects/compbio/experiments/models.97/scripts2k/target2k: db_size= 509169 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-titanic-29452 /projects/compbio/experiments/models.97/scripts2k/blast-prefilter: E-value thresholds used are For /projects/compbio/tmp/target2k-titanic-29378/prefilter_1.fa, 0.01 For /projects/compbio/tmp/target2k-titanic-29378/prefilter_2.fa, 1 For /projects/compbio/tmp/target2k-titanic-29378/prefilter_3.fa, 10 For /projects/compbio/tmp/target2k-titanic-29378/prefilter_4.fa, 400 /projects/compbio/experiments/models.97/scripts2k/blast-prefilter: db=/projects/compbio/data/nrp/nr @@@@ uniqueseq init -alignfile /projects/compbio/tmp/target2k-titanic-29378/m0.a2m -a protein -a2mdots 0 SAM: uniqueseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:50 Reading alignment file /projects/compbio/tmp/target2k-titanic-29378/m0.a2m (1 sequences, 56 columns) as A2M alignment. Writing sequence output to init.a2m. No sequences have been droppd. @@@@ /projects/compbio/bin/wu-blastp /projects/compbio/data/nrp/nr 1.fasta E=400 V=20000 B=0 -gi > 1.fasta-blast.out Searched 509169 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-titanic-29452 < /dev/null @@@@ modelfromalign unused -alignfile m0.a2m \ -constraints_from_align 1 -constraints_out m0.cst SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file m0.a2m (1 sequences, 56 columns) as A2M alignment. # About to count sequences in m0.a2m # About to count sequences in prefilter_1.fa @@@@ modelfromalign tmp_1-a -alignfile m0.a2m \ -insert /projects/compbio/lib/gap1.5.regularizer \ -aweight_bits 0.8\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 -constraints m0.cst -constraints_out tmp_1-a.cst Reading parameter file /projects/compbio/lib/gap1.5.regularizer /projects/compbio/lib/gap1.5.regularizer(1): Reading REGULARIZER: gap1.5. Weak regularizer, gaps average 1.5, matches average 10 SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file m0.a2m (1 sequences, 56 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 509169 \ -simple_threshold -26.3508756357805 -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.357142857142857 -db prefilter_1.fa -constraints tmp_1-a.cst -constraints_out tmp_1-a-md.cst Reading parameter file tmp_1-a.mod tmp_1-a.mod(23): Reading MODEL -- Model from alignment file m0.a2m SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -32.422516 Database has 513 sequences with 298565 residues. @@@@ sortseq tmp_1-a-sort -alignfile tmp_1-a.mult \ -NLLfile tmp_1-a.mstat -a2mdots 0 SAM: sortseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:38 Reading alignment file tmp_1-a.mult (183 sequences, 56 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=184 nll_thresh=-22.3508756357805 frac_id=0.816906638443881 @@@@ uniqueseq thinned_tmp -alignfile tmp_1-a.train.seq -a protein \ -percent_id 0.816906638443881 -a2mdots 0 SAM: uniqueseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:50 Reading alignment file tmp_1-a.train.seq (184 sequences, 56 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 93 (of 184) duplicate sequences with differing IDs Dropping 50 (of 184) sequences with > 81.7% id 41 sequences left after dropping 143 of 184 sequences. @@@@ modelfromalign junk_for_weights \ -alignfile tmp_1-a.train.seq \ -aweight_bits 0.8\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file tmp_1-a.train.seq (41 sequences, 56 columns) as A2M alignment. @@@@ buildmodel tmp_1-b -insert tmp_1-a.mod\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -modellength 0 \ -internal_weight 0 \ -anneal_length 2 \ -many_files 1 \ -Nseq 20000\ -nsurgery 0 \ -fracinsert 0.2 \ -cutinsert 0.25 -cutmatch 0.2 -initial_noise 0.0001 -anneal_noise 0.000001 \ -internal_weight 0 \ -train tmp_1-a.train.seq \ -nmodels 1\ -insert /projects/compbio/lib/gap1.5.regularizer -sequence_weights tmp_1-a.w -constraints tmp_1-a.train.cst Reading parameter file tmp_1-a.mod tmp_1-a.mod(23): Reading MODEL -- Model from alignment file m0.a2m Command line option(48): internal_weight has already been set. Using 0. Reading parameter file /projects/compbio/lib/gap1.5.regularizer /projects/compbio/lib/gap1.5.regularizer(1): Reading REGULARIZER: gap1.5. Weak regularizer, gaps average 1.5, matches average 10 SAM: buildmodel v3.1b (February 24, 1999) compiled 04/18/00_11:43:59 -154.88 -44.17 -58.52 16.42 7 0 57 @@@@ 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 509169 \ -db m0.a2m -db tmp_1-a.sel \ -select_score 4 \ -select_seq 4 \ -simple_threshold -26.3508756357805 -Emax 0.0001 \ -select_mdalign 4 -select_md 4 \ -adpstyle 1 \ -Motifcutoff 0.357142857142857 \ -align_short 5 -mdEmax 0.0001 -constraints tmp_1-a.train.cst -constraints_out /projects/compbio/experiments/casp4/t99/T0099.t2k_1.cst Reading parameter file tmp_1-b.mod tmp_1-b.mod(20): Reading MODEL -- Final model for run tmp_1-b SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -59.023819 Database has 228 sequences with 116650 residues. @@@@ cp -f tmp_1-b.mult /projects/compbio/experiments/casp4/t99/T0099.t2k_1.a2m # About to count sequences in prefilter_2.fa @@@@ modelfromalign tmp_2-a -alignfile /projects/compbio/experiments/casp4/t99/T0099.t2k_1.a2m \ -insert /projects/compbio/lib/stiff-gap5.regularizer \ -aweight_bits 0.7\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 -constraints /projects/compbio/experiments/casp4/t99/T0099.t2k_1.cst -constraints_out tmp_2-a.cst Reading parameter file /projects/compbio/lib/stiff-gap5.regularizer /projects/compbio/lib/stiff-gap5.regularizer(1): Reading REGULARIZER: gap5. stiff regularizer, gaps averge 5, matches average 10 SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file /projects/compbio/experiments/casp4/t99/T0099.t2k_1.a2m (295 sequences, 56 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 509169 \ -simple_threshold -25.6577284550241 -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.357142857142857 -db prefilter_2.fa -constraints tmp_2-a.cst -constraints_out tmp_2-a-md.cst Reading parameter file tmp_2-a.mod tmp_2-a.mod(23): Reading MODEL -- Model from alignment file /projects/compbio/experiments/casp4/t99/T0099.t2k_1.a2m SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -44.341564 Database has 728 sequences with 466873 residues. @@@@ sortseq tmp_2-a-sort -alignfile tmp_2-a.mult \ -NLLfile tmp_2-a.mstat -a2mdots 0 SAM: sortseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:38 Reading alignment file tmp_2-a.mult (891 sequences, 56 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=892 nll_thresh=-21.6577284550241 frac_id=0.648322769864098 @@@@ uniqueseq thinned_tmp -alignfile tmp_2-a.train.seq -a protein \ -percent_id 0.648322769864098 -a2mdots 0 SAM: uniqueseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:50 Reading alignment file tmp_2-a.train.seq (892 sequences, 56 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 402 (of 892) duplicate sequences with differing IDs Dropping 276 (of 892) sequences with > 64.8% id 214 sequences left after dropping 678 of 892 sequences. @@@@ modelfromalign junk_for_weights \ -alignfile tmp_2-a.train.seq \ -aweight_bits 0.7\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file tmp_2-a.train.seq (214 sequences, 56 columns) as A2M alignment. @@@@ buildmodel tmp_2-b -insert tmp_2-a.mod\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -modellength 0 \ -internal_weight 0 \ -anneal_length 2 \ -many_files 1 \ -Nseq 20000\ -nsurgery 0 \ -fracinsert 0.2 \ -cutinsert 0.25 -cutmatch 0.2 -initial_noise 0.0001 -anneal_noise 0.000001 \ -internal_weight 0 \ -train tmp_2-a.train.seq \ -nmodels 1\ -insert /projects/compbio/lib/stiff-gap5.regularizer -sequence_weights tmp_2-a.w -constraints tmp_2-a.train.cst Reading parameter file tmp_2-a.mod tmp_2-a.mod(23): Reading MODEL -- Model from alignment file /projects/compbio/experiments/casp4/t99/T0099.t2k_1.a2m Command line option(48): internal_weight has already been set. Using 0. Reading parameter file /projects/compbio/lib/stiff-gap5.regularizer /projects/compbio/lib/stiff-gap5.regularizer(1): Reading REGULARIZER: gap5. stiff regularizer, gaps averge 5, matches average 10 SAM: buildmodel v3.1b (February 24, 1999) compiled 04/18/00_11:43:59 -154.77 -32.43 -42.85 9.14 4 0 57 @@@@ 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 509169 \ -db m0.a2m -db tmp_2-a.sel \ -select_score 4 \ -select_seq 4 \ -simple_threshold -25.6577284550241 -Emax 0.0002 \ -select_mdalign 4 -select_md 4 \ -adpstyle 1 \ -Motifcutoff 0.357142857142857 \ -align_short 5 -mdEmax 0.0002 -constraints tmp_2-a.train.cst -constraints_out /projects/compbio/experiments/casp4/t99/T0099.t2k_2.cst Reading parameter file tmp_2-b.mod tmp_2-b.mod(20): Reading MODEL -- Final model for run tmp_2-b SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -45.917561 Database has 692 sequences with 428657 residues. @@@@ cp -f tmp_2-b.mult /projects/compbio/experiments/casp4/t99/T0099.t2k_2.a2m # About to count sequences in prefilter_3.fa @@@@ modelfromalign tmp_3-a -alignfile /projects/compbio/experiments/casp4/t99/T0099.t2k_2.a2m \ -insert /projects/compbio/lib/stiff-gap5.regularizer \ -aweight_bits 0.6\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 -constraints /projects/compbio/experiments/casp4/t99/T0099.t2k_2.cst -constraints_out tmp_3-a.cst Reading parameter file /projects/compbio/lib/stiff-gap5.regularizer /projects/compbio/lib/stiff-gap5.regularizer(1): Reading REGULARIZER: gap5. stiff regularizer, gaps averge 5, matches average 10 SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file /projects/compbio/experiments/casp4/t99/T0099.t2k_2.a2m (912 sequences, 56 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 509169 \ -simple_threshold -24.0482905410188 -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.357142857142857 -db prefilter_3.fa -constraints tmp_3-a.cst -constraints_out tmp_3-a-md.cst Reading parameter file tmp_3-a.mod tmp_3-a.mod(23): Reading MODEL -- Model from alignment file /projects/compbio/experiments/casp4/t99/T0099.t2k_2.a2m SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -37.048569 Database has 954 sequences with 578210 residues. @@@@ sortseq tmp_3-a-sort -alignfile tmp_3-a.mult \ -NLLfile tmp_3-a.mstat -a2mdots 0 SAM: sortseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:38 Reading alignment file tmp_3-a.mult (984 sequences, 56 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=985 nll_thresh=-20.0482905410188 frac_id=0.64056723767318 @@@@ uniqueseq thinned_tmp -alignfile tmp_3-a.train.seq -a protein \ -percent_id 0.64056723767318 -a2mdots 0 SAM: uniqueseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:50 Reading alignment file tmp_3-a.train.seq (985 sequences, 56 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 431 (of 985) duplicate sequences with differing IDs Dropping 302 (of 985) sequences with > 64.1% id 252 sequences left after dropping 733 of 985 sequences. @@@@ modelfromalign junk_for_weights \ -alignfile tmp_3-a.train.seq \ -aweight_bits 0.6\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file tmp_3-a.train.seq (252 sequences, 56 columns) as A2M alignment. @@@@ buildmodel tmp_3-b -insert tmp_3-a.mod\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -modellength 0 \ -internal_weight 0 \ -anneal_length 2 \ -many_files 1 \ -Nseq 20000\ -nsurgery 0 \ -fracinsert 0.2 \ -cutinsert 0.25 -cutmatch 0.2 -initial_noise 0.0001 -anneal_noise 0.000001 \ -internal_weight 0 \ -train tmp_3-a.train.seq \ -nmodels 1\ -insert /projects/compbio/lib/stiff-gap5.regularizer -sequence_weights tmp_3-a.w -constraints tmp_3-a.train.cst Reading parameter file tmp_3-a.mod tmp_3-a.mod(23): Reading MODEL -- Model from alignment file /projects/compbio/experiments/casp4/t99/T0099.t2k_2.a2m Command line option(48): internal_weight has already been set. Using 0. Reading parameter file /projects/compbio/lib/stiff-gap5.regularizer /projects/compbio/lib/stiff-gap5.regularizer(1): Reading REGULARIZER: gap5. stiff regularizer, gaps averge 5, matches average 10 SAM: buildmodel v3.1b (February 24, 1999) compiled 04/18/00_11:43:59 -154.96 -24.63 -38.65 9.05 2 0 57 @@@@ 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 509169 \ -db m0.a2m -db tmp_3-a.sel \ -select_score 4 \ -select_seq 4 \ -simple_threshold -24.0482905410188 -Emax 0.001 \ -select_mdalign 4 -select_md 4 \ -adpstyle 1 \ -Motifcutoff 0.357142857142857 \ -align_short 5 -mdEmax 0.001 -constraints tmp_3-a.train.cst -constraints_out /projects/compbio/experiments/casp4/t99/T0099.t2k_3.cst Reading parameter file tmp_3-b.mod tmp_3-b.mod(20): Reading MODEL -- Final model for run tmp_3-b SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -42.057308 Database has 761 sequences with 501680 residues. @@@@ cp -f tmp_3-b.mult /projects/compbio/experiments/casp4/t99/T0099.t2k_3.a2m # About to count sequences in prefilter_4.fa @@@@ modelfromalign tmp_4-a -alignfile /projects/compbio/experiments/casp4/t99/T0099.t2k_3.a2m \ -insert /projects/compbio/lib/fssp-trained.regularizer \ -aweight_bits 0.5\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 -constraints /projects/compbio/experiments/casp4/t99/T0099.t2k_3.cst -constraints_out tmp_4-a.cst Reading parameter file /projects/compbio/lib/fssp-trained.regularizer /projects/compbio/lib/fssp-trained.regularizer(1): Reading REGULARIZER: transitions from FSSP and pseudocounts for diverse weighted HSSP SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file /projects/compbio/experiments/casp4/t99/T0099.t2k_3.a2m (987 sequences, 56 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 509169 \ -simple_threshold -22.4388526207288 -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.357142857142857 -db prefilter_4.fa -constraints tmp_4-a.cst -constraints_out tmp_4-a-md.cst Reading parameter file tmp_4-a.mod tmp_4-a.mod(23): Reading MODEL -- Model from alignment file /projects/compbio/experiments/casp4/t99/T0099.t2k_3.a2m SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) ..... Average NLL-Simple NULL score: -12.988875 Database has 4827 sequences with 1895713 residues. @@@@ sortseq tmp_4-a-sort -alignfile tmp_4-a.mult \ -NLLfile tmp_4-a.mstat -a2mdots 0 SAM: sortseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:38 Reading alignment file tmp_4-a.mult (1023 sequences, 56 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=1024 nll_thresh=-18.4388526207288 frac_id=0.635916390055988 @@@@ uniqueseq thinned_tmp -alignfile tmp_4-a.train.seq -a protein \ -percent_id 0.635916390055988 -a2mdots 0 SAM: uniqueseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:50 Reading alignment file tmp_4-a.train.seq (1024 sequences, 56 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 450 (of 1024) duplicate sequences with differing IDs Dropping 310 (of 1024) sequences with > 63.6% id 264 sequences left after dropping 760 of 1024 sequences. @@@@ modelfromalign junk_for_weights \ -alignfile tmp_4-a.train.seq \ -aweight_bits 0.5\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file tmp_4-a.train.seq (264 sequences, 56 columns) as A2M alignment. @@@@ buildmodel tmp_4-b -insert tmp_4-a.mod\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1\ -modellength 0 \ -internal_weight 0 \ -anneal_length 2 \ -many_files 1 \ -Nseq 20000\ -nsurgery 0 \ -fracinsert 0.2 \ -cutinsert 0.25 -cutmatch 0.2 -initial_noise 0.0001 -anneal_noise 0.000001 \ -internal_weight 0 \ -train tmp_4-a.train.seq \ -nmodels 1\ -insert /projects/compbio/lib/fssp-trained.regularizer -sequence_weights tmp_4-a.w -constraints tmp_4-a.train.cst Reading parameter file tmp_4-a.mod tmp_4-a.mod(23): Reading MODEL -- Model from alignment file /projects/compbio/experiments/casp4/t99/T0099.t2k_3.a2m Command line option(48): internal_weight has already been set. Using 0. Reading parameter file /projects/compbio/lib/fssp-trained.regularizer /projects/compbio/lib/fssp-trained.regularizer(1): Reading REGULARIZER: transitions from FSSP and pseudocounts for diverse weighted HSSP SAM: buildmodel v3.1b (February 24, 1999) compiled 04/18/00_11:43:59 -155.36 -24.35 -38.18 9.10 4 0 57 @@@@ 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 509169 \ -db m0.a2m -db tmp_4-a.sel \ -select_score 4 \ -select_seq 4 \ -simple_threshold -22.4388526207288 -Emax 0.005 \ -select_mdalign 4 -select_md 4 \ -adpstyle 1 \ -Motifcutoff 0.357142857142857 \ -align_short 5 -mdEmax 0.005 -constraints tmp_4-a.train.cst -constraints_out /projects/compbio/experiments/casp4/t99/T0099.t2k_4.cst Reading parameter file tmp_4-b.mod tmp_4-b.mod(20): Reading MODEL -- Final model for run tmp_4-b SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) . Average NLL-Simple NULL score: -40.000755 Database has 813 sequences with 542870 residues. @@@@ cp -f tmp_4-b.mult /projects/compbio/experiments/casp4/t99/T0099.t2k_4.a2m @@@@ modelfromalign tmp_5-a -alignfile /projects/compbio/experiments/casp4/t99/T0099.t2k_4.a2m \ -insert /projects/compbio/lib/fssp-trained.regularizer \ -aweight_bits 0.5\ -fimtrans -1\ -fimstrength 1\ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -binary_output 1\ -prior_library /projects/compbio/lib/recode3.20comp \ -a2mdots 0 \ -a protein \ -sw 2 -jump_in_prob 0.2 -jump_out_prob 1 \ -aweight_method 1\ -aweight_exponent 10 -constraints /projects/compbio/experiments/casp4/t99/T0099.t2k_4.cst -constraints_out tmp_5-a.cst Reading parameter file /projects/compbio/lib/fssp-trained.regularizer /projects/compbio/lib/fssp-trained.regularizer(1): Reading REGULARIZER: transitions from FSSP and pseudocounts for diverse weighted HSSP SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file /projects/compbio/experiments/casp4/t99/T0099.t2k_4.a2m (1032 sequences, 56 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 509169 \ -simple_threshold -22.4388526207288 -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.357142857142857 -db prefilter_4.fa -constraints tmp_5-a.cst -constraints_out tmp_5-a-md.cst Reading parameter file tmp_5-a.mod tmp_5-a.mod(23): Reading MODEL -- Model from alignment file /projects/compbio/experiments/casp4/t99/T0099.t2k_4.a2m SAM: hmmscore v3.1b (February 24, 1999) compiled 05/16/00_17:09:35 Calculating Distances (. = 1000 sequences) ..... Average NLL-Simple NULL score: -13.039355 Database has 4827 sequences with 1895713 residues. @@@@ sortseq tmp_5-a-sort -alignfile tmp_5-a.mult \ -NLLfile tmp_5-a.mstat -a2mdots 0 SAM: sortseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:38 Reading alignment file tmp_5-a.mult (870 sequences, 56 columns) as A2M alignment. Reading scores from file tmp_5-a.mstat Writing sequence output to tmp_5-a-sort.a2m. @@@@ cat m0.a2m tmp_5-a-sort.a2m > tmp_5-a.train.seq @@@@ cat m0.cst tmp_5-a-md.cst > tmp_5-a.train.cst @@@@ cp tmp_5-a.train.seq /projects/compbio/experiments/casp4/t99/T0099.t2k.a2m @@@@ cp tmp_5-a.train.cst /projects/compbio/experiments/casp4/t99/T0099.t2k.cst @@@@ rm -rf /projects/compbio/tmp/target2k-titanic-29378 < /dev/null gzip -f T0099.t2k.a2m T0099.t2k_*.a2m prettyalign T0099.t2k.a2m.gz -m8 -i -n -L3333 > T0099.t2k.pa SAM: prettyalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:15 gunzip -c T0099.t2k.a2m.gz > tmp.a2m a2m2html -a2m_in tmp.a2m -retrieve > T0099.t2k.pa.html SAM: /projects/compbio/bin/alpha/prettyalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:15 rm tmp.a2m echo ReadNeuralNet /projects/compbio/usr/karplus/predict-2nd/testing/networks/overrep-2500-IDaa13-5-9-7-10-5-10-11-ehl2-seeded4-stride-trained.net > tmp.script echo ReadA2M T0099.t2k.a2m.gz >> tmp.script echo PrintPrediction T0099.t2k.2d 3670-4530-6947 >> tmp.script echo PrintPredictionFasta T0099.t2k.2d.seq >> tmp.script echo PrintRDB T0099.t2k.2d.rdb >> tmp.script predict-2nd < tmp.script # command:# Neural network set to overrep-2500-IDaa13-5-9-7-10-5-10-11-ehl2-seeded4-stride-trained.net # command:# Reading A2M format from T0099.t2k.a2m.gz # Using SequenceWeight EntropyWeight(1.4, 10) # Using regularizer /projects/compbio/lib/recode3.20comp for sequence weight regularizer. # T0099.t2k.a2m with 871 sequences, total weight= 76.5082 avg weight= 0.0878395 46 iterations # AdjustWeights couldn't save exactly 1.4 bits/position, saving 1.3988 bits. Regularizing alignment for T0099.t2k.a2m.gz # After reading T0099.t2k.a2m.gz, have 56 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 T0099.t2k.2d # command:# Network initialization done # Printing prediction in FASTA format to T0099.t2k.2d.seq # command:# Network initialization done # Printing prediction to T0099.t2k.2d.rdb # command:rm tmp.script gunzip -f T0099.t2k.a2m.gz /projects/compbio/usr/karplus/src/phytree/phytree -o -i -r flat T0099.t2k T0099.t2k.a2m /projects/compbio/lib/recode2.20comp Prior library /projects/compbio/lib/recode2.20comp read. Reading alignment file T0099.t2k.a2m (871 sequences) as A2M alignment. Alignment T0099.t2k.a2m read. Adding sequences:0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210 220 230 240 250 260 270 280 290 300 310 320 330 340 350 360 370 380 390 400 410 420 430 440 450 460 470 480 490 500 510 520 530 540 550 560 570 580 590 600 610 620 630 640 650 660 670 680 690 700 710 720 730 740 750 760 770 780 790 800 810 820 830 840 850 860 870 One dot for every 10 merges: ....................................................................................... Merges done. CPU = user 0: 5:48 system 0: 0: 0 gzip -f T0099.t2k.a2m gzip -f T0099.t2k_sorted.a2m rm T0099.t2k.phytrace T0099.t2k.tree_weight prettyalign T0099.t2k_sorted.a2m.gz -m8 -i -n -L3333 > T0099.t2k_sorted.pa SAM: prettyalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:15 /projects/compbio/usr/karplus/src/phytree/dg T0099.t2k.tree /projects/compbio/usr/karplus/src/phytree/dg: 29723 Memory fault /projects/compbio/usr/karplus/src/phytree/dtree T0099.t2k.tree /projects/compbio/usr/karplus/src/phytree/dtree: 29567 Memory fault make[1]: Leaving directory `/auto/projects/compbio/experiments/casp4/t99' /projects/compbio/experiments/models.97/scripts99/score-target -alignment T0099.t2k.a2m.gz -build w0.5 \ -db /projects/compbio/data/pdb/all-protein \ -Emax 40.0 \ -viterbi 0 -sw_score 2 \ -output T0099.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-29879 @@@@ mkdir /projects/compbio/tmp/score-target-29879/for-scorescript @@@@ /projects/compbio/experiments/models.97/scripts99/score-mods-vs-seqs.pl -root_tmp /projects/compbio/tmp/score-target-29879/for-scorescript -a2msfile /projects/compbio/tmp/score-target-29879/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-29879/score-target.op @@@@ chgrp protein /projects/compbio/tmp/score-target-29879/for-scorescript/split-tmp @@@@ chgrp protein /projects/compbio/tmp/score-target-29879/for-scorescript/split-tmp @@@@ /projects/compbio/bin/scripts/w0.5 /projects/compbio/experiments/casp4/t99/T0099.t2k.a2m.gz /projects/compbio/tmp/score-target-29879/for-scorescript/split-tmp/titanic-29717/T0099.t2k.mod Reading /projects/compbio/bin/scripts/sam-t99.conf @@@@ chgrp protein /projects/compbio/tmp/tmp-build-weighted-model-titanic-29557 @@@@ uniqueseq /projects/compbio/tmp/tmp-build-weighted-model-titanic-29557/tmp -alignfile /projects/compbio/experiments/casp4/t99/T0099.t2k.a2m.gz \ -a protein -percent_id 0.8 SAM: uniqueseq v3.1b (February 24, 1999) compiled 04/18/00_11:44:50 Reading alignment file /projects/compbio/experiments/casp4/t99/T0099.t2k.a2m.gz (871 sequences, 56 columns) as A2M alignment. Writing sequence output to /projects/compbio/tmp/tmp-build-weighted-model-titanic-29557/tmp.a2m. Dropping 377 (of 871) duplicate sequences with differing IDs Dropping 215 (of 871) sequences with > 80.0% id 279 sequences left after dropping 592 of 871 sequences. @@@@ modelfromalign /projects/compbio/tmp/tmp-build-weighted-model-titanic-29557/tmp -alignfile /projects/compbio/tmp/tmp-build-weighted-model-titanic-29557/tmp.a2m \ -insert /projects/compbio/lib/fssp-trained.regularizer \ -ins_jump_conf 1 \ -match_jump_conf 1 \ -del_jump_conf 1 \ -prior_library /projects/compbio/lib/recode3.20comp \ -binary_output 1 -a protein -aweight_method 1 -aweight_bits 0.5 -aweight_exponent 10 Reading parameter file /projects/compbio/lib/fssp-trained.regularizer /projects/compbio/lib/fssp-trained.regularizer(1): Reading REGULARIZER: transitions from FSSP and pseudocounts for diverse weighted HSSP SAM: modelfromalign v3.1b (February 24, 1999) compiled 04/18/00_11:44:13 Reading alignment file /projects/compbio/tmp/tmp-build-weighted-model-titanic-29557/tmp.a2m (279 sequences, 56 columns) as A2M alignment. @@@@ cp -f /projects/compbio/tmp/tmp-build-weighted-model-titanic-29557/tmp.mod /projects/compbio/tmp/score-target-29879/for-scorescript/split-tmp/titanic-29717/T0099.t2k.mod @@@@ rm -rf /projects/compbio/tmp/tmp-build-weighted-model-titanic-29557 tmp.script echo ReadA2M T0099.t2k.a2m.gz >> tmp.script echo PrintPrediction T0099.t2k.2d 3670-4530-6947 >> tmp.script echo PrintPredictionFasta T0099.t2k.2d.seq >> tmp.script echo PrintRDB T0099.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 T0099.t2k.a2m.gz # Using SequenceWeight EntropyWeight(1.4, 10) # Using regularizer /projects/compbio/lib/recode3.20comp for sequence weight regularizer. # T0099.t2k.a2m with 871 sequences, total weight= 76.5082 avg weight= 0.0878395 46 iterations # AdjustWeights couldn't save exactly 1.4 bits/position, saving 1.3988 bits. Regularizing alignment for T0099.t2k.a2m.gz # After reading T0099.t2k.a2m.gz, have 56 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 T0099.t2k.2d # command:# Network initialization done # Printing prediction in FASTA format to T0099.t2k.2d.seq # command:# Network initialization done # Printing prediction to T0099.t2k.2d.rdb # command:rm tmp.script