X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=forester%2Fjava%2Fsrc%2Forg%2Fforester%2Futil%2FSequenceIdParser.java;h=8def26026d0fbecc168a4d847cea97205a2ba237;hb=34e555e6895d06fa5b7de6a646c37b0817efafc6;hp=b8bebe1d90ae13009860e193670c28000b3f05a6;hpb=cb4ea01dc92bafecf44c8e4fa85d07528b80943f;p=jalview.git

diff --git a/forester/java/src/org/forester/util/SequenceIdParser.java b/forester/java/src/org/forester/util/SequenceIdParser.java
index b8bebe1..8def260 100644
--- a/forester/java/src/org/forester/util/SequenceIdParser.java
+++ b/forester/java/src/org/forester/util/SequenceIdParser.java
@@ -28,26 +28,107 @@
 
 package org.forester.util;
 
+import java.util.regex.Matcher;
+import java.util.regex.Pattern;
+
 import org.forester.phylogeny.data.Identifier;
+import org.forester.ws.uniprot.DatabaseTools;
 
 public final class SequenceIdParser {
 
     
-    // ref_XP_002434188_1_mites___ticks_
-    // ref_NP_001121530_1_frogs___toads_
+   
     // gb_ADF31344_1_segmented_worms_
     // gb_AAA96518_1
     // gb_EHB07727_1_rodents_
     // dbj_BAF37827_1_turtles_
     // emb_CAA73223_1_primates_
     // lcl_91970_unknown_
+    // mites|ref_XP_002434188_1
+    // ref_XP_002434188_1_mites___ticks_
+    // ref_NP_001121530_1_frogs___toads_
     
+    //The format for GenBank Accession numbers are:
+    //Nucleotide: 1 letter + 5 numerals OR 2 letters + 6 numerals
+    //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)" );
+    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)" );
+    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)" );
+    private final static boolean DEBUG                           = true;
+
     
+    /**
+     * Returns null if no match.
+     * 
+     */
     public final static Identifier parse( final String s ) {
-        
+        String v = parseGenbankAccessor( s );
+        if ( !ForesterUtil.isEmpty( v ) ) {
+            return new Identifier( v, "ncbi" );
+        }
+        v = parseRefSeqAccessor( s );
+        if ( !ForesterUtil.isEmpty( v ) ) {
+            return new Identifier( v, "ncbi" );
+        }
         return null;
     }
     
+    /**
+     * Returns null if no match.
+     * 
+     * @param query
+     * @param db 
+     * @return
+     */
+    static public String parseGenbankAccessor( final String query ) {
+        Matcher m = GENBANK_NUCLEOTIDE_AC_PATTERN_1.matcher( query );
+        if ( m.lookingAt() ) {
+            return m.group( 1 );
+        }
+        else {
+            m = GENBANK_NUCLEOTIDE_AC_PATTERN_2.matcher( query );
+            if ( m.lookingAt() ) {
+                return m.group( 1 );
+            }
+            else {
+                m = GENBANK_PROTEIN_AC_PATTERN.matcher( query );
+                if ( m.lookingAt() ) {
+                    return m.group( 1 );
+                }
+                else {
+                    return null;
+                }
+            }
+        }
+    }
+    
+    public final static String parseRefSeqAccessor( final String query ) {
+        Matcher m = GENBANK_NUCLEOTIDE_AC_PATTERN_1.matcher( query );
+        if ( m.lookingAt() ) {
+            return m.group( 1 );
+        }
+        else {
+            m = GENBANK_NUCLEOTIDE_AC_PATTERN_2.matcher( query );
+            if ( m.lookingAt() ) {
+                return m.group( 1 );
+            }
+            else {
+                m = GENBANK_PROTEIN_AC_PATTERN.matcher( query );
+                if ( m.lookingAt() ) {
+                    return m.group( 1 );
+                }
+                else {
+                    return null;
+                }
+            }
+        }
+    }
+    
+    
     
     private SequenceIdParser() {
         // Hiding the constructor.