JAL-2743 fall back to parsing description if chromosome lookup fails;
authorgmungoc <g.m.carstairs@dundee.ac.uk>
Wed, 25 Oct 2017 14:37:53 +0000 (15:37 +0100)
committergmungoc <g.m.carstairs@dundee.ac.uk>
Wed, 25 Oct 2017 14:37:53 +0000 (15:37 +0100)
code formatting

src/jalview/ext/ensembl/EnsemblGene.java
src/jalview/io/vcf/VCFLoader.java

index 11322c8..cdcfa96 100644 (file)
@@ -166,25 +166,25 @@ public class EnsemblGene extends EnsemblSeqProxy
 
   /**
    * Calls the /lookup/id REST service, parses the response for gene
-   * coordinates, and if successful, adds these to the sequence
+   * coordinates, and if successful, adds these to the sequence. If this fails,
+   * fall back on trying to parse the sequence description in case it is in
+   * Ensembl-gene format e.g. chromosome:GRCh38:17:45051610:45109016:1.
    * 
    * @param seq
    * @param geneId
    */
   void findGeneLoci(SequenceI seq, String geneId)
   {
-    /*
-     * if sequence description is in gene loci format, parse this
-     * - but try remote lookup anyway as 'better' if available
-     */
-    parseChromosomeLocations(seq);
-
     GeneLociI geneLoci = new EnsemblLookup(getDomain()).getGeneLoci(geneId);
     if (geneLoci != null)
     {
       seq.setGeneLoci(geneLoci.getSpeciesId(), geneLoci.getAssemblyId(),
               geneLoci.getChromosomeId(), geneLoci.getMap());
     }
+    else
+    {
+      parseChromosomeLocations(seq);
+    }
   }
 
   /**
index af7bb0b..85bf7ef 100644 (file)
@@ -895,8 +895,7 @@ public class VCFLoader
      */
     EnsemblMap mapper = new EnsemblMap();
     int[] mapping = mapper.getAssemblyMapping(species, chromosome, fromRef,
-            toRef,
-            queryRange);
+            toRef, queryRange);
 
     if (mapping == null)
     {