make -k T0128-blast.rdb T0128-double-blast.rdb make[1]: Entering directory `/export/projects/compbio/experiments/casp4/t128' /projects/compbio/experiments/casp4/scripts/single-blast -q T0128.seq -interdb /projects/compbio/data/nrp/nr -db /projects/compbio/data/pdb/all-protein > T0128-blast.rdb starting T0128 with E=10 /projects/compbio/experiments/casp4/scripts/double-blast -q T0128.seq -interdb /projects/compbio/data/nrp/nr -db /projects/compbio/data/pdb/all-protein > T0128-double-blast.rdb starting T0128 with E=0.000050 T0128 T0128 blasting: T0128 gi|134672|sp|P11796|SODM_NICPL_28:216 blasting: T0128 gi|279456|pir||DSPMN_40:237 blasting: T0128 gi|134668|sp|P09233|SODM_MAIZE_32:223 blasting: T0128 gi|34707|emb|CAA32502.1|_24:216 blasting: T0128 gi|34709|emb|CAA42066.1|_24:216 blasting: T0128 gi|34711|emb|CAA68533.1|_24:216 blasting: T0128 gi|36518|emb|CAA30687.1|_24:216 blasting: T0128 gi|134649|sp|P28760|SODM_BACCA_3:201 blasting: T0128 gi|134638|sp|P19685|SODF_COXBU_3:192 blasting: T0128 gi|42727|emb|CAA43108.1|_3:89 blasting: T0128 gi|134659|sp|P00448|SODM_ECOLI_3:201 blasting: T0128 gi|134666|sp|P28763|SODM_LISIV_3:201 blasting: T0128 gi|56691|emb|CAA68549.1|_25:216 blasting: T0128 gi|134663|sp|P09224|SODM_HALHA_3:188 blasting: T0128 gi|77595|pir||S07147_4:203 blasting: T0128 gi|134645|sp|P09223|SODF_PSEPU_2:183 blasting: T0128 gi|401109|sp|P31108|SODF_LEGPN_3:192 blasting: T0128 gi|134643|sp|P09213|SODF_PHOLE_2:191 blasting: T0128 gi|134644|sp|P19665|SODF_PORGI_3:190 blasting: T0128 gi|134647|sp|P18655|SODF_SYNP7_3:190 blasting: T0128 gi|134670|sp|P13367|SODM_MYCLE_4:197 blasting: T0128 gi|134641|sp|P18868|SODF_METTH_7:205 blasting: T0128 gi|82728|pir||S03839_32:223 blasting: T0128 gi|134640|sp|P23744|SODF_METJ_2:201 blasting: T0128 gi|97324|pir||A43585_3:190 blasting: T0128 gi|134667|sp|P28764|SODM_LISMO_3:201 blasting: T0128 gi|134671|sp|P17670|SODF_MYCTU_4:197 blasting: T0128 gi|134662|sp|P09737|SOD1_HALCU_4:195 blasting: T0128 gi|99759|pir||B39267_19:221 blasting: T0128 gi|100284|pir||A39267_26:225 blasting: T0128 gi|134648|sp|P19666|SODF_TETPY_1:195 blasting: T0128 gi|134673|sp|P28765|SODM_PALVU_1:144 blasting: T0128 gi|103414|pir||S23657_2:143 blasting: T0128 gi|134674|sp|P28766|SODM_PARCL_2:144 blasting: T0128 gi|134651|sp|P28761|SODM_BRAFL_2:144 blasting: T0128 gi|134676|sp|P28767|SODM_PETMA_1:144 blasting: T0128 gi|134661|sp|P28762|SODM_EPTST_1:144 blasting: T0128 gi|134677|sp|P28768|SODM_PIG_1:144 blasting: T0128 gi|112660|pir||S12882_3:190 blasting: T0128 gi|134639|sp|P09157|SODF_ECOLI_3:191 blasting: T0128 gi|134642|sp|P22302|SODF_NICPL_1:200 blasting: T0128 gi|134646|sp|P28759|SODF_SOYBN_20:221 blasting: T0128 gi|134646|sp|P28759|SODF_SOYBN_103:241 blasting: T0128 gi|134650|sp|P00449|SODM_BACST_3:201 blasting: T0128 gi|134665|sp|P04179|SODM_HUMAN_24:216 blasting: T0128 gi|1711432|sp|P53638|SODF_BACFR_3:193 blasting: T0128 gi|586007|sp|P37369|SODF_BORPE_3:185 blasting: T0128 gi|417791|sp|Q03302|SODM_HALMA_9:200 blasting: T0128 gi|148814|gb|AAA72704.1|_3:188 blasting: T0128 gi|548947|sp|Q03299|SOD2_HALSG_3:188 blasting: T0128 gi|548946|sp|Q03303|SOD1_HALSG_4:195 blasting: T0128 gi|417789|sp|Q03300|SOD1_HALVO_4:196 blasting: T0128 gi|417790|sp|Q03301|SOD2_HALVO_3:195 blasting: T0128 gi|1351082|sp|P21276|SODF_ARATH_19:221 blasting: T0128 gi|1711455|sp|P53653|SODM_THEAQ_5:198 blasting: T0128 gi|401110|sp|P31161|SODM_CAEEL_25:219 blasting: T0128 gi|223570|prf||0901224A_1:201 blasting: T0128 gi|227056|prf||1613421A_3:191 blasting: T0128 gi|2507403|sp|P09738|SODM_STRMU_5:201 blasting: T0128 gi|3041732|sp|P09671|SODM_MOUSE_25:216 blasting: T0128 gi|464774|sp|P34107|SODF_ENTHI_3:190 blasting: T0128 gi|464775|sp|P35017|SODM_HEVBR_31:219 blasting: T0128 gi|7433348|pir||T04312_28:219 blasting: T0128 gi|349942|pdb|3SDP|A_2:183 blasting: T0128 gi|586008|sp|Q08713|SODF_SULAC_5:210 blasting: T0128 gi|406218|emb|CAA51195.1|_2:188 blasting: T0128 gi|1174382|sp|Q00637|SODM_DROME_18:212 blasting: T0128 gi|1711436|sp|P53641|SODF_PSEAE_3:188 blasting: T0128 gi|628856|pir||S39871_3:206 blasting: T0128 gi|442613|pdb|1ABM|A_1:192 blasting: T0128 gi|448837|prf||1918164A_3:202 blasting: T0128 gi|1174386|sp|P41982|SODM_RABIT_6:197 blasting: T0128 gi|464776|sp|P80293|SODM_PROFR_3:198 blasting: T0128 gi|8928537|sp|P53640|SODF_CAMJE_2:202 blasting: T0128 gi|1711453|sp|P53651|SODM_NOCAS_4:199 blasting: T0128 gi|498260|gb|AAA30655.1|_33:225 blasting: T0128 gi|1174379|sp|P43312|SODF_HELPY_2:188 blasting: T0128 gi|514917|gb|AAA28694.1|_18:212 blasting: T0128 gi|1711433|sp|P53639|SODF_CAMCO_2:202 blasting: T0128 gi|1072933|pir||S47167_2:202 blasting: T0128 gi|531712|emb|CAA56951.1|_38:91 blasting: T0128 gi|533307|dbj|BAA02902.1|_3:193 blasting: T0128 gi|1174391|sp|P41980|SODP_MAIZE_30:221 blasting: T0128 gi|1174390|sp|P41979|SODO_MAIZE_30:221 blasting: T0128 gi|1174389|sp|P41978|SODN_MAIZE_30:220 blasting: T0128 gi|541031|pir||PN0615_2:169 blasting: T0128 gi|541376|pir||PN0614_1:162 blasting: T0128 gi|542489|pir||S41320_25:225 blasting: T0128 gi|1351085|sp|P47201|SODM_MYCAV_4:197 blasting: T0128 gi|1174385|sp|P41981|SODM_ONCVO_25:221 blasting: T0128 gi|630796|pir||S48832_25:221 blasting: T0128 gi|730761|sp|P40726|SODF_SALTY_3:40 blasting: T0128 gi|7433346|pir||T04072_28:219 blasting: T0128 gi|7433347|pir||T04075_28:219 blasting: T0128 gi|1078796|pir||S51750_1:129 blasting: T0128 gi|1711422|sp|P50058|SOD1_PLEBO_46:245 blasting: T0128 gi|1711423|sp|P50059|SOD2_PLEBO_3:205 blasting: T0128 gi|1711424|sp|P50060|SOD3_PLEBO_28:236 blasting: T0128 gi|1711435|sp|P50061|SODF_PLEBO_3:197 blasting: T0128 gi|1174381|sp|P42821|SODM_CORDI_1:162 blasting: T0128 gi|2148070|pir||I40858_1:163 blasting: T0128 gi|620090|emb|CAA57147.1|_1:163 blasting: T0128 gi|620094|emb|CAA57148.1|_1:163 blasting: T0128 gi|620096|emb|CAA57149.1|_1:163 blasting: T0128 gi|1711450|sp|P50913|SODM_MYCSC_1:161 blasting: T0128 gi|620114|emb|CAA57153.1|_1:163 blasting: T0128 gi|6094319|sp|Q60036|SODF_METTM_4:197 blasting: T0128 gi|632260|pir||A54001_4:34 blasting: T0128 gi|2117678|pir||I39745_1:163 blasting: T0128 gi|640118|pdb|1ISC|A_2:190 blasting: T0128 gi|1711444|sp|P46728|SODM_MYCIT_1:163 blasting: T0128 gi|1711445|sp|P50912|SODM_MYCKA_1:163 blasting: T0128 gi|641945|emb|CAA57154.1|_1:146 blasting: T0128 gi|1711441|sp|P53643|SODM_MYCCE_1:138 blasting: T0128 gi|1711442|sp|P53644|SODM_MYCCH_2:138 blasting: T0128 gi|667022|emb|CAA88242.1|_1:138 blasting: T0128 gi|667024|emb|CAA88250.1|_2:138 blasting: T0128 gi|667026|emb|CAA88251.1|_1:138 blasting: T0128 gi|1711446|sp|P53645|SODM_MYCMA_1:138 blasting: T0128 gi|1711447|sp|P53646|SODM_MYCMR_2:138 blasting: T0128 gi|1711448|sp|P53647|SODM_MYCPA_1:138 blasting: T0128 gi|1711449|sp|P53648|SODM_MYCPH_2:138 blasting: T0128 gi|1711451|sp|P53649|SODM_MYCSM_1:138 blasting: T0128 gi|1711452|sp|P53650|SODM_MYCSZ_2:138 blasting: T0128 gi|2117676|pir||JC4611_30:232 blasting: T0128 gi|743362|prf||2012263A_3:195 blasting: T0128 gi|1351086|sp|P41977|SODN_CAEEL_25:216 blasting: T0128 gi|1711440|sp|P50911|SODM_LACLC_3:205 blasting: T0128 gi|2117672|pir||I40871_2:117 blasting: T0128 gi|1075707|pir||S54789_1:120 blasting: T0128 gi|1075699|pir||S54801_1:118 blasting: T0128 gi|1075756|pir||S54792_1:126 blasting: T0128 gi|1075676|pir||S54793_1:115 blasting: T0128 gi|1075775|pir||S54794_4:145 blasting: T0128 gi|1075766|pir||S54795_1:116 blasting: T0128 gi|1174387|sp|P43019|SODM_SALTY_3:201 blasting: T0128 gi|809037|emb|CAA89215.1|_7:123 blasting: T0128 gi|832851|gb|AAB60902.1|_25:216 blasting: T0128 gi|1078795|pir||S55499_4:196 blasting: T0128 gi|913892|gb|AAB34099.1|_4:32 blasting: T0128 gi|1084048|pir||S50999_2:39 blasting: T0128 gi|3287977|sp|P27084|SODM_PEA_40:230 blasting: T0128 gi|2500829|sp|Q59519|SODM_MYCFO_4:197 blasting: T0128 gi|3334336|sp|Q42684|SODM_CHLRE_5:218 blasting: T0128 gi|1711438|sp|P53642|SODM_BORPE_3:203 blasting: T0128 gi|1061202|emb|CAA62838.1|_4:189 blasting: T0128 gi|1351083|sp|P49114|SODM_CAVPO_24:208 blasting: T0128 gi|1066332|gb|AAB39949.1|_3:195 blasting: T0128 gi|1174383|sp|P43725|SODM_HAEIN_3:214 blasting: T0128 gi|1075742|pir||S54791_7:123 blasting: T0128 gi|1076021|pir||S52372_2:137 blasting: T0128 gi|1076033|pir||S52370_1:137 blasting: T0128 gi|1076034|pir||S52379_2:137 blasting: T0128 gi|1076035|pir||S52366_1:137 blasting: T0128 gi|1084032|pir||A53294_3:202 blasting: T0128 gi|1086008|pir||S52360_1:137 blasting: T0128 gi|1086009|pir||S52362_1:137 blasting: T0128 gi|1086010|pir||S52361_2:137 blasting: T0128 gi|1086011|pir||S52364_1:137 blasting: T0128 gi|1086012|pir||S52365_2:137 blasting: T0128 gi|1086013|pir||S52376_2:137 blasting: T0128 gi|1129044|gb|AAA83933.1|_24:71 blasting: T0128 gi|1711457|sp|P53654|SODM_XANCP_3:202 blasting: T0128 gi|7433339|pir||T10943_9:99 blasting: T0128 gi|2500830|sp|Q59679|SODM_PASHA_3:204 blasting: T0128 gi|1174380|sp|P41976|SODM_BOVIN_24:216 blasting: T0128 gi|1203805|gb|AAC43639.1|_3:208 blasting: T0128 gi|1711458|sp|P53655|SODM_YEREN_3:206 blasting: T0128 gi|1263909|emb|CAA65734.1|_3:44 blasting: T0128 gi|1711437|sp|P54375|SODM_BACSU_3:204 blasting: T0128 gi|6685985|sp|Q59094|SODM_ACICA_29:226 blasting: T0128 gi|7433340|pir||T09788_29:215 blasting: T0128 gi|1399534|gb|AAB03279.1|_3:204 blasting: T0128 gi|1401274|gb|AAB03399.1|_1:87 blasting: T0128 gi|2500825|sp|P77928|SODG_PSEPU_3:192 blasting: T0128 gi|2500831|sp|P77929|SODM_PSEPU_3:202 blasting: T0128 gi|1515346|gb|AAB17024.1|_4:183 blasting: T0128 gi|1518972|gb|AAB16766.1|_1:142 blasting: T0128 gi|1518974|gb|AAB16767.1|_1:142 blasting: T0128 gi|1518976|gb|AAB16768.1|_1:142 blasting: T0128 gi|1518978|gb|AAB16769.1|_1:142 blasting: T0128 gi|1518980|gb|AAB16770.1|_1:142 blasting: T0128 gi|1518982|gb|AAB16771.1|_1:142 blasting: T0128 gi|1518984|gb|AAB16772.1|_1:142 blasting: T0128 gi|1518986|gb|AAB16773.1|_1:142 blasting: T0128 gi|1518988|gb|AAB16774.1|_1:142 blasting: T0128 gi|1518990|gb|AAB16775.1|_1:142 blasting: T0128 gi|1518992|gb|AAB16776.1|_1:142 blasting: T0128 gi|1518994|gb|AAB16777.1|_1:142 blasting: T0128 gi|1518996|gb|AAB16778.1|_1:142 blasting: T0128 gi|1518998|gb|AAB16779.1|_1:142 blasting: T0128 gi|1519000|gb|AAB16780.1|_1:142 blasting: T0128 gi|1519002|gb|AAB16781.1|_1:142 blasting: T0128 gi|1519006|gb|AAB16783.1|_1:142 blasting: T0128 gi|1519008|gb|AAB16784.1|_1:142 blasting: T0128 gi|1519010|gb|AAB16785.1|_1:142 blasting: T0128 gi|1519012|gb|AAB16786.1|_1:142 blasting: T0128 gi|1519016|gb|AAB16788.1|_1:142 blasting: T0128 gi|1519018|gb|AAB16789.1|_1:142 blasting: T0128 gi|1519020|gb|AAB16790.1|_1:142 blasting: T0128 gi|1519022|gb|AAB16791.1|_1:142 blasting: T0128 gi|1519024|gb|AAB16792.1|_1:142 blasting: T0128 gi|1519028|gb|AAB16794.1|_1:142 blasting: T0128 gi|1519030|gb|AAB16795.1|_1:142 blasting: T0128 gi|1519034|gb|AAB16797.1|_1:142 blasting: T0128 gi|2500828|sp|Q92429|SODM_GANMI_3:198 blasting: T0128 gi|7433342|pir||T06258_31:230 blasting: T0128 gi|2507402|sp|P53652|SODM_PSEAE_3:202 blasting: T0128 gi|1633498|pdb|1VAR|A_2:192 blasting: T0128 gi|3334338|sp|Q92450|SODM_ASPFU_13:207 blasting: T0128 gi|1649008|gb|AAB17667.1|_2:98 blasting: T0128 gi|1649010|gb|AAB17668.1|_59:126 blasting: T0128 gi|1649010|gb|AAB17668.1|_2:36 blasting: T0128 gi|2500827|sp|P77968|SODF_SYNY3_3:199 blasting: T0128 gi|7433345|pir||T06801_31:230 blasting: T0128 gi|1669755|emb|CAA94532.1|_1:93 blasting: T0128 gi|1669761|emb|CAA94533.1|_1:124 blasting: T0128 gi|1669765|emb|CAA94527.1|_1:105 blasting: T0128 gi|1669767|emb|CAA94528.1|_1:105 blasting: T0128 gi|1669769|emb|CAA94529.1|_1:93 blasting: T0128 gi|1669771|emb|CAA94530.1|_1:93 blasting: T0128 gi|1669773|emb|CAA94531.1|_1:105 blasting: T0128 gi|1711456|sp|P09214|SODM_THETH_5:204 blasting: T0128 gi|1763271|gb|AAB39705.1|_1:115 blasting: T0128 gi|1835757|gb|AAB46719.1|_1:115 blasting: T0128 gi|1835759|gb|AAB46720.1|_1:115 blasting: T0128 gi|1835761|gb|AAB46721.1|_1:115 blasting: T0128 gi|1835763|gb|AAB46722.1|_1:115 blasting: T0128 gi|1835765|gb|AAB46723.1|_1:115 blasting: T0128 gi|1835767|gb|AAB46724.1|_1:115 blasting: T0128 gi|1835769|gb|AAB46725.1|_1:115 blasting: T0128 gi|1915999|gb|AAC47734.1|_2:192 blasting: T0128 gi|1916001|gb|AAC47735.1|_2:192 blasting: T0128 gi|2072628|emb|CAA73188.1|_1:116 blasting: T0128 gi|2117680|pir||JC4351_4:199 blasting: T0128 gi|2117681|pir||JC4396_4:199 blasting: T0128 gi|2144131|pir||S65795_17:201 blasting: T0128 gi|2149610|gb|AAC47548.1|_36:205 blasting: T0128 gi|2149612|gb|AAC47549.1|_3:195 blasting: T0128 gi|7674339|sp|O08459|SODM_HALHI_1:143 blasting: T0128 gi|7674340|sp|O08460|SOD2_HALME_1:145 blasting: T0128 gi|7674341|sp|O08461|SOD1_HALME_1:142 blasting: T0128 gi|2290275|gb|AAB65095.1|_1:136 blasting: T0128 gi|3915001|sp|O30563|SODM_BORBU_2:203 blasting: T0128 gi|6226150|sp|O33603|SODM_STRAI_1:119 blasting: T0128 gi|2293018|emb|CAB09346.1|_1:126 blasting: T0128 gi|6226151|sp|O33605|SODM_STRAY_1:116 blasting: T0128 gi|2293022|emb|CAB09348.1|_1:116 blasting: T0128 gi|2293024|emb|CAB09349.1|_1:116 blasting: T0128 gi|2293026|emb|CAB09350.1|_1:116 blasting: T0128 gi|2293028|emb|CAB09351.1|_1:116 blasting: T0128 gi|2293030|emb|CAB09352.1|_1:145 blasting: T0128 gi|2293032|emb|CAB09353.1|_1:116 blasting: T0128 gi|2293034|emb|CAB09354.1|_1:145 blasting: T0128 gi|2293036|emb|CAB09355.1|_1:145 blasting: T0128 gi|2293038|emb|CAB09356.1|_1:116 blasting: T0128 gi|2293042|emb|CAB09358.1|_1:116 blasting: T0128 gi|2293044|emb|CAB09359.1|_1:116 blasting: T0128 gi|2293046|emb|CAB09360.1|_1:116 blasting: T0128 gi|2293048|emb|CAB09361.1|_1:116 blasting: T0128 gi|3915015|sp|Q59949|SODM_STRPN_1:116 blasting: T0128 gi|2293052|emb|CAB09363.1|_1:145 blasting: T0128 gi|2293054|emb|CAB09364.1|_1:116 blasting: T0128 gi|2293056|emb|CAB09365.1|_1:116 blasting: T0128 gi|6226773|sp|O33756|SODM_STRPA_1:113 blasting: T0128 gi|2293085|emb|CAB09369.1|_1:145 blasting: T0128 gi|6226774|sp|O33783|SODM_STRSL_1:145 blasting: T0128 gi|2293089|emb|CAB09371.1|_1:116 blasting: T0128 gi|2293091|emb|CAB09372.1|_1:116 blasting: T0128 gi|2293093|emb|CAB09373.1|_1:145 blasting: T0128 gi|2293095|emb|CAB09374.1|_1:145 blasting: T0128 gi|2293097|emb|CAB09375.1|_1:145 blasting: T0128 gi|7433331|pir||E64568_2:188 blasting: T0128 gi|2352423|gb|AAC38829.1|_3:195 blasting: T0128 gi|2352425|gb|AAC38830.1|_35:228 blasting: T0128 gi|2384558|gb|AAB69755.1|_3:193 blasting: T0128 gi|6094318|sp|O15905|SODF_BABBO_3:193 blasting: T0128 gi|2398587|emb|CAA72335.1|_24:223 blasting: T0128 gi|6094321|sp|O30826|SODM_HAEDU_3:208 blasting: T0128 gi|2500826|sp|P80857|SODF_SULSO_7:208 blasting: T0128 gi|2511749|gb|AAC64207.1|_3:200 blasting: T0128 gi|6647843|sp|O35023|SODF_BACSU_79:273 blasting: T0128 gi|3334330|sp|O13401|SODM_CANAL_32:233 blasting: T0128 gi|2654175|gb|AAB87697.1|_1:196 blasting: T0128 gi|3915002|sp|O49066|SODM_CAPAN_28:216 blasting: T0128 gi|7433333|pir||A70119_13:214 blasting: T0128 gi|3915003|sp|O51917|SODF_STRCO_4:202 blasting: T0128 gi|2738984|gb|AAC48291.1|_2:192 blasting: T0128 gi|2765187|emb|CAA72899.1|_5:200 blasting: T0128 gi|6226152|sp|O54210|SODM_STRCB_4:145 blasting: T0128 gi|6226153|sp|O54233|SODM_STRIN_5:116 blasting: T0128 gi|6226154|sp|O54266|SODM_STRPO_5:145 blasting: T0128 gi|2765853|emb|CAB16323.1|_1:126 blasting: T0128 gi|2765855|emb|CAB16324.1|_1:126 blasting: T0128 gi|2765857|emb|CAB16325.1|_1:119 blasting: T0128 gi|2765859|emb|CAB16326.1|_1:116 blasting: T0128 gi|2765861|emb|CAB16327.1|_1:116 blasting: T0128 gi|2765863|emb|CAB16328.1|_1:116 blasting: T0128 gi|2765869|emb|CAB16331.1|_1:116 blasting: T0128 gi|2765877|emb|CAB16335.1|_1:116 blasting: T0128 gi|2765879|emb|CAB16336.1|_1:116 blasting: T0128 gi|2765885|emb|CAB16339.1|_1:116 blasting: T0128 gi|2765887|emb|CAB16340.1|_1:110 blasting: T0128 gi|2765895|emb|CAB16344.1|_1:116 blasting: T0128 gi|2765903|emb|CAB16348.1|_1:116 blasting: T0128 gi|2780818|pdb|1QNM|A_1:192 blasting: T0128 gi|2808462|emb|CAA11227.1|_5:200 blasting: T0128 gi|2914416|pdb|1AP6|A_1:192 blasting: T0128 gi|8478425|sp|O67470|SODF_AQUAE_20:209 blasting: T0128 gi|7433350|pir||T08181_24:219 blasting: T0128 gi|3114705|gb|AAC15842.1|_5:207 blasting: T0128 gi|7245492|pdb|1D5N|A_2:200 blasting: T0128 gi|7484734|pir||T09799_6:186 blasting: T0128 gi|8134711|sp|O86165|SODM_MYCLP_4:197 blasting: T0128 gi|3273560|gb|AAC24764.1|_24:220 blasting: T0128 gi|3273751|gb|AAC24832.1|_31:231 blasting: T0128 gi|3273757|gb|AAC24834.1|_50:249 blasting: T0128 gi|3320377|gb|AAC26483.1|_3:201 blasting: T0128 gi|7433324|pir||H71531_6:199 blasting: T0128 gi|3355882|dbj|BAA31974.1|_3:201 blasting: T0128 gi|7493016|pir||T40411_69:213 blasting: T0128 gi|3599469|gb|AAC35356.1|_1:175 blasting: T0128 gi|7388256|sp|O75007|SODM_PENCH_1:210 blasting: T0128 gi|3676839|gb|AAC62115.1|_25:223 blasting: T0128 gi|3688568|emb|CAA75583.1|_8:209 blasting: T0128 gi|3719455|gb|AAC63378.1|_48:177 blasting: T0128 gi|3719455|gb|AAC63378.1|_122:260 blasting: T0128 gi|3719457|gb|AAC63379.1|_40:228 blasting: T0128 gi|3746358|gb|AAC63943.1|_3:200 blasting: T0128 gi|3790114|emb|CAA05291.1|_3:194 blasting: T0128 gi|7674342|sp|Q9ZD15|SODF_RICPR_14:207 blasting: T0128 gi|4090354|gb|AAC98799.1|_2:139 blasting: T0128 gi|4090356|gb|AAC98800.1|_2:139 blasting: T0128 gi|4102997|gb|AAD01640.1|_20:191 blasting: T0128 gi|4106403|gb|AAD02836.1|_3:187 blasting: T0128 gi|4139274|gb|AAD03746.1|_4:196 blasting: T0128 gi|4139304|gb|AAD03761.1|_4:196 blasting: T0128 gi|4139312|gb|AAD03765.1|_4:196 blasting: T0128 gi|4139314|gb|AAD03766.1|_4:196 blasting: T0128 gi|4139316|gb|AAD03767.1|_4:196 blasting: T0128 gi|4139318|gb|AAD03768.1|_4:196 blasting: T0128 gi|4139320|gb|AAD03769.1|_4:196 blasting: T0128 gi|4139322|gb|AAD03770.1|_4:196 blasting: T0128 gi|4139324|gb|AAD03771.1|_4:196 blasting: T0128 gi|9789808|sp|Q9ZKE6|SODF_HELPJ_2:188 blasting: T0128 gi|7433337|pir||JG0179_43:241 blasting: T0128 gi|4218195|emb|CAA10729.1|_2:188 blasting: T0128 gi|4218197|emb|CAA10730.1|_2:188 blasting: T0128 gi|4218199|emb|CAA10731.1|_2:188 blasting: T0128 gi|4321778|gb|AAD15825.1|_4:197 blasting: T0128 gi|4325247|gb|AAD17309.1|_3:192 blasting: T0128 gi|7433344|pir||B72124_6:207 blasting: T0128 gi|7433338|pir||T05531_27:229 blasting: T0128 gi|4456296|emb|CAB36929.1|_2:188 blasting: T0128 gi|6226706|sp|P17550|SODM_ALCEU_18:184 blasting: T0128 gi|4581922|gb|AAD24796.1|AF121078_2_3:202 blasting: T0128 gi|4581924|gb|AAD24797.1|AF121079_1_3:190 blasting: T0128 gi|4585263|gb|AAD25353.1|AF114848_1_4:203 blasting: T0128 gi|9789802|sp|Q9Y8H8|SODF_AERPE_2:211 blasting: T0128 gi|4589878|dbj|BAA76922.1|_18:216 blasting: T0128 gi|4731567|gb|AAD28503.1|AF118809_1_34:231 blasting: T0128 gi|4809202|gb|AAD30139.1|AF141866_1_4:202 blasting: T0128 gi|4927474|gb|AAD33130.1|AF099015_2_4:200 blasting: T0128 gi|8478444|sp|Q9X6W9|SODF_AQUPY_20:209 blasting: T0128 gi|4930136|pdb|1SSS|A_4:205 blasting: T0128 gi|7388259|sp|Q9XD74|SODF_RHIME_3:191 blasting: T0128 gi|5453298|gb|AAD43522.1|_4:196 blasting: T0128 gi|5453299|gb|AAD43523.1|_4:196 blasting: T0128 gi|5542211|pdb|1COJ|A_19:208 blasting: T0128 gi|5758312|gb|AAD50778.1|AF162664_1_5:194 blasting: T0128 gi|5771529|gb|AAD51417.1|AF173990_1_7:200 blasting: T0128 gi|5777414|emb|CAB53458.1|_3:191 blasting: T0128 gi|5803267|dbj|BAA83577.1|_43:245 blasting: T0128 gi|5880381|gb|AAD54651.1|AF085191_1_8:199 blasting: T0128 gi|5902898|dbj|BAA84480.1|_1:151 blasting: T0128 gi|5902904|dbj|BAA84483.1|_1:151 blasting: T0128 gi|5902906|dbj|BAA84484.1|_1:151 blasting: T0128 gi|5902908|dbj|BAA84485.1|_1:151 blasting: T0128 gi|5902912|dbj|BAA84487.1|_1:151 blasting: T0128 gi|5902914|dbj|BAA84488.1|_1:151 blasting: T0128 gi|5902916|dbj|BAA84489.1|_1:151 blasting: T0128 gi|6006619|emb|CAB56851.1|_28:216 blasting: T0128 gi|6016702|gb|AAF01529.1|AC009991_25_31:231 blasting: T0128 gi|6069869|dbj|BAA85490.1|_1:101 blasting: T0128 gi|6069876|dbj|BAA85491.1|_1:101 blasting: T0128 gi|6069890|dbj|BAA85493.1|_1:101 blasting: T0128 gi|6069904|dbj|BAA85495.1|_3:101 blasting: T0128 gi|6069911|dbj|BAA85496.1|_1:101 blasting: T0128 gi|6069917|dbj|BAA85497.1|_3:101 blasting: T0128 gi|6069924|dbj|BAA85498.1|_1:101 blasting: T0128 gi|6069930|dbj|BAA85499.1|_1:101 blasting: T0128 gi|6070008|dbj|BAA85504.1|_1:101 blasting: T0128 gi|6070041|dbj|BAA85506.1|_1:101 blasting: T0128 gi|6070132|dbj|BAA85509.1|_1:101 blasting: T0128 gi|6070138|dbj|BAA85510.1|_1:101 blasting: T0128 gi|6070151|dbj|BAA85512.1|_1:101 blasting: T0128 gi|6070157|dbj|BAA85513.1|_1:101 blasting: T0128 gi|6070163|dbj|BAA85514.1|_1:101 blasting: T0128 gi|6070168|dbj|BAA85515.1|_1:101 blasting: T0128 gi|6070185|dbj|BAA85518.1|_1:101 blasting: T0128 gi|6070192|dbj|BAA85519.1|_1:101 blasting: T0128 gi|6070198|dbj|BAA85520.1|_1:101 blasting: T0128 gi|6070204|dbj|BAA85521.1|_1:101 blasting: T0128 gi|6070257|dbj|BAA85523.1|_1:101 blasting: T0128 gi|6070259|dbj|BAA85524.1|_1:101 blasting: T0128 gi|6070261|dbj|BAA85525.1|_1:101 blasting: T0128 gi|6070265|dbj|BAA85527.1|_1:101 blasting: T0128 gi|6070279|dbj|BAA85534.1|_1:101 blasting: T0128 gi|6070281|dbj|BAA85535.1|_1:101 blasting: T0128 gi|6070287|dbj|BAA85538.1|_1:101 blasting: T0128 gi|6070293|dbj|BAA85541.1|_1:101 blasting: T0128 gi|6070295|dbj|BAA85542.1|_1:101 blasting: T0128 gi|6070301|dbj|BAA85545.1|_1:101 blasting: T0128 gi|6070303|dbj|BAA85546.1|_1:101 blasting: T0128 gi|6070309|dbj|BAA85549.1|_1:101 blasting: T0128 gi|6070311|dbj|BAA85550.1|_1:101 blasting: T0128 gi|6070313|dbj|BAA85551.1|_1:101 blasting: T0128 gi|6070315|dbj|BAA85552.1|_1:101 blasting: T0128 gi|6070327|dbj|BAA85558.1|_1:101 blasting: T0128 gi|6070331|dbj|BAA85560.1|_1:101 blasting: T0128 gi|6070333|dbj|BAA85561.1|_1:101 blasting: T0128 gi|6070348|dbj|BAA85564.1|_1:101 blasting: T0128 gi|6070356|dbj|BAA85568.1|_1:101 blasting: T0128 gi|6070368|dbj|BAA85574.1|_1:101 blasting: T0128 gi|6070370|dbj|BAA85575.1|_1:101 blasting: T0128 gi|6070372|dbj|BAA85576.1|_1:101 blasting: T0128 gi|6070374|dbj|BAA85577.1|_3:101 blasting: T0128 gi|6070380|dbj|BAA85580.1|_3:101 blasting: T0128 gi|6070382|dbj|BAA85581.1|_1:101 blasting: T0128 gi|6070384|dbj|BAA85582.1|_1:101 blasting: T0128 gi|6070386|dbj|BAA85583.1|_1:101 blasting: T0128 gi|6077108|dbj|BAA85501.1|_1:101 blasting: T0128 gi|6077110|dbj|BAA85505.1|_1:101 blasting: T0128 gi|6077112|dbj|BAA85507.1|_1:101 blasting: T0128 gi|6226770|sp|O54086|SODM_STRAG_5:200 blasting: T0128 gi|6226772|sp|O54263|SODM_STROR_1:116 blasting: T0128 gi|6429239|dbj|BAA86881.1|_31:223 blasting: T0128 gi|6453337|emb|CAB61430.1|_2:188 blasting: T0128 gi|7473947|pir||B75415_3:210 blasting: T0128 gi|6523771|emb|CAB62411.1|_22:214 blasting: T0128 gi|6562170|emb|CAB62521.1|_49:108 blasting: T0128 gi|6572119|emb|CAA71559.2|_20:221 blasting: T0128 gi|6573494|pdb|1B06|A_4:209 blasting: T0128 gi|6647851|sp|O93724|SODF_PYRAE_2:211 blasting: T0128 gi|6653269|gb|AAF22650.1|AF201088_1_1:145 blasting: T0128 gi|6653271|gb|AAF22651.1|AF201089_1_1:140 blasting: T0128 gi|6653273|gb|AAF22652.1|AF201090_1_3:168 blasting: T0128 gi|6687499|emb|CAB64955.1|_1:118 blasting: T0128 gi|6687501|emb|CAB64959.1|_1:118 blasting: T0128 gi|6687503|emb|CAB64961.1|_2:120 blasting: T0128 gi|6687505|emb|CAB64964.1|_2:120 blasting: T0128 gi|6687507|emb|CAB64966.1|_1:118 blasting: T0128 gi|6687509|emb|CAB64967.1|_1:120 blasting: T0128 gi|6687511|emb|CAB64979.1|_1:118 blasting: T0128 gi|6687513|emb|CAB64980.1|_1:120 blasting: T0128 gi|6687515|emb|CAB64981.1|_1:118 blasting: T0128 gi|6687525|emb|CAB64970.1|_1:120 blasting: T0128 gi|6687527|emb|CAB64960.1|_1:118 blasting: T0128 gi|6687529|emb|CAB64984.1|_1:120 blasting: T0128 gi|6687533|emb|CAB64986.1|_1:120 blasting: T0128 gi|6687537|emb|CAB64988.1|_1:120 blasting: T0128 gi|6687570|emb|CAB64989.1|_1:118 blasting: T0128 gi|6687572|emb|CAB64991.1|_1:118 blasting: T0128 gi|6687574|emb|CAB64993.1|_1:120 blasting: T0128 gi|6687576|emb|CAB64996.1|_1:118 blasting: T0128 gi|6687578|emb|CAB64997.1|_1:120 blasting: T0128 gi|6687580|emb|CAB65005.1|_1:119 blasting: T0128 gi|6687582|emb|CAB65009.1|_1:117 blasting: T0128 gi|6687584|emb|CAB65011.1|_1:120 blasting: T0128 gi|6687586|emb|CAB65006.1|_1:118 blasting: T0128 gi|6687588|emb|CAB65010.1|_1:120 blasting: T0128 gi|6687590|emb|CAB65012.1|_2:118 blasting: T0128 gi|6687596|emb|CAB64992.1|_2:118 blasting: T0128 gi|6693856|gb|AAF25009.1|AF177945_1_7:200 blasting: T0128 gi|6707930|gb|AAF25723.1|AF146751_1_1:174 blasting: T0128 gi|6707932|gb|AAF25724.1|AF146752_1_1:172 blasting: T0128 gi|6707934|gb|AAF25725.1|AF146753_1_1:172 blasting: T0128 gi|6707936|gb|AAF25726.1|AF146754_1_2:165 blasting: T0128 gi|6707938|gb|AAF25727.1|AF146755_1_3:173 blasting: T0128 gi|6840824|gb|AAF28773.1|AF077224_1_19:220 blasting: T0128 gi|6840824|gb|AAF28773.1|AF077224_1_102:240 blasting: T0128 gi|7110573|gb|AAF36989.1|AF236110_1_8:207 blasting: T0128 gi|7190607|gb|AAF39404.1|_8:201 blasting: T0128 gi|7226122|gb|AAF41295.1|_3:190 blasting: T0128 gi|7248414|dbj|BAA92737.1|_28:231 blasting: T0128 gi|7290015|gb|AAD00533.2|_2:211 blasting: T0128 gi|7328179|emb|CAB82579.1|_6:200 blasting: T0128 gi|7379798|emb|CAB84366.1|_3:190 blasting: T0128 gi|7413493|emb|CAB85688.1|_1:124 blasting: T0128 gi|7433329|pir||S20019_3:201 blasting: T0128 gi|7546411|pdb|1EM1|A_1:192 blasting: T0128 gi|7555818|gb|AAC60522.2|_24:216 blasting: T0128 gi|7572933|emb|CAB87434.1|_34:224 blasting: T0128 gi|7576316|dbj|BAA94550.1|_1:120 blasting: T0128 gi|7576318|dbj|BAA94551.1|_1:120 blasting: T0128 gi|7576320|dbj|BAA94552.1|_2:107 blasting: T0128 gi|7576322|dbj|BAA94553.1|_1:120 blasting: T0128 gi|7576324|dbj|BAA94554.1|_2:107 blasting: T0128 gi|7576326|dbj|BAA94555.1|_2:107 blasting: T0128 gi|7576328|dbj|BAA94556.1|_1:120 blasting: T0128 gi|7576330|dbj|BAA94557.1|_2:107 blasting: T0128 gi|7576334|dbj|BAA94559.1|_1:120 blasting: T0128 gi|7576370|dbj|BAA94560.1|_1:120 blasting: T0128 gi|7576374|dbj|BAA94562.1|_2:107 blasting: T0128 gi|7576376|dbj|BAA94563.1|_2:107 blasting: T0128 gi|7576380|dbj|BAA94565.1|_1:120 blasting: T0128 gi|7576382|dbj|BAA94566.1|_1:120 blasting: T0128 gi|7576386|dbj|BAA94568.1|_2:107 blasting: T0128 gi|7576396|dbj|BAA94573.1|_2:107 blasting: T0128 gi|7576398|dbj|BAA94574.1|_1:120 blasting: T0128 gi|7576400|dbj|BAA94575.1|_2:107 blasting: T0128 gi|7576408|dbj|BAA94579.1|_2:107 blasting: T0128 gi|7576410|dbj|BAA94580.1|_2:107 blasting: T0128 gi|7576414|dbj|BAA94582.1|_1:120 blasting: T0128 gi|7576416|dbj|BAA94583.1|_2:107 blasting: T0128 gi|7576430|dbj|BAA94590.1|_2:107 blasting: T0128 gi|7576432|dbj|BAA94591.1|_2:107 blasting: T0128 gi|7649157|gb|AAF65768.1|AF242310_1_34:223 blasting: T0128 gi|7677029|gb|AAF66995.1|AF062654_1_26:221 blasting: T0128 gi|7798571|dbj|BAA95631.1|_3:201 blasting: T0128 gi|8347758|gb|AAF74770.1|AF264029_1_22:217 blasting: T0128 gi|8347760|gb|AAF74771.1|AF264030_1_86:280 blasting: T0128 gi|8573766|emb|CAB94731.1|_3:197 blasting: T0128 gi|8843846|dbj|BAA97372.1|_48:267 blasting: T0128 gi|8918742|dbj|BAA97805.1|_3:201 blasting: T0128 gi|8918744|dbj|BAA97806.1|_3:201 blasting: T0128 gi|8918746|dbj|BAA97807.1|_3:201 blasting: T0128 gi|8977980|emb|CAB95744.1|_3:192 blasting: T0128 gi|9107021|gb|AAF84727.1|AE004012_1_34:237 blasting: T0128 gi|9107839|gb|AAF85411.1|AE004068_7_3:203 blasting: T0128 gi|9255849|gb|AAF86329.1|_16:124 blasting: T0128 gi|9656592|gb|AAF95193.1|_3:194 blasting: T0128 gi|9657290|gb|AAF95835.1|_17:218 blasting: make[1]: Leaving directory `/export/projects/compbio/experiments/casp4/t128' make -k T0128.t2k.a2m.gz \ T0128.t2k.pa \ T0128.t2k.pa.html \ T0128.t2k.2d \ T0128.t2k.tree \ T0128.t2k_sorted.pa \ T0128.t2k.tree.ps \ T0128.t2k.tree-unroot.ps make[1]: Entering directory `/export/projects/compbio/experiments/casp4/t128' /projects/compbio/experiments/models.97/scripts2k/target2k -out T0128.t2k \ -seed T0128.seq \ -all @@@@ chgrp protein /projects/compbio/tmp/target2k-beta-20339 @@@@ cp /export/projects/compbio/experiments/casp4/t128/T0128.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-23599 /projects/compbio/experiments/models.97/scripts2k/blast-prefilter: E-value thresholds used are For /export/projects/compbio/tmp/target2k-beta-20339/prefilter_1.fa, 0.01 For /export/projects/compbio/tmp/target2k-beta-20339/prefilter_2.fa, 1 For /export/projects/compbio/tmp/target2k-beta-20339/prefilter_3.fa, 10 For /export/projects/compbio/tmp/target2k-beta-20339/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-20339/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-20339/m0.a2m (1 sequences, 222 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-23599 < /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, 222 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, 222 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.0900900900900901 -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: -90.749657 Database has 538 sequences with 92985 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 (523 sequences, 222 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=524 nll_thresh=-22.4115278685381 frac_id=0.684259627508335 @@@@ uniqueseq thinned_tmp -alignfile tmp_1-a.train.seq -a protein \ -percent_id 0.684259627508335 -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file tmp_1-a.train.seq (524 sequences, 222 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 39 (of 524) duplicate sequences with differing IDs Dropping 385 (of 524) sequences with > 68.4% id 100 sequences left after dropping 424 of 524 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 (100 sequences, 222 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 -643.81 26.94 -205.92 70.90 13 0 223 @@@@ 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.0900900900900901 \ -align_short 5 -mdEmax 0.0001 -constraints tmp_1-a.train.cst -constraints_out /export/projects/compbio/experiments/casp4/t128/T0128.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: -190.309280 Database has 533 sequences with 92396 residues. @@@@ cp -f tmp_1-b.mult /export/projects/compbio/experiments/casp4/t128/T0128.t2k_1.a2m # About to count sequences in prefilter_2.fa @@@@ modelfromalign tmp_2-a -alignfile /export/projects/compbio/experiments/casp4/t128/T0128.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/t128/T0128.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/t128/T0128.t2k_1.a2m (532 sequences, 222 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.0900900900900901 -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/t128/T0128.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: -171.961624 Database has 563 sequences with 95872 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 (539 sequences, 222 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=540 nll_thresh=-21.7183806877933 frac_id=0.679371403184734 @@@@ uniqueseq thinned_tmp -alignfile tmp_2-a.train.seq -a protein \ -percent_id 0.679371403184734 -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file tmp_2-a.train.seq (540 sequences, 222 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 12 (of 540) duplicate sequences with differing IDs Dropping 434 (of 540) sequences with > 67.9% id 94 sequences left after dropping 446 of 540 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 (94 sequences, 222 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/t128/T0128.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 -641.82 4.96 -187.49 64.01 9 0 223 @@@@ 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.0900900900900901 \ -align_short 5 -mdEmax 0.0002 -constraints tmp_2-a.train.cst -constraints_out /export/projects/compbio/experiments/casp4/t128/T0128.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: -166.307663 Database has 543 sequences with 93780 residues. @@@@ cp -f tmp_2-b.mult /export/projects/compbio/experiments/casp4/t128/T0128.t2k_2.a2m # About to count sequences in prefilter_3.fa @@@@ modelfromalign tmp_3-a -alignfile /export/projects/compbio/experiments/casp4/t128/T0128.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/t128/T0128.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/t128/T0128.t2k_2.a2m (540 sequences, 222 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.0900900900900901 -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/t128/T0128.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: -152.463608 Database has 606 sequences with 103476 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 (540 sequences, 222 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=541 nll_thresh=-20.1089427738805 frac_id=0.673495658274573 @@@@ uniqueseq thinned_tmp -alignfile tmp_3-a.train.seq -a protein \ -percent_id 0.673495658274573 -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file tmp_3-a.train.seq (541 sequences, 222 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 11 (of 541) duplicate sequences with differing IDs Dropping 440 (of 541) sequences with > 67.3% id 90 sequences left after dropping 451 of 541 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 (90 sequences, 222 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/t128/T0128.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 -641.57 -54.74 -183.02 57.94 7 0 223 @@@@ 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.0900900900900901 \ -align_short 5 -mdEmax 0.001 -constraints tmp_3-a.train.cst -constraints_out /export/projects/compbio/experiments/casp4/t128/T0128.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: -157.160461 Database has 545 sequences with 93925 residues. @@@@ cp -f tmp_3-b.mult /export/projects/compbio/experiments/casp4/t128/T0128.t2k_3.a2m # About to count sequences in prefilter_4.fa @@@@ modelfromalign tmp_4-a -alignfile /export/projects/compbio/experiments/casp4/t128/T0128.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/t128/T0128.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/t128/T0128.t2k_3.a2m (541 sequences, 222 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.0900900900900901 -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/t128/T0128.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: -86.798233 Database has 1043 sequences with 200016 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 (540 sequences, 222 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=541 nll_thresh=-18.4995048540528 frac_id=0.667731334062566 @@@@ uniqueseq thinned_tmp -alignfile tmp_4-a.train.seq -a protein \ -percent_id 0.667731334062566 -a2mdots 0 SAM: uniqueseq v3.2 (July 31, 2000) compiled 07/31/00_14:08:59 Reading alignment file tmp_4-a.train.seq (541 sequences, 222 columns) as A2M alignment. Writing sequence output to thinned_tmp.a2m. Dropping 9 (of 541) duplicate sequences with differing IDs Dropping 441 (of 541) sequences with > 66.8% id 91 sequences left after dropping 450 of 541 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 (91 sequences, 222 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/t128/T0128.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 -641.06 -54.00 -174.58 58.09 8 0 223 @@@@ 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.0900900900900901 \ -align_short 5 -mdEmax 0.005 -constraints tmp_4-a.train.cst -constraints_out /export/projects/compbio/experiments/casp4/t128/T0128.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: -147.381470 Database has 547 sequences with 93967 residues. @@@@ cp -f tmp_4-b.mult /export/projects/compbio/experiments/casp4/t128/T0128.t2k_4.a2m @@@@ modelfromalign tmp_5-a -alignfile /export/projects/compbio/experiments/casp4/t128/T0128.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/t128/T0128.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/t128/T0128.t2k_4.a2m (541 sequences, 222 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.0900900900900901 -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/t128/T0128.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: -87.045021 Database has 1043 sequences with 200016 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 (534 sequences, 222 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/t128/T0128.t2k.a2m @@@@ cp tmp_5-a.train.cst /export/projects/compbio/experiments/casp4/t128/T0128.t2k.cst @@@@ rm -rf /projects/compbio/tmp/target2k-beta-20339 < /dev/null gzip -f T0128.t2k.a2m T0128.t2k_*.a2m prettyalign T0128.t2k.a2m.gz -m8 -i -n -L3333 > T0128.t2k.pa SAM: prettyalign v3.2 (July 31, 2000) compiled 08/11/00_16:27:51 gunzip -c T0128.t2k.a2m.gz > tmp.a2m a2m2html -a2m_in tmp.a2m -retrieve > T0128.t2k.pa.html SAM: /projects/compbio/bin/alpha/prettyalign v3.2 (July 31, 2000) compiled 08/11/00_16:27:51 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 T0128.t2k.a2m.gz >> tmp.script echo PrintPrediction T0128.t2k.2d 3670-4530-6947 >> tmp.script echo PrintPredictionFasta T0128.t2k.2d.seq >> tmp.script echo PrintRDB T0128.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 T0128.t2k.a2m.gz # Using SequenceWeight EntropyWeight(1.4, 10) # Using regularizer /projects/compbio/lib/recode3.20comp for sequence weight regularizer. # T0128.t2k.a2m with 535 sequences, total weight= 38.6986 avg weight= 0.0723337 46 iterations Regularizing alignment for T0128.t2k.a2m.gz # After reading T0128.t2k.a2m.gz, have 222 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 T0128.t2k.2d # command:# Network initialization done # Printing prediction in FASTA format to T0128.t2k.2d.seq # command:# Network initialization done # Printing prediction to T0128.t2k.2d.rdb # command:rm tmp.script gunzip -f T0128.t2k.a2m.gz /projects/compbio/usr/karplus/src/phytree/phytree -o -i -r flat T0128.t2k T0128.t2k.a2m /projects/compbio/lib/recode2.20comp Prior library /projects/compbio/lib/recode2.20comp read. Reading alignment file T0128.t2k.a2m (535 sequences) as A2M alignment. Alignment T0128.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 One dot for every 10 merges: ..................................................... Merges done. CPU = user 0:21:14 system 0: 0: 1 gzip -f T0128.t2k.a2m gzip -f T0128.t2k_sorted.a2m rm T0128.t2k.phytrace T0128.t2k.tree_weight prettyalign T0128.t2k_sorted.a2m.gz -m8 -i -n -L3333 > T0128.t2k_sorted.pa SAM: prettyalign v3.2 (July 31, 2000) compiled 08/11/00_16:27:51 /projects/compbio/usr/karplus/src/phytree/dg T0128.t2k.tree /projects/compbio/usr/karplus/src/phytree/dtree T0128.t2k.tree /projects/compbio/usr/karplus/src/phytree/dtree: 16336 Memory fault make[1]: Leaving directory `/export/projects/compbio/experiments/casp4/t128' echo "WeightingRegularizer /projects/compbio/lib/recode3.20comp" > tmp.script echo "SequenceWeight EntropyWeight 0.5 10.0" >> tmp.script echo "ReadA2M T0128.t2k.a2m.gz" >> tmp.script echo "SavingsRegularizer /projects/compbio/lib/recode3.20comp" >> tmp.script echo "PrintSavings T0128.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 T0128.t2k.a2m.gz % using regularizer recode3.20comp % T0128.t2k.a2m.gz with 535 sequences, total weight= 3.74441 avg weight= 0.00699889 After reading T0128.t2k.a2m.gz, have 222 columns, totalling 637.686 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 222 columns to T0128.t2k.saves # command:rm tmp.script if test -e T0128.t2k-logo.labels ;\ then makelogo -b 60 -t T0128.t2k -o T0128.t2k-logo.eps -l T0128.t2k-logo.labels T0128.t2k.saves ;\ else makelogo -b 60 -t T0128.t2k -o T0128.t2k-logo.eps -lf T0128.t2k.2d.seq T0128.t2k.saves ;\ fi /projects/compbio/experiments/models.97/scripts99/score-target -alignment T0128.t2k.a2m.gz -build w0.5 \ -db /projects/compbio/data/pdb/all-protein \ -Emax 40.0 \ -viterbi 0 -sw_score 2 \ -output T0128.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-31397 @@@@ mkdir /projects/compbio/tmp/score-target-31397/for-scorescript @@@@ /projects/compbio/experiments/models.97/scripts99/score-mods-vs-seqs.pl -root_tmp /projects/compbio/tmp/score-target-31397/for-scorescript -a2msfile /projects/compbio/tmp/score-target-31397/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-31397/score-target.op @@@@ chgrp protein /projects/compbio/tmp/score-target-31397/for-scorescript/split-tmp @@@@ chgrp protein /projects/compbio/tmp/score-target-31397/for-scorescript/split-tmp @@@@ /projects/compbio/bin/scripts/w0.5 /export/projects/compbio/experiments/casp4/t128/T0128.t2k.a2m.gz /projects/compbio/tmp/score-target-31397/for-scorescript/split-tmp/beta-32552/T0128.t2k.mod Reading /projects/compbio/bin/scripts/sam-t99.conf @@@@ chgrp protein /projects/compbio/tmp/tmp-build-weighted-model-beta-28666 @@@@ uniqueseq /projects/compbio/tmp/tmp-build-weighted-model-beta-28666/tmp -alignfile /export/projects/compbio/experiments/casp4/t128/T0128.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/t128/T0128.t2k.a2m.gz (535 sequences, 222 columns) as A2M alignment. Writing sequence output to /projects/compbio/tmp/tmp-build-weighted-model-beta-28666/tmp.a2m. Dropping 16 (of 535) duplicate sequences with differing IDs Dropping 379 (of 535) sequences with > 80.0% id 140 sequences left after dropping 395 of 535 sequences. @@@@ modelfromalign /projects/compbio/tmp/tmp-build-weighted-model-beta-28666/tmp -alignfile /projects/compbio/tmp/tmp-build-weighted-model-beta-28666/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-28666/tmp.a2m (140 sequences, 222 columns) as A2M alignment. @@@@ cp -f /projects/compbio/tmp/tmp-build-weighted-model-beta-28666/tmp.mod /projects/compbio/tmp/score-target-31397/for-scorescript/split-tmp/beta-32552/T0128.t2k.mod @@@@ rm -rf /projects/compbio/tmp/tmp-build-weighted-model-beta-28666 80.0% id 140 sequences left after dropping 395 of 535 sequences. @@@@ modelfromalign /projects/compbio/tmp/tmp-build-weighted-model-beta-16210/tmp -alignfile /projects/compbio/tmp/tmp-build-weighted-model-beta-16210/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-16210/tmp.a2m (140 sequences, 222 columns) as A2M alignment. @@@@ cp -f /projects/compbio/tmp/tmp-build-weighted-model-beta-16210/tmp.mod T0128.t2k-w0.5.mod @@@@ rm -rf /projects/compbio/tmp/tmp-build-weighted-model-beta-16210