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