-Clazz.declarePackage ("jalview.analysis");
-Clazz.load (null, "jalview.analysis.Conservation", ["jalview.datamodel.Annotation", "$.Sequence", "jalview.schemes.ResidueProperties", "jalview.util.Comparison", "java.awt.Color", "java.lang.Character", "$.Double", "$.StringBuffer", "java.util.Hashtable", "$.Vector"], function () {
-c$ = Clazz.decorateAsClass (function () {
-this.sequences = null;
-this.start = 0;
-this.end = 0;
-this.seqNums = null;
-this.maxLength = 0;
-this.seqNumsChanged = false;
-this.total = null;
-this.canonicaliseAa = true;
-this.quality = null;
-this.qualityRange = null;
-this.consString = "";
-this.consSequence = null;
-this.propHash = null;
-this.threshold = 0;
-this.name = "";
-this.cons2 = null;
-this.consSymbs = null;
-Clazz.instantialize (this, arguments);
-}, jalview.analysis, "Conservation");
-Clazz.prepareFields (c$, function () {
-this.qualityRange = new Array (2);
-});
-Clazz.makeConstructor (c$,
-function (name, propHash, threshold, sequences, start, end) {
-this.name = name;
-this.propHash = propHash;
-this.threshold = threshold;
-this.start = start;
-this.end = end;
-this.maxLength = end - start + 1;
-var s;
-var sSize = sequences.size ();
-var sarray = new Array (sSize);
-this.sequences = sarray;
-try {
-for (s = 0; s < sSize; s++) {
-sarray[s] = sequences.get (s);
-if (sarray[s].getLength () > this.maxLength) {
-this.maxLength = sarray[s].getLength ();
-}}
-} catch (ex) {
-if (Clazz.exceptionOf (ex, ArrayIndexOutOfBoundsException)) {
-this.sequences = new Array (0);
-this.maxLength = 0;
-} else {
-throw ex;
-}
-}
-}, "~S,java.util.Hashtable,~N,java.util.List,~N,~N");
-Clazz.defineMethod (c$, "calcSeqNum",
-($fz = function (i) {
-var sq = null;
-var sqnum = null;
-var sSize = this.sequences.length;
-if ((i > -1) && (i < sSize)) {
-sq = this.sequences[i].getSequenceAsString ();
-if (this.seqNums.size () <= i) {
-this.seqNums.addElement ( Clazz.newIntArray (sq.length + 1, 0));
-}if (sq.hashCode () != (this.seqNums.elementAt (i))[0]) {
-var j;
-var len;
-this.seqNumsChanged = true;
-len = sq.length;
-if (this.maxLength < len) {
-this.maxLength = len;
-}sqnum = Clazz.newIntArray (len + 1, 0);
-sqnum[0] = sq.hashCode ();
-for (j = 1; j <= len; j++) {
-sqnum[j] = jalview.schemes.ResidueProperties.aaIndex[sq.charCodeAt (j - 1)];
-}
-this.seqNums.setElementAt (sqnum, i);
-} else {
-System.out.println ("SEQUENCE HAS BEEN DELETED!!!");
-}} else {
-System.err.println ("ERROR: calcSeqNum called with out of range sequence index for Alignment\n");
-}}, $fz.isPrivate = true, $fz), "~N");
-Clazz.defineMethod (c$, "calculate",
-function () {
-var resultHash;
-var ht;
-var thresh;
-var j;
-var jSize = this.sequences.length;
-var values;
-var type;
-var res = null;
-var c;
-var enumeration2;
-this.total = new Array (this.maxLength);
-for (var i = this.start; i <= this.end; i++) {
-values = Clazz.newIntArray (255, 0);
-for (j = 0; j < jSize; j++) {
-if (this.sequences[j].getLength () > i) {
-c = this.sequences[j].getCharAt (i);
-if (this.canonicaliseAa) {
-c = String.fromCharCode (jalview.schemes.ResidueProperties.aaIndex[this.sequences[j].getCharAt (i).charCodeAt (0)]);
-if (c.charCodeAt (0) > 20) {
-c = '-';
-} else {
-c = jalview.schemes.ResidueProperties.aa[c.charCodeAt (0)].charAt (0);
-}} else {
-if (c == '.' || c == ' ') {
-c = '-';
-}if (!this.canonicaliseAa && 'a' <= c && c <= 'z') {
-c = String.fromCharCode (c.charCodeAt (0) - (32));
-}}values[c.charCodeAt (0)]++;
-} else {
-values['-'.charCodeAt (0)]++;
-}}
-thresh = Clazz.doubleToInt ((this.threshold * (jSize)) / 100);
-resultHash = new java.util.Hashtable ();
-for (var v = '-'; v < 'Z'; v = String.fromCharCode (v.charCodeAt (0) + 1)) {
-if (values[v.charCodeAt (0)] > thresh) {
-res = String.valueOf (v);
-enumeration2 = this.propHash.keys ();
-while (enumeration2.hasMoreElements ()) {
-type = enumeration2.nextElement ();
-ht = this.propHash.get (type);
-if (!resultHash.containsKey (type)) {
-if (ht.containsKey (res)) {
-resultHash.put (type, ht.get (res));
-} else {
-resultHash.put (type, ht.get ("-"));
-}} else if ((resultHash.get (type)).equals (ht.get (res)) == false) {
-resultHash.put (type, new Integer (-1));
-}}
-}}
-if (this.total.length > 0) {
-this.total[i - this.start] = resultHash;
-}}
-});
-Clazz.defineMethod (c$, "countConsNGaps",
-function (j) {
-var count = 0;
-var cons = 0;
-var nres = 0;
-var r = Clazz.newIntArray (2, 0);
-var f = '$';
-var i;
-var iSize = this.sequences.length;
-var c;
-for (i = 0; i < iSize; i++) {
-if (j >= this.sequences[i].getLength ()) {
-count++;
-continue;
-}c = this.sequences[i].getCharAt (j);
-if (jalview.util.Comparison.isGap ((c))) {
-count++;
-} else {
-nres++;
-if (nres == 1) {
-f = c;
-cons++;
-} else if (f == c) {
-cons++;
-}}}
-r[0] = (nres == cons) ? 1 : 0;
-r[1] = count;
-return r;
-}, "~N");
-Clazz.defineMethod (c$, "verdict",
-function (consflag, percentageGaps) {
-var consString = new StringBuffer ();
-var type;
-var result;
-var gapcons;
-var totGaps;
-var count;
-var pgaps;
-var resultHash;
-var enumeration;
-for (var i = 0; i < this.start; i++) {
-consString.append ('-');
-}
-this.consSymbs = new Array (this.end - this.start + 1);
-for (var i = this.start; i <= this.end; i++) {
-gapcons = this.countConsNGaps (i);
-totGaps = gapcons[1];
-pgaps = (totGaps * 100) / this.sequences.length;
-this.consSymbs[i - this.start] = String.instantialize ();
-if (percentageGaps > pgaps) {
-resultHash = this.total[i - this.start];
-count = 0;
-enumeration = resultHash.keys ();
-while (enumeration.hasMoreElements ()) {
-type = enumeration.nextElement ();
-result = resultHash.get (type);
-if (consflag) {
-if (result.intValue () == 1) {
-this.consSymbs[i - this.start] = type + " " + this.consSymbs[i - this.start];
-count++;
-}} else {
-if (result.intValue () != -1) {
-{
-if (result.intValue () == 0) {
-this.consSymbs[i - this.start] = this.consSymbs[i - this.start] + " !" + type;
-} else {
-this.consSymbs[i - this.start] = type + " " + this.consSymbs[i - this.start];
-}}count++;
-}}}
-if (count < 10) {
-consString.append ("" + count);
-} else {
-consString.append ((gapcons[0] == 1) ? "*" : "+");
-}} else {
-consString.append ('-');
-}}
-this.consSequence = new jalview.datamodel.Sequence (this.name, consString.toString (), this.start, this.end);
-}, "~B,~N");
-Clazz.defineMethod (c$, "getConsSequence",
-function () {
-return this.consSequence;
-});
-Clazz.defineMethod (c$, "findQuality",
-function () {
-this.findQuality (0, this.maxLength - 1);
-});
-Clazz.defineMethod (c$, "percentIdentity2",
-($fz = function () {
-this.seqNums = new java.util.Vector ();
-var i = 0;
-var iSize = this.sequences.length;
-for (i = 0; i < iSize; i++) {
-this.calcSeqNum (i);
-}
-if ((this.cons2 == null) || this.seqNumsChanged) {
-this.cons2 = Clazz.newIntArray (this.maxLength, 24, 0);
-for (var j = 0; j < 24; j++) {
-for (i = 0; i < this.maxLength; i++) {
-this.cons2[i][j] = 0;
-}
-}
-var sqnum;
-var j = 0;
-while (j < this.sequences.length) {
-sqnum = this.seqNums.elementAt (j);
-for (i = 1; i < sqnum.length; i++) {
-this.cons2[i - 1][sqnum[i]]++;
-}
-for (i = sqnum.length - 1; i < this.maxLength; i++) {
-this.cons2[i][23]++;
-}
-j++;
-}
-}}, $fz.isPrivate = true, $fz));
-Clazz.defineMethod (c$, "findQuality",
-function (start, end) {
-this.quality = new java.util.Vector ();
-var max = -10000;
-var BLOSUM62 = jalview.schemes.ResidueProperties.getBLOSUM62 ();
-this.percentIdentity2 ();
-var size = this.seqNums.size ();
-var lengths = Clazz.newIntArray (size, 0);
-var tot;
-var bigtot;
-var sr;
-var tmp;
-var x;
-var xx;
-var l;
-var j;
-var i;
-var ii;
-var i2;
-var k;
-var seqNum;
-for (l = 0; l < size; l++) {
-lengths[l] = (this.seqNums.elementAt (l)).length - 1;
-}
-for (j = start; j <= end; j++) {
-bigtot = 0;
-x = Clazz.newDoubleArray (24, 0);
-for (ii = 0; ii < 24; ii++) {
-x[ii] = 0;
-for (i2 = 0; i2 < 24; i2++) {
-x[ii] += ((this.cons2[j][i2] * BLOSUM62[ii][i2]) + 4);
-}
-x[ii] /= size;
-}
-for (k = 0; k < size; k++) {
-tot = 0;
-xx = Clazz.newDoubleArray (24, 0);
-seqNum = (j < lengths[k]) ? (this.seqNums.elementAt (k))[j + 1] : 23;
-for (i = 0; i < 23; i++) {
-sr = 0;
-sr = BLOSUM62[i][seqNum] + 4;
-xx[i] = x[i] - sr;
-tot += (xx[i] * xx[i]);
-}
-bigtot += Math.sqrt (tot);
-}
-if (max < bigtot) {
-max = bigtot;
-}this.quality.addElement ( new Double (bigtot));
-}
-var newmax = -10000;
-for (j = start; j <= end; j++) {
-tmp = (this.quality.elementAt (j)).doubleValue ();
-tmp = ((max - tmp) * (size - this.cons2[j][23])) / size;
-this.quality.setElementAt ( new Double (tmp), j);
-if (tmp > newmax) {
-newmax = tmp;
-}}
-this.qualityRange[0] = new Double (0);
-this.qualityRange[1] = new Double (newmax);
-}, "~N,~N");
-Clazz.defineMethod (c$, "completeAnnotations",
-function (conservation, quality2, istart, alWidth) {
-var sequence = this.getConsSequence ().getSequence ();
-var minR;
-var minG;
-var minB;
-var maxR;
-var maxG;
-var maxB;
-minR = 0.3;
-minG = 0.0;
-minB = 0;
-maxR = 1.0 - minR;
-maxG = 0.9 - minG;
-maxB = 0 - minB;
-var min = 0;
-var max = 11;
-var qmin = 0;
-var qmax = 0;
-var c;
-if (conservation.annotations != null && conservation.annotations.length < alWidth) {
-conservation.annotations = new Array (alWidth);
-}if (quality2 != null) {
-quality2.graphMax = this.qualityRange[1].floatValue ();
-if (quality2.annotations != null && quality2.annotations.length < alWidth) {
-quality2.annotations = new Array (alWidth);
-}qmin = this.qualityRange[0].floatValue ();
-qmax = this.qualityRange[1].floatValue ();
-}for (var i = 0; i < alWidth; i++) {
-var value = 0;
-c = sequence[i];
-if (Character.isDigit (c)) {
-value = c.charCodeAt (0) - 48;
-} else if (c == '*') {
-value = 11;
-} else if (c == '+') {
-value = 10;
-}var vprop = value - min;
-vprop /= max;
-conservation.annotations[i] = new jalview.datamodel.Annotation (String.valueOf (c), this.consSymbs[i - this.start], ' ', value, new java.awt.Color (minR + (maxR * vprop), minG + (maxG * vprop), minB + (maxB * vprop)));
-if (quality2 != null) {
-value = (this.quality.elementAt (i)).floatValue ();
-vprop = value - qmin;
-vprop /= qmax;
-quality2.annotations[i] = new jalview.datamodel.Annotation (" ", String.valueOf (value), ' ', value, new java.awt.Color (minR + (maxR * vprop), minG + (maxG * vprop), minB + (maxB * vprop)));
-}}
-}, "jalview.datamodel.AlignmentAnnotation,jalview.datamodel.AlignmentAnnotation,~N,~N");
-c$.calculateConservation = Clazz.defineMethod (c$, "calculateConservation",
-function (name, consHash, threshold, seqs, start, end, posOrNeg, consPercGaps, calcQuality) {
-var cons = new jalview.analysis.Conservation (name, consHash, threshold, seqs, start, end);
-return jalview.analysis.Conservation.calculateConservation (cons, posOrNeg, consPercGaps, calcQuality);
-}, "~S,java.util.Hashtable,~N,java.util.List,~N,~N,~B,~N,~B");
-c$.calculateConservation = Clazz.defineMethod (c$, "calculateConservation",
-function (cons, b, consPercGaps, calcQuality) {
-cons.calculate ();
-cons.verdict (b, consPercGaps);
-if (calcQuality) {
-cons.findQuality ();
-}return cons;
-}, "jalview.analysis.Conservation,~B,~N,~B");
-});
+Clazz.declarePackage ("jalview.analysis");\r
+Clazz.load (null, "jalview.analysis.Conservation", ["jalview.datamodel.Annotation", "$.Sequence", "jalview.schemes.ResidueProperties", "jalview.util.Comparison", "java.awt.Color", "java.lang.Character", "$.Double", "$.StringBuffer", "java.util.Hashtable", "$.Vector"], function () {\r
+c$ = Clazz.decorateAsClass (function () {\r
+this.sequences = null;\r
+this.start = 0;\r
+this.end = 0;\r
+this.seqNums = null;\r
+this.maxLength = 0;\r
+this.seqNumsChanged = false;\r
+this.total = null;\r
+this.canonicaliseAa = true;\r
+this.quality = null;\r
+this.qualityRange = null;\r
+this.consString = "";\r
+this.consSequence = null;\r
+this.propHash = null;\r
+this.threshold = 0;\r
+this.name = "";\r
+this.cons2 = null;\r
+this.consSymbs = null;\r
+Clazz.instantialize (this, arguments);\r
+}, jalview.analysis, "Conservation");\r
+Clazz.prepareFields (c$, function () {\r
+this.qualityRange = new Array (2);\r
+});\r
+Clazz.makeConstructor (c$, \r
+function (name, propHash, threshold, sequences, start, end) {\r
+this.name = name;\r
+this.propHash = propHash;\r
+this.threshold = threshold;\r
+this.start = start;\r
+this.end = end;\r
+this.maxLength = end - start + 1;\r
+var s;\r
+var sSize = sequences.size ();\r
+var sarray = new Array (sSize);\r
+this.sequences = sarray;\r
+try {\r
+for (s = 0; s < sSize; s++) {\r
+sarray[s] = sequences.get (s);\r
+if (sarray[s].getLength () > this.maxLength) {\r
+this.maxLength = sarray[s].getLength ();\r
+}}\r
+} catch (ex) {\r
+if (Clazz.exceptionOf (ex, ArrayIndexOutOfBoundsException)) {\r
+this.sequences = new Array (0);\r
+this.maxLength = 0;\r
+} else {\r
+throw ex;\r
+}\r
+}\r
+}, "~S,java.util.Hashtable,~N,java.util.List,~N,~N");\r
+Clazz.defineMethod (c$, "calcSeqNum", \r
+($fz = function (i) {\r
+var sq = null;\r
+var sqnum = null;\r
+var sSize = this.sequences.length;\r
+if ((i > -1) && (i < sSize)) {\r
+sq = this.sequences[i].getSequenceAsString ();\r
+if (this.seqNums.size () <= i) {\r
+this.seqNums.addElement ( Clazz.newIntArray (sq.length + 1, 0));\r
+}if (sq.hashCode () != (this.seqNums.elementAt (i))[0]) {\r
+var j;\r
+var len;\r
+this.seqNumsChanged = true;\r
+len = sq.length;\r
+if (this.maxLength < len) {\r
+this.maxLength = len;\r
+}sqnum = Clazz.newIntArray (len + 1, 0);\r
+sqnum[0] = sq.hashCode ();\r
+for (j = 1; j <= len; j++) {\r
+sqnum[j] = jalview.schemes.ResidueProperties.aaIndex[sq.charCodeAt (j - 1)];\r
+}\r
+this.seqNums.setElementAt (sqnum, i);\r
+} else {\r
+System.out.println ("SEQUENCE HAS BEEN DELETED!!!");\r
+}} else {\r
+System.err.println ("ERROR: calcSeqNum called with out of range sequence index for Alignment\n");\r
+}}, $fz.isPrivate = true, $fz), "~N");\r
+Clazz.defineMethod (c$, "calculate", \r
+function () {\r
+var resultHash;\r
+var ht;\r
+var thresh;\r
+var j;\r
+var jSize = this.sequences.length;\r
+var values;\r
+var type;\r
+var res = null;\r
+var c;\r
+var enumeration2;\r
+this.total = new Array (this.maxLength);\r
+for (var i = this.start; i <= this.end; i++) {\r
+values = Clazz.newIntArray (255, 0);\r
+for (j = 0; j < jSize; j++) {\r
+if (this.sequences[j].getLength () > i) {\r
+c = this.sequences[j].getCharAt (i);\r
+if (this.canonicaliseAa) {\r
+c = String.fromCharCode (jalview.schemes.ResidueProperties.aaIndex[this.sequences[j].getCharAt (i).charCodeAt (0)]);\r
+if (c.charCodeAt (0) > 20) {\r
+c = '-';\r
+} else {\r
+c = jalview.schemes.ResidueProperties.aa[c.charCodeAt (0)].charAt (0);\r
+}} else {\r
+if (c == '.' || c == ' ') {\r
+c = '-';\r
+}if (!this.canonicaliseAa && 'a' <= c && c <= 'z') {\r
+c = String.fromCharCode (c.charCodeAt (0) - (32));\r
+}}values[c.charCodeAt (0)]++;\r
+} else {\r
+values['-'.charCodeAt (0)]++;\r
+}}\r
+thresh = Clazz.doubleToInt ((this.threshold * (jSize)) / 100);\r
+resultHash = new java.util.Hashtable ();\r
+for (var v = '-'; v < 'Z'; v = String.fromCharCode (v.charCodeAt (0) + 1)) {\r
+if (values[v.charCodeAt (0)] > thresh) {\r
+res = String.valueOf (v);\r
+enumeration2 = this.propHash.keys ();\r
+while (enumeration2.hasMoreElements ()) {\r
+type = enumeration2.nextElement ();\r
+ht = this.propHash.get (type);\r
+if (!resultHash.containsKey (type)) {\r
+if (ht.containsKey (res)) {\r
+resultHash.put (type, ht.get (res));\r
+} else {\r
+resultHash.put (type, ht.get ("-"));\r
+}} else if ((resultHash.get (type)).equals (ht.get (res)) == false) {\r
+resultHash.put (type, new Integer (-1));\r
+}}\r
+}}\r
+if (this.total.length > 0) {\r
+this.total[i - this.start] = resultHash;\r
+}}\r
+});\r
+Clazz.defineMethod (c$, "countConsNGaps", \r
+function (j) {\r
+var count = 0;\r
+var cons = 0;\r
+var nres = 0;\r
+var r = Clazz.newIntArray (2, 0);\r
+var f = '$';\r
+var i;\r
+var iSize = this.sequences.length;\r
+var c;\r
+for (i = 0; i < iSize; i++) {\r
+if (j >= this.sequences[i].getLength ()) {\r
+count++;\r
+continue;\r
+}c = this.sequences[i].getCharAt (j);\r
+if (jalview.util.Comparison.isGap ((c))) {\r
+count++;\r
+} else {\r
+nres++;\r
+if (nres == 1) {\r
+f = c;\r
+cons++;\r
+} else if (f == c) {\r
+cons++;\r
+}}}\r
+r[0] = (nres == cons) ? 1 : 0;\r
+r[1] = count;\r
+return r;\r
+}, "~N");\r
+Clazz.defineMethod (c$, "verdict", \r
+function (consflag, percentageGaps) {\r
+var consString = new StringBuffer ();\r
+var type;\r
+var result;\r
+var gapcons;\r
+var totGaps;\r
+var count;\r
+var pgaps;\r
+var resultHash;\r
+var enumeration;\r
+for (var i = 0; i < this.start; i++) {\r
+consString.append ('-');\r
+}\r
+this.consSymbs = new Array (this.end - this.start + 1);\r
+for (var i = this.start; i <= this.end; i++) {\r
+gapcons = this.countConsNGaps (i);\r
+totGaps = gapcons[1];\r
+pgaps = (totGaps * 100) / this.sequences.length;\r
+this.consSymbs[i - this.start] = String.instantialize ();\r
+if (percentageGaps > pgaps) {\r
+resultHash = this.total[i - this.start];\r
+count = 0;\r
+enumeration = resultHash.keys ();\r
+while (enumeration.hasMoreElements ()) {\r
+type = enumeration.nextElement ();\r
+result = resultHash.get (type);\r
+if (consflag) {\r
+if (result.intValue () == 1) {\r
+this.consSymbs[i - this.start] = type + " " + this.consSymbs[i - this.start];\r
+count++;\r
+}} else {\r
+if (result.intValue () != -1) {\r
+{\r
+if (result.intValue () == 0) {\r
+this.consSymbs[i - this.start] = this.consSymbs[i - this.start] + " !" + type;\r
+} else {\r
+this.consSymbs[i - this.start] = type + " " + this.consSymbs[i - this.start];\r
+}}count++;\r
+}}}\r
+if (count < 10) {\r
+consString.append ("" + count);\r
+} else {\r
+consString.append ((gapcons[0] == 1) ? "*" : "+");\r
+}} else {\r
+consString.append ('-');\r
+}}\r
+this.consSequence = new jalview.datamodel.Sequence (this.name, consString.toString (), this.start, this.end);\r
+}, "~B,~N");\r
+Clazz.defineMethod (c$, "getConsSequence", \r
+function () {\r
+return this.consSequence;\r
+});\r
+Clazz.defineMethod (c$, "findQuality", \r
+function () {\r
+this.findQuality (0, this.maxLength - 1);\r
+});\r
+Clazz.defineMethod (c$, "percentIdentity2", \r
+($fz = function () {\r
+this.seqNums = new java.util.Vector ();\r
+var i = 0;\r
+var iSize = this.sequences.length;\r
+for (i = 0; i < iSize; i++) {\r
+this.calcSeqNum (i);\r
+}\r
+if ((this.cons2 == null) || this.seqNumsChanged) {\r
+this.cons2 = Clazz.newIntArray (this.maxLength, 24, 0);\r
+for (var j = 0; j < 24; j++) {\r
+for (i = 0; i < this.maxLength; i++) {\r
+this.cons2[i][j] = 0;\r
+}\r
+}\r
+var sqnum;\r
+var j = 0;\r
+while (j < this.sequences.length) {\r
+sqnum = this.seqNums.elementAt (j);\r
+for (i = 1; i < sqnum.length; i++) {\r
+this.cons2[i - 1][sqnum[i]]++;\r
+}\r
+for (i = sqnum.length - 1; i < this.maxLength; i++) {\r
+this.cons2[i][23]++;\r
+}\r
+j++;\r
+}\r
+}}, $fz.isPrivate = true, $fz));\r
+Clazz.defineMethod (c$, "findQuality", \r
+function (start, end) {\r
+this.quality = new java.util.Vector ();\r
+var max = -10000;\r
+var BLOSUM62 = jalview.schemes.ResidueProperties.getBLOSUM62 ();\r
+this.percentIdentity2 ();\r
+var size = this.seqNums.size ();\r
+var lengths = Clazz.newIntArray (size, 0);\r
+var tot;\r
+var bigtot;\r
+var sr;\r
+var tmp;\r
+var x;\r
+var xx;\r
+var l;\r
+var j;\r
+var i;\r
+var ii;\r
+var i2;\r
+var k;\r
+var seqNum;\r
+for (l = 0; l < size; l++) {\r
+lengths[l] = (this.seqNums.elementAt (l)).length - 1;\r
+}\r
+for (j = start; j <= end; j++) {\r
+bigtot = 0;\r
+x = Clazz.newDoubleArray (24, 0);\r
+for (ii = 0; ii < 24; ii++) {\r
+x[ii] = 0;\r
+for (i2 = 0; i2 < 24; i2++) {\r
+x[ii] += ((this.cons2[j][i2] * BLOSUM62[ii][i2]) + 4);\r
+}\r
+x[ii] /= size;\r
+}\r
+for (k = 0; k < size; k++) {\r
+tot = 0;\r
+xx = Clazz.newDoubleArray (24, 0);\r
+seqNum = (j < lengths[k]) ? (this.seqNums.elementAt (k))[j + 1] : 23;\r
+for (i = 0; i < 23; i++) {\r
+sr = 0;\r
+sr = BLOSUM62[i][seqNum] + 4;\r
+xx[i] = x[i] - sr;\r
+tot += (xx[i] * xx[i]);\r
+}\r
+bigtot += Math.sqrt (tot);\r
+}\r
+if (max < bigtot) {\r
+max = bigtot;\r
+}this.quality.addElement ( new Double (bigtot));\r
+}\r
+var newmax = -10000;\r
+for (j = start; j <= end; j++) {\r
+tmp = (this.quality.elementAt (j)).doubleValue ();\r
+tmp = ((max - tmp) * (size - this.cons2[j][23])) / size;\r
+this.quality.setElementAt ( new Double (tmp), j);\r
+if (tmp > newmax) {\r
+newmax = tmp;\r
+}}\r
+this.qualityRange[0] = new Double (0);\r
+this.qualityRange[1] = new Double (newmax);\r
+}, "~N,~N");\r
+Clazz.defineMethod (c$, "completeAnnotations", \r
+function (conservation, quality2, istart, alWidth) {\r
+var sequence = this.getConsSequence ().getSequence ();\r
+var minR;\r
+var minG;\r
+var minB;\r
+var maxR;\r
+var maxG;\r
+var maxB;\r
+minR = 0.3;\r
+minG = 0.0;\r
+minB = 0;\r
+maxR = 1.0 - minR;\r
+maxG = 0.9 - minG;\r
+maxB = 0 - minB;\r
+var min = 0;\r
+var max = 11;\r
+var qmin = 0;\r
+var qmax = 0;\r
+var c;\r
+if (conservation.annotations != null && conservation.annotations.length < alWidth) {\r
+conservation.annotations = new Array (alWidth);\r
+}if (quality2 != null) {\r
+quality2.graphMax = this.qualityRange[1].floatValue ();\r
+if (quality2.annotations != null && quality2.annotations.length < alWidth) {\r
+quality2.annotations = new Array (alWidth);\r
+}qmin = this.qualityRange[0].floatValue ();\r
+qmax = this.qualityRange[1].floatValue ();\r
+}for (var i = 0; i < alWidth; i++) {\r
+var value = 0;\r
+c = sequence[i];\r
+if (Character.isDigit (c)) {\r
+value = c.charCodeAt (0) - 48;\r
+} else if (c == '*') {\r
+value = 11;\r
+} else if (c == '+') {\r
+value = 10;\r
+}var vprop = value - min;\r
+vprop /= max;\r
+conservation.annotations[i] = new jalview.datamodel.Annotation (String.valueOf (c), this.consSymbs[i - this.start], ' ', value, new java.awt.Color (minR + (maxR * vprop), minG + (maxG * vprop), minB + (maxB * vprop)));\r
+if (quality2 != null) {\r
+value = (this.quality.elementAt (i)).floatValue ();\r
+vprop = value - qmin;\r
+vprop /= qmax;\r
+quality2.annotations[i] = new jalview.datamodel.Annotation (" ", String.valueOf (value), ' ', value, new java.awt.Color (minR + (maxR * vprop), minG + (maxG * vprop), minB + (maxB * vprop)));\r
+}}\r
+}, "jalview.datamodel.AlignmentAnnotation,jalview.datamodel.AlignmentAnnotation,~N,~N");\r
+c$.calculateConservation = Clazz.defineMethod (c$, "calculateConservation", \r
+function (name, consHash, threshold, seqs, start, end, posOrNeg, consPercGaps, calcQuality) {\r
+var cons = new jalview.analysis.Conservation (name, consHash, threshold, seqs, start, end);\r
+return jalview.analysis.Conservation.calculateConservation (cons, posOrNeg, consPercGaps, calcQuality);\r
+}, "~S,java.util.Hashtable,~N,java.util.List,~N,~N,~B,~N,~B");\r
+c$.calculateConservation = Clazz.defineMethod (c$, "calculateConservation", \r
+function (cons, b, consPercGaps, calcQuality) {\r
+cons.calculate ();\r
+cons.verdict (b, consPercGaps);\r
+if (calcQuality) {\r
+cons.findQuality ();\r
+}return cons;\r
+}, "jalview.analysis.Conservation,~B,~N,~B");\r
+});\r
git://source.jalview.org / jalviewjs.git / blobdiff