echo "WeightingRegularizer /projects/compbio/lib/recode3.20comp" > tmp.script echo "SequenceWeight EntropyWeight 0.5 10.0" >> tmp.script echo "ReadA2M T0111.t2k.a2m.gz" >> tmp.script echo "SavingsRegularizer /projects/compbio/lib/recode3.20comp" >> tmp.script echo "PrintSavings T0111.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 T0111.t2k.a2m.gz % using regularizer recode3.20comp % T0111.t2k.a2m.gz with 198 sequences, total weight= 2.91204 avg weight= 0.0147073 After reading T0111.t2k.a2m.gz, have 431 columns, totalling 1013.43 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 431 columns to T0111.t2k.saves # command:rm tmp.script if test -e T0111.t2k-logo.labels ;\ then makelogo -b 60 -t T0111.t2k -o T0111.t2k-logo.eps -l T0111.t2k-logo.labels T0111.t2k.saves ;\ else makelogo -b 60 -t T0111.t2k -o T0111.t2k-logo.eps -lf T0111.t2k.2d.seq T0111.t2k.saves ;\ fi rm T0111.t2k.saves