X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=forester%2Fjava%2Fsrc%2Forg%2Fforester%2Futil%2FSequenceIdParser.java;h=d828a6a01bd8158b1d2350d2385823feded1470d;hb=aec065f948a075773794133f102ea19eb1d59f64;hp=0856f34d947a7648bc8d33c9da7d2e2beef151d0;hpb=e174d62534f7fc6f3de133d523a402a87735b27f;p=jalview.git

diff --git a/forester/java/src/org/forester/util/SequenceIdParser.java b/forester/java/src/org/forester/util/SequenceIdParser.java
index 0856f34..d828a6a 100644
--- a/forester/java/src/org/forester/util/SequenceIdParser.java
+++ b/forester/java/src/org/forester/util/SequenceIdParser.java
@@ -24,7 +24,7 @@
 // Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA
 //
 // Contact: phylosoft @ gmail . com
-// WWW: www.phylosoft.org/forester
+// WWW: https://sites.google.com/site/cmzmasek/home/software/forester
 
 package org.forester.util;
 
@@ -49,16 +49,21 @@ public final class SequenceIdParser {
     //Protein:    3 letters + 5 numerals
     //http://www.ncbi.nlm.nih.gov/Sequin/acc.html
     private final static Pattern GENBANK_NUCLEOTIDE_AC_PATTERN_1 = Pattern
-                                                                         .compile( "(?:\\A|.*[^a-zA-Z0-9])([A-Z]\\d{5})(?:[^a-zA-Z0-9]|\\Z)" );
+                                                                         .compile( "(?:\\A|.*[^a-zA-Z0-9])([A-Z]\\d{5}(?:\\.\\d+)?)(?:[^a-zA-Z0-9]|\\Z)" );
     private final static Pattern GENBANK_NUCLEOTIDE_AC_PATTERN_2 = Pattern
-                                                                         .compile( "(?:\\A|.*[^a-zA-Z0-9])([A-Z]{2}\\d{6})(?:[^a-zA-Z0-9]|\\Z)" );
+                                                                         .compile( "(?:\\A|.*[^a-zA-Z0-9])([A-Z]{2}\\d{6}(?:\\.\\d+)?)(?:[^a-zA-Z0-9]|\\Z)" );
     private final static Pattern GENBANK_PROTEIN_AC_PATTERN      = Pattern
-                                                                         .compile( "(?:\\A|.*[^a-zA-Z0-9])([A-Z]{3}\\d{5})(?:[^a-zA-Z0-9]|\\Z)" );
+                                                                         .compile( "(?:\\A|.*[^a-zA-Z0-9])([A-Z]{3}\\d{5}(?:\\.\\d+)?)(?:[^a-zA-Z0-9]|\\Z)" );
     // RefSeq accession numbers can be distinguished from GenBank accessions 
     // by their distinct prefix format of 2 characters followed by an
     // underscore character ('_'). For example, a RefSeq protein accession is NP_015325. 
     private final static Pattern REFSEQ_PATTERN                  = Pattern
                                                                          .compile( "(?:\\A|.*[^a-zA-Z0-9])([A-Z]{2}_\\d{6,})(?:[^a-zA-Z0-9]|\\Z)" );
+    // See: http://web.expasy.org/docs/userman.html#ID_line
+    private final static Pattern TREMBL_PATTERN                  = Pattern
+                                                                         .compile( "(?:\\A|.*[^a-zA-Z0-9])([A-Z][0-9][A-Z0-9]{3}[0-9])(?:[^a-zA-Z0-9]|\\Z)" );
+    private final static Pattern GI_PATTERN                      = Pattern
+                                                                         .compile( "(?:\\b|_)(?:GI|gi)[|_=:](\\d+)(?:\\b|_)" );
 
     /**
      * Returns null if no match.
@@ -73,9 +78,25 @@ public final class SequenceIdParser {
         if ( !ForesterUtil.isEmpty( v ) ) {
             return new Identifier( v, Identifier.REFSEQ );
         }
+        v = parseTrEMBLAccessor( s );
+        if ( !ForesterUtil.isEmpty( v ) ) {
+            return new Identifier( v, Identifier.SP );
+        }
         return null;
     }
 
+    public final static boolean isProtein( final String query ) {
+        final String r1 = parseRefSeqAccessor( query );
+        if ( !ForesterUtil.isEmpty( r1 ) && ( r1.charAt( 1 ) == 'P' ) ) {
+            return true;
+        }
+        final String r2 = parseTrEMBLAccessor( query );
+        if ( !ForesterUtil.isEmpty( r2 ) ) {
+            return true;
+        }
+        return GENBANK_PROTEIN_AC_PATTERN.matcher( query ).lookingAt();
+    }
+
     /**
      * Returns null if no match.
      * 
@@ -106,7 +127,7 @@ public final class SequenceIdParser {
      * Returns null if no match.
      * 
      */
-    private final static String parseRefSeqAccessor( final String query ) {
+    public final static String parseRefSeqAccessor( final String query ) {
         final Matcher m = REFSEQ_PATTERN.matcher( query );
         if ( m.lookingAt() ) {
             return m.group( 1 );
@@ -114,7 +135,27 @@ public final class SequenceIdParser {
         return null;
     }
 
+    /**
+     * Returns null if no match.
+     * 
+     */
+    private final static String parseTrEMBLAccessor( final String query ) {
+        final Matcher m = TREMBL_PATTERN.matcher( query );
+        if ( m.lookingAt() ) {
+            return m.group( 1 );
+        }
+        return null;
+    }
+
     private SequenceIdParser() {
         // Hiding the constructor.
     }
+
+    public static String parseGInumber( final String query ) {
+        final Matcher m = GI_PATTERN.matcher( query );
+        if ( m.find() ) {
+            return m.group( 1 );
+        }
+        return null;
+    }
 }