phylotastic hackathon at NESCENT 120606

[jalview.git] / forester / java / src / org / forester / util / SequenceIdParser.java

diff --git a/forester/java/src/org/forester/util/SequenceIdParser.java b/forester/java/src/org/forester/util/SequenceIdParser.java
index b8bebe1..8def260 100644 (file)
--- a/forester/java/src/org/forester/util/SequenceIdParser.java
+++ b/forester/java/src/org/forester/util/SequenceIdParser.java
@@ -28,26 +28,107 @@
 \r
 package org.forester.util;\r
 \r
+import java.util.regex.Matcher;\r
+import java.util.regex.Pattern;\r
+\r
 import org.forester.phylogeny.data.Identifier;\r
+import org.forester.ws.uniprot.DatabaseTools;\r
 \r
 public final class SequenceIdParser {\r
 \r
     \r
-    // ref_XP_002434188_1_mites___ticks_\r
-    // ref_NP_001121530_1_frogs___toads_\r
+   \r
     // gb_ADF31344_1_segmented_worms_\r
     // gb_AAA96518_1\r
     // gb_EHB07727_1_rodents_\r
     // dbj_BAF37827_1_turtles_\r
     // emb_CAA73223_1_primates_\r
     // lcl_91970_unknown_\r
+    // mites|ref_XP_002434188_1\r
+    // ref_XP_002434188_1_mites___ticks_\r
+    // ref_NP_001121530_1_frogs___toads_\r
     \r
+    //The format for GenBank Accession numbers are:\r
+    //Nucleotide: 1 letter + 5 numerals OR 2 letters + 6 numerals\r
+    //Protein:    3 letters + 5 numerals\r
+    //http://www.ncbi.nlm.nih.gov/Sequin/acc.html\r
+    private final static Pattern GENBANK_NUCLEOTIDE_AC_PATTERN_1 = Pattern\r
+                                                                         .compile( "(?:\\A|.*[^a-zA-Z0-9])([A-Z]\\d{5})(?:[^a-zA-Z0-9]|\\Z)" );\r
+    private final static Pattern GENBANK_NUCLEOTIDE_AC_PATTERN_2 = Pattern\r
+                                                                         .compile( "(?:\\A|.*[^a-zA-Z0-9])([A-Z]{2}\\d{6})(?:[^a-zA-Z0-9]|\\Z)" );\r
+    private final static Pattern GENBANK_PROTEIN_AC_PATTERN      = Pattern\r
+                                                                         .compile( "(?:\\A|.*[^a-zA-Z0-9])([A-Z]{3}\\d{5})(?:[^a-zA-Z0-9]|\\Z)" );\r
+    private final static boolean DEBUG                           = true;\r
+\r
     \r
+    /**\r
+     * Returns null if no match.\r
+     * \r
+     */\r
     public final static Identifier parse( final String s ) {\r
-        \r
+        String v = parseGenbankAccessor( s );\r
+        if ( !ForesterUtil.isEmpty( v ) ) {\r
+            return new Identifier( v, "ncbi" );\r
+        }\r
+        v = parseRefSeqAccessor( s );\r
+        if ( !ForesterUtil.isEmpty( v ) ) {\r
+            return new Identifier( v, "ncbi" );\r
+        }\r
         return null;\r
     }\r
     \r
+    /**\r
+     * Returns null if no match.\r
+     * \r
+     * @param query\r
+     * @param db \r
+     * @return\r
+     */\r
+    static public String parseGenbankAccessor( final String query ) {\r
+        Matcher m = GENBANK_NUCLEOTIDE_AC_PATTERN_1.matcher( query );\r
+        if ( m.lookingAt() ) {\r
+            return m.group( 1 );\r
+        }\r
+        else {\r
+            m = GENBANK_NUCLEOTIDE_AC_PATTERN_2.matcher( query );\r
+            if ( m.lookingAt() ) {\r
+                return m.group( 1 );\r
+            }\r
+            else {\r
+                m = GENBANK_PROTEIN_AC_PATTERN.matcher( query );\r
+                if ( m.lookingAt() ) {\r
+                    return m.group( 1 );\r
+                }\r
+                else {\r
+                    return null;\r
+                }\r
+            }\r
+        }\r
+    }\r
+    \r
+    public final static String parseRefSeqAccessor( final String query ) {\r
+        Matcher m = GENBANK_NUCLEOTIDE_AC_PATTERN_1.matcher( query );\r
+        if ( m.lookingAt() ) {\r
+            return m.group( 1 );\r
+        }\r
+        else {\r
+            m = GENBANK_NUCLEOTIDE_AC_PATTERN_2.matcher( query );\r
+            if ( m.lookingAt() ) {\r
+                return m.group( 1 );\r
+            }\r
+            else {\r
+                m = GENBANK_PROTEIN_AC_PATTERN.matcher( query );\r
+                if ( m.lookingAt() ) {\r
+                    return m.group( 1 );\r
+                }\r
+                else {\r
+                    return null;\r
+                }\r
+            }\r
+        }\r
+    }\r
+    \r
+    \r
     \r
     private SequenceIdParser() {\r
         // Hiding the constructor.\r