-Clazz.declarePackage ("jalview.analysis");\r
-Clazz.load (null, "jalview.analysis.AlignmentUtils", ["jalview.analysis.CodonComparator", "jalview.datamodel.AlignedCodonFrame", "$.Alignment", "$.AlignmentAnnotation", "$.DBRefEntry", "$.DBRefSource", "$.FeatureProperties", "$.SearchResults", "$.Sequence", "jalview.schemes.ResidueProperties", "jalview.util.DBRefUtils", "$.MapList", "$.MappingUtils", "java.lang.StringBuilder", "java.util.ArrayList", "$.Arrays", "$.HashMap", "$.HashSet", "$.LinkedHashMap", "$.LinkedHashSet", "$.TreeMap"], function () {\r
-c$ = Clazz.declareType (jalview.analysis, "AlignmentUtils");\r
-c$.expandContext = Clazz.defineMethod (c$, "expandContext", \r
-function (core, flankSize) {\r
-var sq = new java.util.ArrayList ();\r
-var maxoffset = 0;\r
-for (var s, $s = core.getSequences ().iterator (); $s.hasNext () && ((s = $s.next ()) || true);) {\r
-var newSeq = s.deriveSequence ();\r
-var newSeqStart = newSeq.getStart () - 1;\r
-if (newSeqStart > maxoffset && newSeq.getDatasetSequence ().getStart () < s.getStart ()) {\r
-maxoffset = newSeqStart;\r
-}sq.add (newSeq);\r
-}\r
-if (flankSize > -1) {\r
-maxoffset = Math.min (maxoffset, flankSize);\r
-}for (var s, $s = sq.iterator (); $s.hasNext () && ((s = $s.next ()) || true);) {\r
-var ds = s;\r
-while (ds.getDatasetSequence () != null) {\r
-ds = ds.getDatasetSequence ();\r
-}\r
-var s_end = s.findPosition (s.getStart () + s.getLength ());\r
-var ustream_ds = s.getStart () - ds.getStart ();\r
-var dstream_ds = ds.getEnd () - s_end;\r
-var offset = maxoffset - ustream_ds;\r
-if (flankSize >= 0) {\r
-if (flankSize < ustream_ds) {\r
-offset = maxoffset - flankSize;\r
-ustream_ds = flankSize;\r
-}if (flankSize <= dstream_ds) {\r
-dstream_ds = flankSize - 1;\r
-}}var upstream = String.instantialize (ds.getSequence (s.getStart () - 1 - ustream_ds, s.getStart () - 1)).toLowerCase ().toCharArray ();\r
-var downstream = String.instantialize (ds.getSequence (s_end - 1, s_end + dstream_ds)).toLowerCase ().toCharArray ();\r
-var coreseq = s.getSequence ();\r
-var nseq = Clazz.newCharArray (offset + upstream.length + downstream.length + coreseq.length, '\0');\r
-var c = core.getGapCharacter ();\r
-var p = 0;\r
-for (; p < offset; p++) {\r
-nseq[p] = c;\r
-}\r
-System.arraycopy (upstream, 0, nseq, p, upstream.length);\r
-System.arraycopy (coreseq, 0, nseq, p + upstream.length, coreseq.length);\r
-System.arraycopy (downstream, 0, nseq, p + coreseq.length + upstream.length, downstream.length);\r
-s.setSequence ( String.instantialize (nseq));\r
-s.setStart (s.getStart () - ustream_ds);\r
-s.setEnd (s_end + downstream.length);\r
-}\r
-var newAl = new jalview.datamodel.Alignment (sq.toArray ( new Array (0)));\r
-for (var s, $s = sq.iterator (); $s.hasNext () && ((s = $s.next ()) || true);) {\r
-if (s.getAnnotation () != null) {\r
-for (var aa, $aa = 0, $$aa = s.getAnnotation (); $aa < $$aa.length && ((aa = $$aa[$aa]) || true); $aa++) {\r
-aa.adjustForAlignment ();\r
-newAl.addAnnotation (aa);\r
-}\r
-}}\r
-newAl.setDataset (core.getDataset ());\r
-return newAl;\r
-}, "jalview.datamodel.AlignmentI,~N");\r
-c$.getSequenceIndex = Clazz.defineMethod (c$, "getSequenceIndex", \r
-function (al, seq) {\r
-var result = -1;\r
-var pos = 0;\r
-for (var alSeq, $alSeq = al.getSequences ().iterator (); $alSeq.hasNext () && ((alSeq = $alSeq.next ()) || true);) {\r
-if (alSeq === seq) {\r
-result = pos;\r
-break;\r
-}pos++;\r
-}\r
-return result;\r
-}, "jalview.datamodel.AlignmentI,jalview.datamodel.SequenceI");\r
-c$.getSequencesByName = Clazz.defineMethod (c$, "getSequencesByName", \r
-function (al) {\r
-var theMap = new java.util.LinkedHashMap ();\r
-for (var seq, $seq = al.getSequences ().iterator (); $seq.hasNext () && ((seq = $seq.next ()) || true);) {\r
-var name = seq.getName ();\r
-if (name != null) {\r
-var seqs = theMap.get (name);\r
-if (seqs == null) {\r
-seqs = new java.util.ArrayList ();\r
-theMap.put (name, seqs);\r
-}seqs.add (seq);\r
-}}\r
-return theMap;\r
-}, "jalview.datamodel.AlignmentI");\r
-c$.mapProteinToCdna = Clazz.defineMethod (c$, "mapProteinToCdna", \r
-function (proteinAlignment, cdnaAlignment) {\r
-if (proteinAlignment == null || cdnaAlignment == null) {\r
-return false;\r
-}var mappedDna = new java.util.HashSet ();\r
-var mappedProtein = new java.util.HashSet ();\r
-var mappingPerformed = jalview.analysis.AlignmentUtils.mapProteinToCdna (proteinAlignment, cdnaAlignment, mappedDna, mappedProtein, true);\r
-mappingPerformed = new Boolean (mappingPerformed | jalview.analysis.AlignmentUtils.mapProteinToCdna (proteinAlignment, cdnaAlignment, mappedDna, mappedProtein, false)).valueOf ();\r
-return mappingPerformed;\r
-}, "jalview.datamodel.AlignmentI,jalview.datamodel.AlignmentI");\r
-c$.mapProteinToCdna = Clazz.defineMethod (c$, "mapProteinToCdna", \r
-function (proteinAlignment, cdnaAlignment, mappedDna, mappedProtein, xrefsOnly) {\r
-var mappingPerformed = false;\r
-var thisSeqs = proteinAlignment.getSequences ();\r
-for (var aaSeq, $aaSeq = thisSeqs.iterator (); $aaSeq.hasNext () && ((aaSeq = $aaSeq.next ()) || true);) {\r
-var proteinMapped = false;\r
-var acf = new jalview.datamodel.AlignedCodonFrame ();\r
-for (var cdnaSeq, $cdnaSeq = cdnaAlignment.getSequences ().iterator (); $cdnaSeq.hasNext () && ((cdnaSeq = $cdnaSeq.next ()) || true);) {\r
-if (xrefsOnly && !jalview.analysis.AlignmentUtils.haveCrossRef (aaSeq, cdnaSeq)) {\r
-continue;\r
-}if (!xrefsOnly && (mappedProtein.contains (aaSeq) || mappedDna.contains (cdnaSeq))) {\r
-continue;\r
-}if (!jalview.analysis.AlignmentUtils.mappingExists (proteinAlignment.getCodonFrames (), aaSeq.getDatasetSequence (), cdnaSeq.getDatasetSequence ())) {\r
-var map = jalview.analysis.AlignmentUtils.mapProteinToCdna (aaSeq, cdnaSeq);\r
-if (map != null) {\r
-acf.addMap (cdnaSeq, aaSeq, map);\r
-mappingPerformed = true;\r
-proteinMapped = true;\r
-mappedDna.add (cdnaSeq);\r
-mappedProtein.add (aaSeq);\r
-}}}\r
-if (proteinMapped) {\r
-proteinAlignment.addCodonFrame (acf);\r
-}}\r
-return mappingPerformed;\r
-}, "jalview.datamodel.AlignmentI,jalview.datamodel.AlignmentI,java.util.Set,java.util.Set,~B");\r
-c$.mappingExists = Clazz.defineMethod (c$, "mappingExists", \r
-function (set, aaSeq, cdnaSeq) {\r
-if (set != null) {\r
-for (var acf, $acf = set.iterator (); $acf.hasNext () && ((acf = $acf.next ()) || true);) {\r
-if (cdnaSeq === acf.getDnaForAaSeq (aaSeq)) {\r
-return true;\r
-}}\r
-}return false;\r
-}, "java.util.Set,jalview.datamodel.SequenceI,jalview.datamodel.SequenceI");\r
-c$.mapProteinToCdna = Clazz.defineMethod (c$, "mapProteinToCdna", \r
-function (proteinSeq, cdnaSeq) {\r
-var proteinDataset = proteinSeq.getDatasetSequence ();\r
-var aaSeqChars = proteinDataset != null ? proteinDataset.getSequence () : proteinSeq.getSequence ();\r
-var cdnaDataset = cdnaSeq.getDatasetSequence ();\r
-var cdnaSeqChars = cdnaDataset != null ? cdnaDataset.getSequence () : cdnaSeq.getSequence ();\r
-if (aaSeqChars == null || cdnaSeqChars == null) {\r
-return null;\r
-}var mappedLength = 3 * aaSeqChars.length;\r
-var cdnaLength = cdnaSeqChars.length;\r
-var cdnaStart = 1;\r
-var cdnaEnd = cdnaLength;\r
-var proteinStart = 1;\r
-var proteinEnd = aaSeqChars.length;\r
-if (cdnaLength != mappedLength && cdnaLength > 2) {\r
-var lastCodon = String.valueOf (cdnaSeqChars, cdnaLength - 3, 3).toUpperCase ();\r
-for (var stop, $stop = jalview.schemes.ResidueProperties.STOP.iterator (); $stop.hasNext () && ((stop = $stop.next ()) || true);) {\r
-if (lastCodon.equals (stop)) {\r
-cdnaEnd -= 3;\r
-cdnaLength -= 3;\r
-break;\r
-}}\r
-}if (cdnaLength != mappedLength && cdnaLength > 2 && String.valueOf (cdnaSeqChars, 0, 3).toUpperCase ().equals (jalview.schemes.ResidueProperties.START)) {\r
-cdnaStart += 3;\r
-cdnaLength -= 3;\r
-}if (cdnaLength != mappedLength) {\r
-return null;\r
-}if (!jalview.analysis.AlignmentUtils.translatesAs (cdnaSeqChars, cdnaStart - 1, aaSeqChars)) {\r
-return null;\r
-}var map = new jalview.util.MapList ( Clazz.newIntArray (-1, [cdnaStart, cdnaEnd]), Clazz.newIntArray (-1, [1, proteinEnd]), 3, 1);\r
-return map;\r
-}, "jalview.datamodel.SequenceI,jalview.datamodel.SequenceI");\r
-c$.translatesAs = Clazz.defineMethod (c$, "translatesAs", \r
-function (cdnaSeqChars, cdnaStart, aaSeqChars) {\r
-var aaResidue = 0;\r
-for (var i = cdnaStart; i < cdnaSeqChars.length - 2 && aaResidue < aaSeqChars.length; i += 3, aaResidue++) {\r
-var codon = String.valueOf (cdnaSeqChars, i, 3);\r
-var translated = jalview.schemes.ResidueProperties.codonTranslate (codon);\r
-var aaRes = aaSeqChars[aaResidue];\r
-if ((translated == null || "STOP".equals (translated)) && aaRes == 'X') {\r
-continue;\r
-}if (translated == null || !(aaRes == translated.charAt (0))) {\r
-return false;\r
-}}\r
-return (aaResidue == aaSeqChars.length);\r
-}, "~A,~N,~A");\r
-c$.alignSequenceAs = Clazz.defineMethod (c$, "alignSequenceAs", \r
-function (seq, al, gap, preserveMappedGaps, preserveUnmappedGaps) {\r
-var mappings = al.getCodonFrame (seq);\r
-if (mappings == null || mappings.isEmpty ()) {\r
-return false;\r
-}var alignFrom = null;\r
-var mapping = null;\r
-for (var mp, $mp = mappings.iterator (); $mp.hasNext () && ((mp = $mp.next ()) || true);) {\r
-alignFrom = mp.findAlignedSequence (seq.getDatasetSequence (), al);\r
-if (alignFrom != null) {\r
-mapping = mp;\r
-break;\r
-}}\r
-if (alignFrom == null) {\r
-return false;\r
-}jalview.analysis.AlignmentUtils.alignSequenceAs (seq, alignFrom, mapping, gap, al.getGapCharacter (), preserveMappedGaps, preserveUnmappedGaps);\r
-return true;\r
-}, "jalview.datamodel.SequenceI,jalview.datamodel.AlignmentI,~S,~B,~B");\r
-c$.alignSequenceAs = Clazz.defineMethod (c$, "alignSequenceAs", \r
-function (alignTo, alignFrom, mapping, myGap, sourceGap, preserveMappedGaps, preserveUnmappedGaps) {\r
-var thisSeq = alignTo.getSequence ();\r
-var thatAligned = alignFrom.getSequence ();\r
-var thisAligned = new StringBuilder (2 * thisSeq.length);\r
-var thisSeqPos = 0;\r
-var sourceDsPos = 0;\r
-var basesWritten = 0;\r
-var myGapChar = myGap.charAt (0);\r
-var ratio = myGap.length;\r
-var sourceGapMappedLength = 0;\r
-var inExon = false;\r
-for (var sourceChar, $sourceChar = 0, $$sourceChar = thatAligned; $sourceChar < $$sourceChar.length && ((sourceChar = $$sourceChar[$sourceChar]) || true); $sourceChar++) {\r
-if (sourceChar == sourceGap) {\r
-sourceGapMappedLength += ratio;\r
-continue;\r
-}sourceDsPos++;\r
-var mappedPos = mapping.getMappedRegion (alignTo, alignFrom, sourceDsPos);\r
-if (mappedPos == null) {\r
-System.err.println ("Can't align: no codon mapping to residue " + sourceDsPos + "(" + sourceChar + ")");\r
-return;\r
-}var mappedCodonStart = mappedPos[0];\r
-var mappedCodonEnd = mappedPos[mappedPos.length - 1];\r
-var trailingCopiedGap = new StringBuilder ();\r
-var intronLength = 0;\r
-while (basesWritten < mappedCodonEnd && thisSeqPos < thisSeq.length) {\r
-var c = thisSeq[thisSeqPos++];\r
-if (c != myGapChar) {\r
-basesWritten++;\r
-if (basesWritten < mappedCodonStart) {\r
-if (preserveUnmappedGaps && trailingCopiedGap.length () > 0) {\r
-thisAligned.append (trailingCopiedGap.toString ());\r
-intronLength += trailingCopiedGap.length ();\r
-trailingCopiedGap = new StringBuilder ();\r
-}intronLength++;\r
-inExon = false;\r
-} else {\r
-var startOfCodon = basesWritten == mappedCodonStart;\r
-var gapsToAdd = jalview.analysis.AlignmentUtils.calculateGapsToInsert (preserveMappedGaps, preserveUnmappedGaps, sourceGapMappedLength, inExon, trailingCopiedGap.length (), intronLength, startOfCodon);\r
-for (var i = 0; i < gapsToAdd; i++) {\r
-thisAligned.append (myGapChar);\r
-}\r
-sourceGapMappedLength = 0;\r
-inExon = true;\r
-}thisAligned.append (c);\r
-trailingCopiedGap = new StringBuilder ();\r
-} else {\r
-if (inExon && preserveMappedGaps) {\r
-trailingCopiedGap.append (myGapChar);\r
-} else if (!inExon && preserveUnmappedGaps) {\r
-trailingCopiedGap.append (myGapChar);\r
-}}}\r
-}\r
-while (thisSeqPos < thisSeq.length) {\r
-var c = thisSeq[thisSeqPos++];\r
-if (c != myGapChar || preserveUnmappedGaps) {\r
-thisAligned.append (c);\r
-}}\r
-alignTo.setSequence ( String.instantialize (thisAligned));\r
-}, "jalview.datamodel.SequenceI,jalview.datamodel.SequenceI,jalview.datamodel.AlignedCodonFrame,~S,~S,~B,~B");\r
-c$.calculateGapsToInsert = Clazz.defineMethod (c$, "calculateGapsToInsert", \r
-function (preserveMappedGaps, preserveUnmappedGaps, sourceGapMappedLength, inExon, trailingGapLength, intronLength, startOfCodon) {\r
-var gapsToAdd = 0;\r
-if (startOfCodon) {\r
-if (inExon && !preserveMappedGaps) {\r
-trailingGapLength = 0;\r
-}if (!inExon && !(preserveMappedGaps && preserveUnmappedGaps)) {\r
-trailingGapLength = 0;\r
-}if (inExon) {\r
-gapsToAdd = Math.max (sourceGapMappedLength, trailingGapLength);\r
-} else {\r
-if (intronLength + trailingGapLength <= sourceGapMappedLength) {\r
-gapsToAdd = sourceGapMappedLength - intronLength;\r
-} else {\r
-gapsToAdd = Math.min (intronLength + trailingGapLength - sourceGapMappedLength, trailingGapLength);\r
-}}} else {\r
-if (!preserveMappedGaps) {\r
-trailingGapLength = 0;\r
-}gapsToAdd = Math.max (sourceGapMappedLength, trailingGapLength);\r
-}return gapsToAdd;\r
-}, "~B,~B,~N,~B,~N,~N,~B");\r
-c$.getAlignedTranslation = Clazz.defineMethod (c$, "getAlignedTranslation", \r
-function (sequences, gapCharacter, mappings) {\r
-var alignedSeqs = new java.util.ArrayList ();\r
-for (var seq, $seq = sequences.iterator (); $seq.hasNext () && ((seq = $seq.next ()) || true);) {\r
-var mapped = jalview.analysis.AlignmentUtils.getAlignedTranslation (seq, gapCharacter, mappings);\r
-alignedSeqs.addAll (mapped);\r
-}\r
-return alignedSeqs;\r
-}, "java.util.List,~S,java.util.Set");\r
-c$.getAlignedTranslation = Clazz.defineMethod (c$, "getAlignedTranslation", \r
-function (seq, gapCharacter, mappings) {\r
-var result = new java.util.ArrayList ();\r
-for (var mapping, $mapping = mappings.iterator (); $mapping.hasNext () && ((mapping = $mapping.next ()) || true);) {\r
-if (mapping.involvesSequence (seq)) {\r
-var mapped = jalview.analysis.AlignmentUtils.getAlignedTranslation (seq, gapCharacter, mapping);\r
-if (mapped != null) {\r
-result.add (mapped);\r
-}}}\r
-return result;\r
-}, "jalview.datamodel.SequenceI,~S,java.util.Set");\r
-c$.getAlignedTranslation = Clazz.defineMethod (c$, "getAlignedTranslation", \r
-function (seq, gapCharacter, mapping) {\r
-var gap = String.valueOf (gapCharacter);\r
-var toDna = false;\r
-var fromRatio = 1;\r
-var mapTo = mapping.getDnaForAaSeq (seq);\r
-if (mapTo != null) {\r
-toDna = true;\r
-gap = String.valueOf ( Clazz.newCharArray (-1, [gapCharacter, gapCharacter, gapCharacter]));\r
-} else {\r
-mapTo = mapping.getAaForDnaSeq (seq);\r
-fromRatio = 3;\r
-}var newseq = new StringBuilder (seq.getLength () * (toDna ? 3 : 1));\r
-var residueNo = 0;\r
-var phrase = Clazz.newIntArray (fromRatio, 0);\r
-var phraseOffset = 0;\r
-var gapWidth = 0;\r
-var first = true;\r
-var alignedSeq = new jalview.datamodel.Sequence ("", "");\r
-for (var c, $c = 0, $$c = seq.getSequence (); $c < $$c.length && ((c = $$c[$c]) || true); $c++) {\r
-if (c == gapCharacter) {\r
-gapWidth++;\r
-if (gapWidth >= fromRatio) {\r
-newseq.append (gap);\r
-gapWidth = 0;\r
-}} else {\r
-phrase[phraseOffset++] = residueNo + 1;\r
-if (phraseOffset == fromRatio) {\r
-var sr = new jalview.datamodel.SearchResults ();\r
-for (var pos, $pos = 0, $$pos = phrase; $pos < $$pos.length && ((pos = $$pos[$pos]) || true); $pos++) {\r
-mapping.markMappedRegion (seq, pos, sr);\r
-}\r
-newseq.append (sr.toString ());\r
-if (first) {\r
-first = false;\r
-var mappedTo = sr.getResultSequence (0);\r
-alignedSeq.setName (mappedTo.getName ());\r
-alignedSeq.setDescription (mappedTo.getDescription ());\r
-alignedSeq.setDatasetSequence (mappedTo);\r
-}phraseOffset = 0;\r
-}residueNo++;\r
-}}\r
-alignedSeq.setSequence (newseq.toString ());\r
-return alignedSeq;\r
-}, "jalview.datamodel.SequenceI,~S,jalview.datamodel.AlignedCodonFrame");\r
-c$.alignProteinAsDna = Clazz.defineMethod (c$, "alignProteinAsDna", \r
-function (protein, dna) {\r
-var mappings = protein.getCodonFrames ();\r
-var alignedCodons = new java.util.TreeMap ( new jalview.analysis.CodonComparator ());\r
-for (var dnaSeq, $dnaSeq = dna.getSequences ().iterator (); $dnaSeq.hasNext () && ((dnaSeq = $dnaSeq.next ()) || true);) {\r
-for (var mapping, $mapping = mappings.iterator (); $mapping.hasNext () && ((mapping = $mapping.next ()) || true);) {\r
-var seqMap = mapping.getMappingForSequence (dnaSeq);\r
-var prot = mapping.findAlignedSequence (dnaSeq.getDatasetSequence (), protein);\r
-if (prot != null) {\r
-jalview.analysis.AlignmentUtils.addCodonPositions (dnaSeq, prot, protein.getGapCharacter (), seqMap, alignedCodons);\r
-}}\r
-}\r
-return jalview.analysis.AlignmentUtils.alignProteinAs (protein, alignedCodons);\r
-}, "jalview.datamodel.AlignmentI,jalview.datamodel.AlignmentI");\r
-c$.alignProteinAs = Clazz.defineMethod (c$, "alignProteinAs", \r
-function (protein, alignedCodons) {\r
-var alignedWidth = alignedCodons.size ();\r
-var gaps = Clazz.newCharArray (alignedWidth, '\0');\r
-java.util.Arrays.fill (gaps, protein.getGapCharacter ());\r
-var allGaps = String.valueOf (gaps);\r
-for (var seq, $seq = protein.getSequences ().iterator (); $seq.hasNext () && ((seq = $seq.next ()) || true);) {\r
-seq.setSequence (allGaps);\r
-}\r
-var column = 0;\r
-for (var codon, $codon = alignedCodons.keySet ().iterator (); $codon.hasNext () && ((codon = $codon.next ()) || true);) {\r
-var columnResidues = alignedCodons.get (codon);\r
-for (var entry, $entry = columnResidues.entrySet ().iterator (); $entry.hasNext () && ((entry = $entry.next ()) || true);) {\r
-entry.getKey ().getSequence ()[column] = entry.getValue ().charAt (0);\r
-}\r
-column++;\r
-}\r
-return 0;\r
-}, "jalview.datamodel.AlignmentI,java.util.Map");\r
-c$.addCodonPositions = Clazz.defineMethod (c$, "addCodonPositions", \r
-function (dna, protein, gapChar, seqMap, alignedCodons) {\r
-var codons = seqMap.getCodonIterator (dna, gapChar);\r
-while (codons.hasNext ()) {\r
-var codon = codons.next ();\r
-var seqProduct = alignedCodons.get (codon);\r
-if (seqProduct == null) {\r
-seqProduct = new java.util.HashMap ();\r
-alignedCodons.put (codon, seqProduct);\r
-}seqProduct.put (protein, codon.product);\r
-}\r
-}, "jalview.datamodel.SequenceI,jalview.datamodel.SequenceI,~S,jalview.datamodel.Mapping,java.util.Map");\r
-c$.isMappable = Clazz.defineMethod (c$, "isMappable", \r
-function (al1, al2) {\r
-if (al1.isNucleotide () == al2.isNucleotide ()) {\r
-return false;\r
-}var dna = al1.isNucleotide () ? al1 : al2;\r
-var protein = dna === al1 ? al2 : al1;\r
-var mappings = protein.getCodonFrames ();\r
-for (var dnaSeq, $dnaSeq = dna.getSequences ().iterator (); $dnaSeq.hasNext () && ((dnaSeq = $dnaSeq.next ()) || true);) {\r
-for (var proteinSeq, $proteinSeq = protein.getSequences ().iterator (); $proteinSeq.hasNext () && ((proteinSeq = $proteinSeq.next ()) || true);) {\r
-if (jalview.analysis.AlignmentUtils.isMappable (dnaSeq, proteinSeq, mappings)) {\r
-return true;\r
-}}\r
-}\r
-return false;\r
-}, "jalview.datamodel.AlignmentI,jalview.datamodel.AlignmentI");\r
-c$.isMappable = Clazz.defineMethod (c$, "isMappable", \r
-function (dnaSeq, proteinSeq, mappings) {\r
-var dnaDs = dnaSeq.getDatasetSequence () == null ? dnaSeq : dnaSeq.getDatasetSequence ();\r
-var proteinDs = proteinSeq.getDatasetSequence () == null ? proteinSeq : proteinSeq.getDatasetSequence ();\r
-for (var mapping, $mapping = mappings.iterator (); $mapping.hasNext () && ((mapping = $mapping.next ()) || true);) {\r
-if (proteinDs === mapping.getAaForDnaSeq (dnaDs)) {\r
-return true;\r
-}}\r
-return jalview.analysis.AlignmentUtils.mapProteinToCdna (proteinDs, dnaDs) != null;\r
-}, "jalview.datamodel.SequenceI,jalview.datamodel.SequenceI,java.util.Set");\r
-c$.findAddableReferenceAnnotations = Clazz.defineMethod (c$, "findAddableReferenceAnnotations", \r
-function (sequenceScope, labelForCalcId, candidates, al) {\r
-if (sequenceScope == null) {\r
-return;\r
-}for (var seq, $seq = sequenceScope.iterator (); $seq.hasNext () && ((seq = $seq.next ()) || true);) {\r
-var dataset = seq.getDatasetSequence ();\r
-if (dataset == null) {\r
-continue;\r
-}var datasetAnnotations = dataset.getAnnotation ();\r
-if (datasetAnnotations == null) {\r
-continue;\r
-}var result = new java.util.ArrayList ();\r
-for (var dsann, $dsann = 0, $$dsann = datasetAnnotations; $dsann < $$dsann.length && ((dsann = $$dsann[$dsann]) || true); $dsann++) {\r
-var matchedAlignmentAnnotations = al.findAnnotations (seq, dsann.getCalcId (), dsann.label);\r
-if (!matchedAlignmentAnnotations.iterator ().hasNext ()) {\r
-result.add (dsann);\r
-if (labelForCalcId != null) {\r
-labelForCalcId.put (dsann.getCalcId (), dsann.label);\r
-}}}\r
-if (!result.isEmpty ()) {\r
-candidates.put (seq, result);\r
-}}\r
-}, "java.util.List,java.util.Map,java.util.Map,jalview.datamodel.AlignmentI");\r
-c$.addReferenceAnnotations = Clazz.defineMethod (c$, "addReferenceAnnotations", \r
-function (annotations, alignment, selectionGroup) {\r
-for (var seq, $seq = annotations.keySet ().iterator (); $seq.hasNext () && ((seq = $seq.next ()) || true);) {\r
-for (var ann, $ann = annotations.get (seq).iterator (); $ann.hasNext () && ((ann = $ann.next ()) || true);) {\r
-var copyAnn = new jalview.datamodel.AlignmentAnnotation (ann);\r
-var startRes = 0;\r
-var endRes = ann.annotations.length;\r
-if (selectionGroup != null) {\r
-startRes = selectionGroup.getStartRes ();\r
-endRes = selectionGroup.getEndRes ();\r
-}copyAnn.restrict (startRes, endRes);\r
-if (!seq.hasAnnotation (ann)) {\r
-seq.addAlignmentAnnotation (copyAnn);\r
-}copyAnn.adjustForAlignment ();\r
-alignment.addAnnotation (copyAnn);\r
-copyAnn.visible = true;\r
-}\r
-}\r
-}, "java.util.Map,jalview.datamodel.AlignmentI,jalview.datamodel.SequenceGroup");\r
-c$.showOrHideSequenceAnnotations = Clazz.defineMethod (c$, "showOrHideSequenceAnnotations", \r
-function (al, types, forSequences, anyType, doShow) {\r
-for (var aa, $aa = 0, $$aa = al.getAlignmentAnnotation (); $aa < $$aa.length && ((aa = $$aa[$aa]) || true); $aa++) {\r
-if (anyType || types.contains (aa.label)) {\r
-if ((aa.sequenceRef != null) && (forSequences == null || forSequences.contains (aa.sequenceRef))) {\r
-aa.visible = doShow;\r
-}}}\r
-}, "jalview.datamodel.AlignmentI,java.util.Collection,java.util.List,~B,~B");\r
-c$.haveCrossRef = Clazz.defineMethod (c$, "haveCrossRef", \r
-function (seq1, seq2) {\r
-return jalview.analysis.AlignmentUtils.hasCrossRef (seq1, seq2) || jalview.analysis.AlignmentUtils.hasCrossRef (seq2, seq1);\r
-}, "jalview.datamodel.SequenceI,jalview.datamodel.SequenceI");\r
-c$.hasCrossRef = Clazz.defineMethod (c$, "hasCrossRef", \r
-function (seq1, seq2) {\r
-if (seq1 == null || seq2 == null) {\r
-return false;\r
-}var name = seq2.getName ();\r
-var xrefs = seq1.getDBRef ();\r
-if (xrefs != null) {\r
-for (var xref, $xref = 0, $$xref = xrefs; $xref < $$xref.length && ((xref = $$xref[$xref]) || true); $xref++) {\r
-var xrefName = xref.getSource () + "|" + xref.getAccessionId ();\r
-if (xrefName.equalsIgnoreCase (name)) {\r
-return true;\r
-}}\r
-}return false;\r
-}, "jalview.datamodel.SequenceI,jalview.datamodel.SequenceI");\r
-c$.makeExonAlignment = Clazz.defineMethod (c$, "makeExonAlignment", \r
-function (dna, mappings) {\r
-var newMappings = new java.util.LinkedHashSet ();\r
-var exonSequences = new java.util.ArrayList ();\r
-for (var dnaSeq, $dnaSeq = 0, $$dnaSeq = dna; $dnaSeq < $$dnaSeq.length && ((dnaSeq = $$dnaSeq[$dnaSeq]) || true); $dnaSeq++) {\r
-var ds = dnaSeq.getDatasetSequence ();\r
-var seqMappings = jalview.util.MappingUtils.findMappingsForSequence (ds, mappings);\r
-for (var acf, $acf = seqMappings.iterator (); $acf.hasNext () && ((acf = $acf.next ()) || true);) {\r
-var newMapping = new jalview.datamodel.AlignedCodonFrame ();\r
-var mappedExons = jalview.analysis.AlignmentUtils.makeExonSequences (ds, acf, newMapping);\r
-if (!mappedExons.isEmpty ()) {\r
-exonSequences.addAll (mappedExons);\r
-newMappings.add (newMapping);\r
-}}\r
-}\r
-var al = new jalview.datamodel.Alignment (exonSequences.toArray ( new Array (exonSequences.size ())));\r
-al.setDataset (null);\r
-mappings.clear ();\r
-mappings.addAll (newMappings);\r
-return al;\r
-}, "~A,java.util.Set");\r
-c$.makeExonSequences = Clazz.defineMethod (c$, "makeExonSequences", \r
-function (dnaSeq, mapping, newMapping) {\r
-var exonSequences = new java.util.ArrayList ();\r
-var seqMappings = mapping.getMappingsForSequence (dnaSeq);\r
-var dna = dnaSeq.getSequence ();\r
-for (var seqMapping, $seqMapping = seqMappings.iterator (); $seqMapping.hasNext () && ((seqMapping = $seqMapping.next ()) || true);) {\r
-var newSequence = new StringBuilder (dnaSeq.getLength ());\r
-var dnaExonRanges = seqMapping.getMap ().getFromRanges ();\r
-for (var range, $range = dnaExonRanges.iterator (); $range.hasNext () && ((range = $range.next ()) || true);) {\r
-for (var pos = range[0]; pos <= range[1]; pos++) {\r
-newSequence.append (dna[pos - 1]);\r
-}\r
-}\r
-var exon = new jalview.datamodel.Sequence (dnaSeq.getName (), newSequence.toString ());\r
-var cdsAccId = jalview.datamodel.FeatureProperties.getCodingFeature (jalview.datamodel.DBRefSource.EMBL);\r
-var cdsRefs = jalview.util.DBRefUtils.selectRefs (seqMapping.getTo ().getDBRef (), jalview.datamodel.DBRefSource.CODINGDBS);\r
-if (cdsRefs != null) {\r
-for (var cdsRef, $cdsRef = 0, $$cdsRef = cdsRefs; $cdsRef < $$cdsRef.length && ((cdsRef = $$cdsRef[$cdsRef]) || true); $cdsRef++) {\r
-exon.addDBRef ( new jalview.datamodel.DBRefEntry (cdsRef));\r
-cdsAccId = cdsRef.getAccessionId ();\r
-}\r
-}exon.setName (exon.getName () + "|" + cdsAccId);\r
-exon.createDatasetSequence ();\r
-var exonRange = new java.util.ArrayList ();\r
-exonRange.add ( Clazz.newIntArray (-1, [1, newSequence.length ()]));\r
-var map = new jalview.util.MapList (exonRange, seqMapping.getMap ().getToRanges (), 3, 1);\r
-newMapping.addMap (exon.getDatasetSequence (), seqMapping.getTo (), map);\r
-var cdsToDnaMap = new jalview.util.MapList (dnaExonRanges, exonRange, 1, 1);\r
-newMapping.addMap (dnaSeq, exon.getDatasetSequence (), cdsToDnaMap);\r
-exonSequences.add (exon);\r
-}\r
-return exonSequences;\r
-}, "jalview.datamodel.SequenceI,jalview.datamodel.AlignedCodonFrame,jalview.datamodel.AlignedCodonFrame");\r
-});\r
+Clazz.declarePackage ("jalview.analysis");
+Clazz.load (null, "jalview.analysis.AlignmentUtils", ["jalview.analysis.CodonComparator", "jalview.datamodel.AlignedCodonFrame", "$.Alignment", "$.AlignmentAnnotation", "$.DBRefEntry", "$.DBRefSource", "$.FeatureProperties", "$.SearchResults", "$.Sequence", "jalview.schemes.ResidueProperties", "jalview.util.DBRefUtils", "$.MapList", "$.MappingUtils", "java.lang.StringBuilder", "java.util.ArrayList", "$.Arrays", "$.HashMap", "$.HashSet", "$.LinkedHashMap", "$.LinkedHashSet", "$.TreeMap"], function () {
+c$ = Clazz.declareType (jalview.analysis, "AlignmentUtils");
+c$.expandContext = Clazz.defineMethod (c$, "expandContext",
+function (core, flankSize) {
+var sq = new java.util.ArrayList ();
+var maxoffset = 0;
+for (var s, $s = core.getSequences ().iterator (); $s.hasNext () && ((s = $s.next ()) || true);) {
+var newSeq = s.deriveSequence ();
+var newSeqStart = newSeq.getStart () - 1;
+if (newSeqStart > maxoffset && newSeq.getDatasetSequence ().getStart () < s.getStart ()) {
+maxoffset = newSeqStart;
+}sq.add (newSeq);
+}
+if (flankSize > -1) {
+maxoffset = Math.min (maxoffset, flankSize);
+}for (var s, $s = sq.iterator (); $s.hasNext () && ((s = $s.next ()) || true);) {
+var ds = s;
+while (ds.getDatasetSequence () != null) {
+ds = ds.getDatasetSequence ();
+}
+var s_end = s.findPosition (s.getStart () + s.getLength ());
+var ustream_ds = s.getStart () - ds.getStart ();
+var dstream_ds = ds.getEnd () - s_end;
+var offset = maxoffset - ustream_ds;
+if (flankSize >= 0) {
+if (flankSize < ustream_ds) {
+offset = maxoffset - flankSize;
+ustream_ds = flankSize;
+}if (flankSize <= dstream_ds) {
+dstream_ds = flankSize - 1;
+}}var upstream = String.instantialize (ds.getSequence (s.getStart () - 1 - ustream_ds, s.getStart () - 1)).toLowerCase ().toCharArray ();
+var downstream = String.instantialize (ds.getSequence (s_end - 1, s_end + dstream_ds)).toLowerCase ().toCharArray ();
+var coreseq = s.getSequence ();
+var nseq = Clazz.newCharArray (offset + upstream.length + downstream.length + coreseq.length, '\0');
+var c = core.getGapCharacter ();
+var p = 0;
+for (; p < offset; p++) {
+nseq[p] = c;
+}
+System.arraycopy (upstream, 0, nseq, p, upstream.length);
+System.arraycopy (coreseq, 0, nseq, p + upstream.length, coreseq.length);
+System.arraycopy (downstream, 0, nseq, p + coreseq.length + upstream.length, downstream.length);
+s.setSequence ( String.instantialize (nseq));
+s.setStart (s.getStart () - ustream_ds);
+s.setEnd (s_end + downstream.length);
+}
+var newAl = new jalview.datamodel.Alignment (sq.toArray ( new Array (0)));
+for (var s, $s = sq.iterator (); $s.hasNext () && ((s = $s.next ()) || true);) {
+if (s.getAnnotation () != null) {
+for (var aa, $aa = 0, $$aa = s.getAnnotation (); $aa < $$aa.length && ((aa = $$aa[$aa]) || true); $aa++) {
+aa.adjustForAlignment ();
+newAl.addAnnotation (aa);
+}
+}}
+newAl.setDataset (core.getDataset ());
+return newAl;
+}, "jalview.datamodel.AlignmentI,~N");
+c$.getSequenceIndex = Clazz.defineMethod (c$, "getSequenceIndex",
+function (al, seq) {
+var result = -1;
+var pos = 0;
+for (var alSeq, $alSeq = al.getSequences ().iterator (); $alSeq.hasNext () && ((alSeq = $alSeq.next ()) || true);) {
+if (alSeq === seq) {
+result = pos;
+break;
+}pos++;
+}
+return result;
+}, "jalview.datamodel.AlignmentI,jalview.datamodel.SequenceI");
+c$.getSequencesByName = Clazz.defineMethod (c$, "getSequencesByName",
+function (al) {
+var theMap = new java.util.LinkedHashMap ();
+for (var seq, $seq = al.getSequences ().iterator (); $seq.hasNext () && ((seq = $seq.next ()) || true);) {
+var name = seq.getName ();
+if (name != null) {
+var seqs = theMap.get (name);
+if (seqs == null) {
+seqs = new java.util.ArrayList ();
+theMap.put (name, seqs);
+}seqs.add (seq);
+}}
+return theMap;
+}, "jalview.datamodel.AlignmentI");
+c$.mapProteinToCdna = Clazz.defineMethod (c$, "mapProteinToCdna",
+function (proteinAlignment, cdnaAlignment) {
+if (proteinAlignment == null || cdnaAlignment == null) {
+return false;
+}var mappedDna = new java.util.HashSet ();
+var mappedProtein = new java.util.HashSet ();
+var mappingPerformed = jalview.analysis.AlignmentUtils.mapProteinToCdna (proteinAlignment, cdnaAlignment, mappedDna, mappedProtein, true);
+mappingPerformed = new Boolean (mappingPerformed | jalview.analysis.AlignmentUtils.mapProteinToCdna (proteinAlignment, cdnaAlignment, mappedDna, mappedProtein, false)).valueOf ();
+return mappingPerformed;
+}, "jalview.datamodel.AlignmentI,jalview.datamodel.AlignmentI");
+c$.mapProteinToCdna = Clazz.defineMethod (c$, "mapProteinToCdna",
+function (proteinAlignment, cdnaAlignment, mappedDna, mappedProtein, xrefsOnly) {
+var mappingPerformed = false;
+var thisSeqs = proteinAlignment.getSequences ();
+for (var aaSeq, $aaSeq = thisSeqs.iterator (); $aaSeq.hasNext () && ((aaSeq = $aaSeq.next ()) || true);) {
+var proteinMapped = false;
+var acf = new jalview.datamodel.AlignedCodonFrame ();
+for (var cdnaSeq, $cdnaSeq = cdnaAlignment.getSequences ().iterator (); $cdnaSeq.hasNext () && ((cdnaSeq = $cdnaSeq.next ()) || true);) {
+if (xrefsOnly && !jalview.analysis.AlignmentUtils.haveCrossRef (aaSeq, cdnaSeq)) {
+continue;
+}if (!xrefsOnly && (mappedProtein.contains (aaSeq) || mappedDna.contains (cdnaSeq))) {
+continue;
+}if (!jalview.analysis.AlignmentUtils.mappingExists (proteinAlignment.getCodonFrames (), aaSeq.getDatasetSequence (), cdnaSeq.getDatasetSequence ())) {
+var map = jalview.analysis.AlignmentUtils.mapProteinToCdna (aaSeq, cdnaSeq);
+if (map != null) {
+acf.addMap (cdnaSeq, aaSeq, map);
+mappingPerformed = true;
+proteinMapped = true;
+mappedDna.add (cdnaSeq);
+mappedProtein.add (aaSeq);
+}}}
+if (proteinMapped) {
+proteinAlignment.addCodonFrame (acf);
+}}
+return mappingPerformed;
+}, "jalview.datamodel.AlignmentI,jalview.datamodel.AlignmentI,java.util.Set,java.util.Set,~B");
+c$.mappingExists = Clazz.defineMethod (c$, "mappingExists",
+function (set, aaSeq, cdnaSeq) {
+if (set != null) {
+for (var acf, $acf = set.iterator (); $acf.hasNext () && ((acf = $acf.next ()) || true);) {
+if (cdnaSeq === acf.getDnaForAaSeq (aaSeq)) {
+return true;
+}}
+}return false;
+}, "java.util.Set,jalview.datamodel.SequenceI,jalview.datamodel.SequenceI");
+c$.mapProteinToCdna = Clazz.defineMethod (c$, "mapProteinToCdna",
+function (proteinSeq, cdnaSeq) {
+var proteinDataset = proteinSeq.getDatasetSequence ();
+var aaSeqChars = proteinDataset != null ? proteinDataset.getSequence () : proteinSeq.getSequence ();
+var cdnaDataset = cdnaSeq.getDatasetSequence ();
+var cdnaSeqChars = cdnaDataset != null ? cdnaDataset.getSequence () : cdnaSeq.getSequence ();
+if (aaSeqChars == null || cdnaSeqChars == null) {
+return null;
+}var mappedLength = 3 * aaSeqChars.length;
+var cdnaLength = cdnaSeqChars.length;
+var cdnaStart = 1;
+var cdnaEnd = cdnaLength;
+var proteinStart = 1;
+var proteinEnd = aaSeqChars.length;
+if (cdnaLength != mappedLength && cdnaLength > 2) {
+var lastCodon = String.valueOf (cdnaSeqChars, cdnaLength - 3, 3).toUpperCase ();
+for (var stop, $stop = jalview.schemes.ResidueProperties.STOP.iterator (); $stop.hasNext () && ((stop = $stop.next ()) || true);) {
+if (lastCodon.equals (stop)) {
+cdnaEnd -= 3;
+cdnaLength -= 3;
+break;
+}}
+}if (cdnaLength != mappedLength && cdnaLength > 2 && String.valueOf (cdnaSeqChars, 0, 3).toUpperCase ().equals (jalview.schemes.ResidueProperties.START)) {
+cdnaStart += 3;
+cdnaLength -= 3;
+}if (cdnaLength != mappedLength) {
+return null;
+}if (!jalview.analysis.AlignmentUtils.translatesAs (cdnaSeqChars, cdnaStart - 1, aaSeqChars)) {
+return null;
+}var map = new jalview.util.MapList ( Clazz.newIntArray (-1, [cdnaStart, cdnaEnd]), Clazz.newIntArray (-1, [1, proteinEnd]), 3, 1);
+return map;
+}, "jalview.datamodel.SequenceI,jalview.datamodel.SequenceI");
+c$.translatesAs = Clazz.defineMethod (c$, "translatesAs",
+function (cdnaSeqChars, cdnaStart, aaSeqChars) {
+var aaResidue = 0;
+for (var i = cdnaStart; i < cdnaSeqChars.length - 2 && aaResidue < aaSeqChars.length; i += 3, aaResidue++) {
+var codon = String.valueOf (cdnaSeqChars, i, 3);
+var translated = jalview.schemes.ResidueProperties.codonTranslate (codon);
+var aaRes = aaSeqChars[aaResidue];
+if ((translated == null || "STOP".equals (translated)) && aaRes == 'X') {
+continue;
+}if (translated == null || !(aaRes == translated.charAt (0))) {
+return false;
+}}
+return (aaResidue == aaSeqChars.length);
+}, "~A,~N,~A");
+c$.alignSequenceAs = Clazz.defineMethod (c$, "alignSequenceAs",
+function (seq, al, gap, preserveMappedGaps, preserveUnmappedGaps) {
+var mappings = al.getCodonFrame (seq);
+if (mappings == null || mappings.isEmpty ()) {
+return false;
+}var alignFrom = null;
+var mapping = null;
+for (var mp, $mp = mappings.iterator (); $mp.hasNext () && ((mp = $mp.next ()) || true);) {
+alignFrom = mp.findAlignedSequence (seq.getDatasetSequence (), al);
+if (alignFrom != null) {
+mapping = mp;
+break;
+}}
+if (alignFrom == null) {
+return false;
+}jalview.analysis.AlignmentUtils.alignSequenceAs (seq, alignFrom, mapping, gap, al.getGapCharacter (), preserveMappedGaps, preserveUnmappedGaps);
+return true;
+}, "jalview.datamodel.SequenceI,jalview.datamodel.AlignmentI,~S,~B,~B");
+c$.alignSequenceAs = Clazz.defineMethod (c$, "alignSequenceAs",
+function (alignTo, alignFrom, mapping, myGap, sourceGap, preserveMappedGaps, preserveUnmappedGaps) {
+var thisSeq = alignTo.getSequence ();
+var thatAligned = alignFrom.getSequence ();
+var thisAligned = new StringBuilder (2 * thisSeq.length);
+var thisSeqPos = 0;
+var sourceDsPos = 0;
+var basesWritten = 0;
+var myGapChar = myGap.charAt (0);
+var ratio = myGap.length;
+var sourceGapMappedLength = 0;
+var inExon = false;
+for (var sourceChar, $sourceChar = 0, $$sourceChar = thatAligned; $sourceChar < $$sourceChar.length && ((sourceChar = $$sourceChar[$sourceChar]) || true); $sourceChar++) {
+if (sourceChar == sourceGap) {
+sourceGapMappedLength += ratio;
+continue;
+}sourceDsPos++;
+var mappedPos = mapping.getMappedRegion (alignTo, alignFrom, sourceDsPos);
+if (mappedPos == null) {
+System.err.println ("Can't align: no codon mapping to residue " + sourceDsPos + "(" + sourceChar + ")");
+return;
+}var mappedCodonStart = mappedPos[0];
+var mappedCodonEnd = mappedPos[mappedPos.length - 1];
+var trailingCopiedGap = new StringBuilder ();
+var intronLength = 0;
+while (basesWritten < mappedCodonEnd && thisSeqPos < thisSeq.length) {
+var c = thisSeq[thisSeqPos++];
+if (c != myGapChar) {
+basesWritten++;
+if (basesWritten < mappedCodonStart) {
+if (preserveUnmappedGaps && trailingCopiedGap.length () > 0) {
+thisAligned.append (trailingCopiedGap.toString ());
+intronLength += trailingCopiedGap.length ();
+trailingCopiedGap = new StringBuilder ();
+}intronLength++;
+inExon = false;
+} else {
+var startOfCodon = basesWritten == mappedCodonStart;
+var gapsToAdd = jalview.analysis.AlignmentUtils.calculateGapsToInsert (preserveMappedGaps, preserveUnmappedGaps, sourceGapMappedLength, inExon, trailingCopiedGap.length (), intronLength, startOfCodon);
+for (var i = 0; i < gapsToAdd; i++) {
+thisAligned.append (myGapChar);
+}
+sourceGapMappedLength = 0;
+inExon = true;
+}thisAligned.append (c);
+trailingCopiedGap = new StringBuilder ();
+} else {
+if (inExon && preserveMappedGaps) {
+trailingCopiedGap.append (myGapChar);
+} else if (!inExon && preserveUnmappedGaps) {
+trailingCopiedGap.append (myGapChar);
+}}}
+}
+while (thisSeqPos < thisSeq.length) {
+var c = thisSeq[thisSeqPos++];
+if (c != myGapChar || preserveUnmappedGaps) {
+thisAligned.append (c);
+}}
+alignTo.setSequence ( String.instantialize (thisAligned));
+}, "jalview.datamodel.SequenceI,jalview.datamodel.SequenceI,jalview.datamodel.AlignedCodonFrame,~S,~S,~B,~B");
+c$.calculateGapsToInsert = Clazz.defineMethod (c$, "calculateGapsToInsert",
+function (preserveMappedGaps, preserveUnmappedGaps, sourceGapMappedLength, inExon, trailingGapLength, intronLength, startOfCodon) {
+var gapsToAdd = 0;
+if (startOfCodon) {
+if (inExon && !preserveMappedGaps) {
+trailingGapLength = 0;
+}if (!inExon && !(preserveMappedGaps && preserveUnmappedGaps)) {
+trailingGapLength = 0;
+}if (inExon) {
+gapsToAdd = Math.max (sourceGapMappedLength, trailingGapLength);
+} else {
+if (intronLength + trailingGapLength <= sourceGapMappedLength) {
+gapsToAdd = sourceGapMappedLength - intronLength;
+} else {
+gapsToAdd = Math.min (intronLength + trailingGapLength - sourceGapMappedLength, trailingGapLength);
+}}} else {
+if (!preserveMappedGaps) {
+trailingGapLength = 0;
+}gapsToAdd = Math.max (sourceGapMappedLength, trailingGapLength);
+}return gapsToAdd;
+}, "~B,~B,~N,~B,~N,~N,~B");
+c$.getAlignedTranslation = Clazz.defineMethod (c$, "getAlignedTranslation",
+function (sequences, gapCharacter, mappings) {
+var alignedSeqs = new java.util.ArrayList ();
+for (var seq, $seq = sequences.iterator (); $seq.hasNext () && ((seq = $seq.next ()) || true);) {
+var mapped = jalview.analysis.AlignmentUtils.getAlignedTranslation (seq, gapCharacter, mappings);
+alignedSeqs.addAll (mapped);
+}
+return alignedSeqs;
+}, "java.util.List,~S,java.util.Set");
+c$.getAlignedTranslation = Clazz.defineMethod (c$, "getAlignedTranslation",
+function (seq, gapCharacter, mappings) {
+var result = new java.util.ArrayList ();
+for (var mapping, $mapping = mappings.iterator (); $mapping.hasNext () && ((mapping = $mapping.next ()) || true);) {
+if (mapping.involvesSequence (seq)) {
+var mapped = jalview.analysis.AlignmentUtils.getAlignedTranslation (seq, gapCharacter, mapping);
+if (mapped != null) {
+result.add (mapped);
+}}}
+return result;
+}, "jalview.datamodel.SequenceI,~S,java.util.Set");
+c$.getAlignedTranslation = Clazz.defineMethod (c$, "getAlignedTranslation",
+function (seq, gapCharacter, mapping) {
+var gap = String.valueOf (gapCharacter);
+var toDna = false;
+var fromRatio = 1;
+var mapTo = mapping.getDnaForAaSeq (seq);
+if (mapTo != null) {
+toDna = true;
+gap = String.valueOf ( Clazz.newCharArray (-1, [gapCharacter, gapCharacter, gapCharacter]));
+} else {
+mapTo = mapping.getAaForDnaSeq (seq);
+fromRatio = 3;
+}var newseq = new StringBuilder (seq.getLength () * (toDna ? 3 : 1));
+var residueNo = 0;
+var phrase = Clazz.newIntArray (fromRatio, 0);
+var phraseOffset = 0;
+var gapWidth = 0;
+var first = true;
+var alignedSeq = new jalview.datamodel.Sequence ("", "");
+for (var c, $c = 0, $$c = seq.getSequence (); $c < $$c.length && ((c = $$c[$c]) || true); $c++) {
+if (c == gapCharacter) {
+gapWidth++;
+if (gapWidth >= fromRatio) {
+newseq.append (gap);
+gapWidth = 0;
+}} else {
+phrase[phraseOffset++] = residueNo + 1;
+if (phraseOffset == fromRatio) {
+var sr = new jalview.datamodel.SearchResults ();
+for (var pos, $pos = 0, $$pos = phrase; $pos < $$pos.length && ((pos = $$pos[$pos]) || true); $pos++) {
+mapping.markMappedRegion (seq, pos, sr);
+}
+newseq.append (sr.toString ());
+if (first) {
+first = false;
+var mappedTo = sr.getResultSequence (0);
+alignedSeq.setName (mappedTo.getName ());
+alignedSeq.setDescription (mappedTo.getDescription ());
+alignedSeq.setDatasetSequence (mappedTo);
+}phraseOffset = 0;
+}residueNo++;
+}}
+alignedSeq.setSequence (newseq.toString ());
+return alignedSeq;
+}, "jalview.datamodel.SequenceI,~S,jalview.datamodel.AlignedCodonFrame");
+c$.alignProteinAsDna = Clazz.defineMethod (c$, "alignProteinAsDna",
+function (protein, dna) {
+var mappings = protein.getCodonFrames ();
+var alignedCodons = new java.util.TreeMap ( new jalview.analysis.CodonComparator ());
+for (var dnaSeq, $dnaSeq = dna.getSequences ().iterator (); $dnaSeq.hasNext () && ((dnaSeq = $dnaSeq.next ()) || true);) {
+for (var mapping, $mapping = mappings.iterator (); $mapping.hasNext () && ((mapping = $mapping.next ()) || true);) {
+var seqMap = mapping.getMappingForSequence (dnaSeq);
+var prot = mapping.findAlignedSequence (dnaSeq.getDatasetSequence (), protein);
+if (prot != null) {
+jalview.analysis.AlignmentUtils.addCodonPositions (dnaSeq, prot, protein.getGapCharacter (), seqMap, alignedCodons);
+}}
+}
+return jalview.analysis.AlignmentUtils.alignProteinAs (protein, alignedCodons);
+}, "jalview.datamodel.AlignmentI,jalview.datamodel.AlignmentI");
+c$.alignProteinAs = Clazz.defineMethod (c$, "alignProteinAs",
+function (protein, alignedCodons) {
+var alignedWidth = alignedCodons.size ();
+var gaps = Clazz.newCharArray (alignedWidth, '\0');
+java.util.Arrays.fill (gaps, protein.getGapCharacter ());
+var allGaps = String.valueOf (gaps);
+for (var seq, $seq = protein.getSequences ().iterator (); $seq.hasNext () && ((seq = $seq.next ()) || true);) {
+seq.setSequence (allGaps);
+}
+var column = 0;
+for (var codon, $codon = alignedCodons.keySet ().iterator (); $codon.hasNext () && ((codon = $codon.next ()) || true);) {
+var columnResidues = alignedCodons.get (codon);
+for (var entry, $entry = columnResidues.entrySet ().iterator (); $entry.hasNext () && ((entry = $entry.next ()) || true);) {
+entry.getKey ().getSequence ()[column] = entry.getValue ().charAt (0);
+}
+column++;
+}
+return 0;
+}, "jalview.datamodel.AlignmentI,java.util.Map");
+c$.addCodonPositions = Clazz.defineMethod (c$, "addCodonPositions",
+function (dna, protein, gapChar, seqMap, alignedCodons) {
+var codons = seqMap.getCodonIterator (dna, gapChar);
+while (codons.hasNext ()) {
+var codon = codons.next ();
+var seqProduct = alignedCodons.get (codon);
+if (seqProduct == null) {
+seqProduct = new java.util.HashMap ();
+alignedCodons.put (codon, seqProduct);
+}seqProduct.put (protein, codon.product);
+}
+}, "jalview.datamodel.SequenceI,jalview.datamodel.SequenceI,~S,jalview.datamodel.Mapping,java.util.Map");
+c$.isMappable = Clazz.defineMethod (c$, "isMappable",
+function (al1, al2) {
+if (al1.isNucleotide () == al2.isNucleotide ()) {
+return false;
+}var dna = al1.isNucleotide () ? al1 : al2;
+var protein = dna === al1 ? al2 : al1;
+var mappings = protein.getCodonFrames ();
+for (var dnaSeq, $dnaSeq = dna.getSequences ().iterator (); $dnaSeq.hasNext () && ((dnaSeq = $dnaSeq.next ()) || true);) {
+for (var proteinSeq, $proteinSeq = protein.getSequences ().iterator (); $proteinSeq.hasNext () && ((proteinSeq = $proteinSeq.next ()) || true);) {
+if (jalview.analysis.AlignmentUtils.isMappable (dnaSeq, proteinSeq, mappings)) {
+return true;
+}}
+}
+return false;
+}, "jalview.datamodel.AlignmentI,jalview.datamodel.AlignmentI");
+c$.isMappable = Clazz.defineMethod (c$, "isMappable",
+function (dnaSeq, proteinSeq, mappings) {
+var dnaDs = dnaSeq.getDatasetSequence () == null ? dnaSeq : dnaSeq.getDatasetSequence ();
+var proteinDs = proteinSeq.getDatasetSequence () == null ? proteinSeq : proteinSeq.getDatasetSequence ();
+for (var mapping, $mapping = mappings.iterator (); $mapping.hasNext () && ((mapping = $mapping.next ()) || true);) {
+if (proteinDs === mapping.getAaForDnaSeq (dnaDs)) {
+return true;
+}}
+return jalview.analysis.AlignmentUtils.mapProteinToCdna (proteinDs, dnaDs) != null;
+}, "jalview.datamodel.SequenceI,jalview.datamodel.SequenceI,java.util.Set");
+c$.findAddableReferenceAnnotations = Clazz.defineMethod (c$, "findAddableReferenceAnnotations",
+function (sequenceScope, labelForCalcId, candidates, al) {
+if (sequenceScope == null) {
+return;
+}for (var seq, $seq = sequenceScope.iterator (); $seq.hasNext () && ((seq = $seq.next ()) || true);) {
+var dataset = seq.getDatasetSequence ();
+if (dataset == null) {
+continue;
+}var datasetAnnotations = dataset.getAnnotation ();
+if (datasetAnnotations == null) {
+continue;
+}var result = new java.util.ArrayList ();
+for (var dsann, $dsann = 0, $$dsann = datasetAnnotations; $dsann < $$dsann.length && ((dsann = $$dsann[$dsann]) || true); $dsann++) {
+var matchedAlignmentAnnotations = al.findAnnotations (seq, dsann.getCalcId (), dsann.label);
+if (!matchedAlignmentAnnotations.iterator ().hasNext ()) {
+result.add (dsann);
+if (labelForCalcId != null) {
+labelForCalcId.put (dsann.getCalcId (), dsann.label);
+}}}
+if (!result.isEmpty ()) {
+candidates.put (seq, result);
+}}
+}, "java.util.List,java.util.Map,java.util.Map,jalview.datamodel.AlignmentI");
+c$.addReferenceAnnotations = Clazz.defineMethod (c$, "addReferenceAnnotations",
+function (annotations, alignment, selectionGroup) {
+for (var seq, $seq = annotations.keySet ().iterator (); $seq.hasNext () && ((seq = $seq.next ()) || true);) {
+for (var ann, $ann = annotations.get (seq).iterator (); $ann.hasNext () && ((ann = $ann.next ()) || true);) {
+var copyAnn = new jalview.datamodel.AlignmentAnnotation (ann);
+var startRes = 0;
+var endRes = ann.annotations.length;
+if (selectionGroup != null) {
+startRes = selectionGroup.getStartRes ();
+endRes = selectionGroup.getEndRes ();
+}copyAnn.restrict (startRes, endRes);
+if (!seq.hasAnnotation (ann)) {
+seq.addAlignmentAnnotation (copyAnn);
+}copyAnn.adjustForAlignment ();
+alignment.addAnnotation (copyAnn);
+copyAnn.visible = true;
+}
+}
+}, "java.util.Map,jalview.datamodel.AlignmentI,jalview.datamodel.SequenceGroup");
+c$.showOrHideSequenceAnnotations = Clazz.defineMethod (c$, "showOrHideSequenceAnnotations",
+function (al, types, forSequences, anyType, doShow) {
+for (var aa, $aa = 0, $$aa = al.getAlignmentAnnotation (); $aa < $$aa.length && ((aa = $$aa[$aa]) || true); $aa++) {
+if (anyType || types.contains (aa.label)) {
+if ((aa.sequenceRef != null) && (forSequences == null || forSequences.contains (aa.sequenceRef))) {
+aa.visible = doShow;
+}}}
+}, "jalview.datamodel.AlignmentI,java.util.Collection,java.util.List,~B,~B");
+c$.haveCrossRef = Clazz.defineMethod (c$, "haveCrossRef",
+function (seq1, seq2) {
+return jalview.analysis.AlignmentUtils.hasCrossRef (seq1, seq2) || jalview.analysis.AlignmentUtils.hasCrossRef (seq2, seq1);
+}, "jalview.datamodel.SequenceI,jalview.datamodel.SequenceI");
+c$.hasCrossRef = Clazz.defineMethod (c$, "hasCrossRef",
+function (seq1, seq2) {
+if (seq1 == null || seq2 == null) {
+return false;
+}var name = seq2.getName ();
+var xrefs = seq1.getDBRef ();
+if (xrefs != null) {
+for (var xref, $xref = 0, $$xref = xrefs; $xref < $$xref.length && ((xref = $$xref[$xref]) || true); $xref++) {
+var xrefName = xref.getSource () + "|" + xref.getAccessionId ();
+if (xrefName.equalsIgnoreCase (name)) {
+return true;
+}}
+}return false;
+}, "jalview.datamodel.SequenceI,jalview.datamodel.SequenceI");
+c$.makeExonAlignment = Clazz.defineMethod (c$, "makeExonAlignment",
+function (dna, mappings) {
+var newMappings = new java.util.LinkedHashSet ();
+var exonSequences = new java.util.ArrayList ();
+for (var dnaSeq, $dnaSeq = 0, $$dnaSeq = dna; $dnaSeq < $$dnaSeq.length && ((dnaSeq = $$dnaSeq[$dnaSeq]) || true); $dnaSeq++) {
+var ds = dnaSeq.getDatasetSequence ();
+var seqMappings = jalview.util.MappingUtils.findMappingsForSequence (ds, mappings);
+for (var acf, $acf = seqMappings.iterator (); $acf.hasNext () && ((acf = $acf.next ()) || true);) {
+var newMapping = new jalview.datamodel.AlignedCodonFrame ();
+var mappedExons = jalview.analysis.AlignmentUtils.makeExonSequences (ds, acf, newMapping);
+if (!mappedExons.isEmpty ()) {
+exonSequences.addAll (mappedExons);
+newMappings.add (newMapping);
+}}
+}
+var al = new jalview.datamodel.Alignment (exonSequences.toArray ( new Array (exonSequences.size ())));
+al.setDataset (null);
+mappings.clear ();
+mappings.addAll (newMappings);
+return al;
+}, "~A,java.util.Set");
+c$.makeExonSequences = Clazz.defineMethod (c$, "makeExonSequences",
+function (dnaSeq, mapping, newMapping) {
+var exonSequences = new java.util.ArrayList ();
+var seqMappings = mapping.getMappingsForSequence (dnaSeq);
+var dna = dnaSeq.getSequence ();
+for (var seqMapping, $seqMapping = seqMappings.iterator (); $seqMapping.hasNext () && ((seqMapping = $seqMapping.next ()) || true);) {
+var newSequence = new StringBuilder (dnaSeq.getLength ());
+var dnaExonRanges = seqMapping.getMap ().getFromRanges ();
+for (var range, $range = dnaExonRanges.iterator (); $range.hasNext () && ((range = $range.next ()) || true);) {
+for (var pos = range[0]; pos <= range[1]; pos++) {
+newSequence.append (dna[pos - 1]);
+}
+}
+var exon = new jalview.datamodel.Sequence (dnaSeq.getName (), newSequence.toString ());
+var cdsAccId = jalview.datamodel.FeatureProperties.getCodingFeature (jalview.datamodel.DBRefSource.EMBL);
+var cdsRefs = jalview.util.DBRefUtils.selectRefs (seqMapping.getTo ().getDBRef (), jalview.datamodel.DBRefSource.CODINGDBS);
+if (cdsRefs != null) {
+for (var cdsRef, $cdsRef = 0, $$cdsRef = cdsRefs; $cdsRef < $$cdsRef.length && ((cdsRef = $$cdsRef[$cdsRef]) || true); $cdsRef++) {
+exon.addDBRef ( new jalview.datamodel.DBRefEntry (cdsRef));
+cdsAccId = cdsRef.getAccessionId ();
+}
+}exon.setName (exon.getName () + "|" + cdsAccId);
+exon.createDatasetSequence ();
+var exonRange = new java.util.ArrayList ();
+exonRange.add ( Clazz.newIntArray (-1, [1, newSequence.length ()]));
+var map = new jalview.util.MapList (exonRange, seqMapping.getMap ().getToRanges (), 3, 1);
+newMapping.addMap (exon.getDatasetSequence (), seqMapping.getTo (), map);
+var cdsToDnaMap = new jalview.util.MapList (dnaExonRanges, exonRange, 1, 1);
+newMapping.addMap (dnaSeq, exon.getDatasetSequence (), cdsToDnaMap);
+exonSequences.add (exon);
+}
+return exonSequences;
+}, "jalview.datamodel.SequenceI,jalview.datamodel.AlignedCodonFrame,jalview.datamodel.AlignedCodonFrame");
+});
git://source.jalview.org / jalviewjs.git / blobdiff