X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;ds=sidebyside;f=forester%2Fjava%2Fsrc%2Forg%2Fforester%2Fws%2Fseqdb%2FSequenceDbWsTools.java;h=9106c38aae62aa1472ad384338ce31384a409a6c;hb=10297bd8b8a4b4ab198a17a42fc6ff24ae2ed49b;hp=1bb2e986ee39bbc6d00d74b280f79308704853df;hpb=cbc5c71b164a57b8ad6c988d015057c7f0972478;p=jalview.git diff --git a/forester/java/src/org/forester/ws/seqdb/SequenceDbWsTools.java b/forester/java/src/org/forester/ws/seqdb/SequenceDbWsTools.java index 1bb2e98..9106c38 100644 --- a/forester/java/src/org/forester/ws/seqdb/SequenceDbWsTools.java +++ b/forester/java/src/org/forester/ws/seqdb/SequenceDbWsTools.java @@ -26,6 +26,7 @@ package org.forester.ws.seqdb; import java.io.BufferedReader; +import java.io.FileNotFoundException; import java.io.IOException; import java.io.InputStreamReader; import java.io.UnsupportedEncodingException; @@ -37,38 +38,38 @@ import java.util.List; import java.util.SortedSet; import java.util.TreeSet; +import org.forester.go.GoTerm; +import org.forester.io.parsers.phyloxml.PhyloXmlDataFormatException; import org.forester.phylogeny.Phylogeny; import org.forester.phylogeny.PhylogenyNode; import org.forester.phylogeny.data.Accession; +import org.forester.phylogeny.data.Accession.Source; +import org.forester.phylogeny.data.Annotation; import org.forester.phylogeny.data.Identifier; import org.forester.phylogeny.data.Sequence; import org.forester.phylogeny.data.Taxonomy; import org.forester.phylogeny.iterators.PhylogenyNodeIterator; +import org.forester.sequence.MolecularSequence.TYPE; import org.forester.util.ForesterUtil; -import org.forester.util.SequenceIdParser; +import org.forester.util.SequenceAccessionTools; public final class SequenceDbWsTools { - public final static String BASE_UNIPROT_URL = "http://www.uniprot.org/"; - public final static String BASE_EMBL_DB_URL = "http://www.ebi.ac.uk/Tools/dbfetch/dbfetch/"; - public final static String EMBL_DBS_EMBL = "embl"; - public final static String EMBL_DBS_REFSEQ_P = "refseqp"; - public final static String EMBL_DBS_REFSEQ_N = "refseqn"; - private final static String URL_ENC = "UTF-8"; - private final static boolean DEBUG = false; - - private static List getTaxonomiesFromCommonName( final String cn, final int max_taxonomies_return ) - throws IOException { - final List result = getTaxonomyStringFromCommonName( cn, max_taxonomies_return ); - if ( result.size() > 0 ) { - return parseUniProtTaxonomy( result ); - } - return null; - } + public final static String BASE_UNIPROT_URL = "http://www.uniprot.org/"; + public final static int DEFAULT_LINES_TO_RETURN = 4000; + public final static String EMBL_DBS_REFSEQ_N = "refseqn"; + public final static String EMBL_DBS_REFSEQ_P = "refseqp"; + public final static String EMBL_GENBANK = "http://www.ebi.ac.uk/Tools/dbfetch/dbfetch?db=GENBANK&style=raw&id="; + public final static String EMBL_REFSEQ = "http://www.ebi.ac.uk/Tools/dbfetch/dbfetch?db=REFSEQ&style=raw&id="; + public final static String EMBL_EMBL = "http://www.ebi.ac.uk/Tools/dbfetch/dbfetch?db=EMBL&style=raw&id="; + private final static boolean DEBUG = true; + private final static String URL_ENC = "UTF-8"; + private final static int SLEEP = 200; + private static final boolean ALLOW_TO_OVERWRITE_MOL_SEQ = false; public static List getTaxonomiesFromCommonNameStrict( final String cn, final int max_taxonomies_return ) - throws IOException { + throws IOException { final List taxonomies = getTaxonomiesFromCommonName( cn, max_taxonomies_return ); if ( ( taxonomies != null ) && ( taxonomies.size() > 0 ) ) { final List filtered_taxonomies = new ArrayList(); @@ -91,25 +92,15 @@ public final class SequenceDbWsTools { return null; } - private static List getTaxonomiesFromScientificName( final String sn, - final int max_taxonomies_return ) - throws IOException { - final List result = getTaxonomyStringFromScientificName( sn, max_taxonomies_return ); - if ( result.size() > 0 ) { - return parseUniProtTaxonomy( result ); - } - return null; - } - /** * Does not return "sub-types". * For example, for "Mus musculus" only returns "Mus musculus" * and not "Mus musculus", "Mus musculus bactrianus", ... - * + * */ public static List getTaxonomiesFromScientificNameStrict( final String sn, final int max_taxonomies_return ) - throws IOException { + throws IOException { final List taxonomies = getTaxonomiesFromScientificName( sn, max_taxonomies_return ); if ( ( taxonomies != null ) && ( taxonomies.size() > 0 ) ) { final List filtered_taxonomies = new ArrayList(); @@ -125,7 +116,7 @@ public final class SequenceDbWsTools { public static List getTaxonomiesFromTaxonomyCode( final String code, final int max_taxonomies_return ) - throws IOException { + throws IOException { final String my_code = new String( code ); final List result = getTaxonomyStringFromTaxonomyCode( my_code, max_taxonomies_return ); if ( result.size() > 0 ) { @@ -134,135 +125,99 @@ public final class SequenceDbWsTools { return null; } - public static SequenceDatabaseEntry obtainEmblEntry( final Identifier id, final int max_lines_to_return ) + public static SequenceDatabaseEntry obtainEmblEntry( final Accession acc ) throws IOException { + return obtainEmblEntry( acc, DEFAULT_LINES_TO_RETURN ); + } + + public static SequenceDatabaseEntry obtainEmblEntry( final Accession acc, final int max_lines_to_return ) throws IOException { - final List lines = queryEmblDb( id, max_lines_to_return ); - return EbiDbEntry.createInstanceFromPlainText( lines ); + final List lines = queryEmblDb( acc, max_lines_to_return ); + return EbiDbEntry.createInstanceFromPlainTextForRefSeq( lines ); } - public static SequenceDatabaseEntry obtainRefSeqEntryFromEmbl( final Identifier id, final int max_lines_to_return ) + public static SequenceDatabaseEntry obtainEntry( final String acc_str ) throws IOException { + if ( ForesterUtil.isEmpty( acc_str ) ) { + throw new IllegalArgumentException( "cannot not extract sequence db accessor from null or empty string" ); + } + final Accession acc = SequenceAccessionTools.parseAccessorFromString( acc_str ); + if ( acc == null ) { + throw new IllegalArgumentException( "could not extract acceptable sequence db accessor from \"" + acc_str + + "\"" ); + } + if ( acc.getSource().equals( Source.REFSEQ.toString() ) || acc.getSource().equals( Source.EMBL.toString() ) + || acc.getSource().equals( Source.NCBI.toString() ) ) { + return obtainEmblEntry( acc, DEFAULT_LINES_TO_RETURN ); + } + else if ( acc.getSource().equals( Source.UNIPROT.toString() ) ) { + return obtainUniProtEntry( acc.getValue(), DEFAULT_LINES_TO_RETURN ); + } + else { + throw new IllegalArgumentException( "don't know how to handle request for source \"" + acc.getSource() + + "\"" ); + } + } + + public static SequenceDatabaseEntry obtainRefSeqEntryFromEmbl( final Accession acc ) throws IOException { + return obtainRefSeqEntryFromEmbl( acc, DEFAULT_LINES_TO_RETURN ); + } + + public static SequenceDatabaseEntry obtainRefSeqEntryFromEmbl( final Accession acc, final int max_lines_to_return ) throws IOException { - final List lines = queryEmblDb( id, max_lines_to_return ); + final List lines = queryEmblDbForRefSeqEntry( acc, max_lines_to_return ); return EbiDbEntry.createInstanceFromPlainTextForRefSeq( lines ); } - public static SortedSet obtainSeqInformation( final Phylogeny phy, - final boolean ext_nodes_only, - final boolean allow_to_set_taxonomic_data, - final int lines_to_return ) throws IOException { + public final static Accession obtainSeqAccession( final PhylogenyNode node ) { + Accession acc = SequenceAccessionTools.obtainFromSeqAccession( node ); + if ( !isAccessionAcceptable( acc ) ) { + acc = SequenceAccessionTools.obtainAccessorFromDataFields( node ); + } + return acc; + } + + public final static void obtainSeqInformation( final boolean allow_to_set_taxonomic_data, + final int lines_to_return, + final SortedSet not_found, + final PhylogenyNode node ) throws IOException { + final Accession acc = obtainSeqAccession( node ); + if ( !isAccessionAcceptable( acc ) ) { + if ( node.isExternal() || !node.isEmpty() ) { + not_found.add( node.toString() ); + } + } + else { + addDataFromDbToNode( allow_to_set_taxonomic_data, lines_to_return, not_found, node, acc ); + } + } + + public final static void obtainSeqInformation( final boolean allow_to_set_taxonomic_data, + final SortedSet not_found, + final PhylogenyNode node ) throws IOException { + obtainSeqInformation( allow_to_set_taxonomic_data, DEFAULT_LINES_TO_RETURN, not_found, node ); + } + + public final static SortedSet obtainSeqInformation( final Phylogeny phy, + final boolean ext_nodes_only, + final boolean allow_to_set_taxonomic_data, + final int lines_to_return ) throws IOException { final SortedSet not_found = new TreeSet(); for( final PhylogenyNodeIterator iter = phy.iteratorPostorder(); iter.hasNext(); ) { final PhylogenyNode node = iter.next(); - if ( ext_nodes_only && node.isInternal() ) { - continue; - } - String query = null; - Identifier id = null; - Db db = Db.NONE; - if ( node.getNodeData().isHasSequence() && ( node.getNodeData().getSequence().getAccession() != null ) - && !ForesterUtil.isEmpty( node.getNodeData().getSequence().getAccession().getSource() ) - && !ForesterUtil.isEmpty( node.getNodeData().getSequence().getAccession().getValue() ) - && node.getNodeData().getSequence().getAccession().getValue().toLowerCase().startsWith( "uniprot" ) ) { - query = node.getNodeData().getSequence().getAccession().getValue(); - db = Db.UNIPROT; - } - else if ( node.getNodeData().isHasSequence() - && ( node.getNodeData().getSequence().getAccession() != null ) - && !ForesterUtil.isEmpty( node.getNodeData().getSequence().getAccession().getSource() ) - && !ForesterUtil.isEmpty( node.getNodeData().getSequence().getAccession().getValue() ) - && ( node.getNodeData().getSequence().getAccession().getValue().toLowerCase().startsWith( "embl" ) || node - .getNodeData().getSequence().getAccession().getValue().toLowerCase().startsWith( "ebi" ) ) ) { - query = node.getNodeData().getSequence().getAccession().getValue(); - db = Db.EMBL; - } - else if ( !ForesterUtil.isEmpty( node.getName() ) ) { - if ( ( query = ForesterUtil.extractUniProtKbProteinSeqIdentifier( node ) ) != null ) { - db = Db.UNIPROT; - } - else if ( ( id = SequenceIdParser.parse( node.getName() ) ) != null ) { - if ( id.getProvider().equalsIgnoreCase( Identifier.NCBI ) ) { - db = Db.NCBI; - } - else if ( id.getProvider().equalsIgnoreCase( Identifier.REFSEQ ) ) { - db = Db.REFSEQ; - } - } - } - if ( db == Db.NONE ) { - not_found.add( node.getName() ); - } - SequenceDatabaseEntry db_entry = null; - if ( !ForesterUtil.isEmpty( query ) ) { - if ( db == Db.UNIPROT ) { - if ( DEBUG ) { - System.out.println( "uniprot: " + query ); - } - db_entry = obtainUniProtEntry( query, lines_to_return ); - } - if ( ( db == Db.EMBL ) || ( ( db == Db.UNIPROT ) && ( db_entry == null ) ) ) { - if ( DEBUG ) { - System.out.println( "embl: " + query ); - } - db_entry = obtainEmblEntry( new Identifier( query ), lines_to_return ); - if ( ( db == Db.UNIPROT ) && ( db_entry != null ) ) { - db = Db.EMBL; - } - } - } - else if ( ( db == Db.REFSEQ ) && ( id != null ) ) { - db_entry = obtainRefSeqEntryFromEmbl( id, lines_to_return ); - } - else if ( ( db == Db.NCBI ) && ( id != null ) ) { - db_entry = obtainEmblEntry( id, lines_to_return ); //TODO ? - } - if ( ( db_entry != null ) && !db_entry.isEmpty() ) { - final Sequence seq = node.getNodeData().isHasSequence() ? node.getNodeData().getSequence() - : new Sequence(); - if ( !ForesterUtil.isEmpty( db_entry.getAccession() ) ) { - String type = null; - if ( db == Db.EMBL ) { - type = "embl"; - } - else if ( db == Db.UNIPROT ) { - type = "uniprot"; - } - else if ( db == Db.NCBI ) { - type = "ncbi"; - } - else if ( db == Db.REFSEQ ) { - type = "refseq"; - } - seq.setAccession( new Accession( db_entry.getAccession(), type ) ); - } - if ( !ForesterUtil.isEmpty( db_entry.getSequenceName() ) ) { - seq.setName( db_entry.getSequenceName() ); - } - if ( !ForesterUtil.isEmpty( db_entry.getSequenceSymbol() ) ) { - seq.setSymbol( db_entry.getSequenceSymbol() ); - } - final Taxonomy tax = node.getNodeData().isHasTaxonomy() ? node.getNodeData().getTaxonomy() - : new Taxonomy(); - if ( !ForesterUtil.isEmpty( db_entry.getTaxonomyScientificName() ) ) { - tax.setScientificName( db_entry.getTaxonomyScientificName() ); - } - if ( allow_to_set_taxonomic_data && !ForesterUtil.isEmpty( db_entry.getTaxonomyIdentifier() ) ) { - tax.setIdentifier( new Identifier( db_entry.getTaxonomyIdentifier(), "uniprot" ) ); - } - node.getNodeData().setTaxonomy( tax ); - node.getNodeData().setSequence( seq ); - } - else if ( db != Db.NONE ) { - not_found.add( node.getName() ); - } - try { - Thread.sleep( 10 );// Sleep for 10 ms - } - catch ( final InterruptedException ie ) { + if ( node.isExternal() || !ext_nodes_only ) { + obtainSeqInformation( allow_to_set_taxonomic_data, lines_to_return, not_found, node ); } } return not_found; } + public final static void obtainSeqInformation( final PhylogenyNode node ) throws IOException { + obtainSeqInformation( true, DEFAULT_LINES_TO_RETURN, new TreeSet(), node ); + } + + public static SequenceDatabaseEntry obtainUniProtEntry( final String query ) throws IOException { + return obtainUniProtEntry( query, DEFAULT_LINES_TO_RETURN ); + } + public static SequenceDatabaseEntry obtainUniProtEntry( final String query, final int max_lines_to_return ) throws IOException { final List lines = queryUniprot( "uniprot/" + query + ".txt", max_lines_to_return ); @@ -296,8 +251,8 @@ public final class SequenceDbWsTools { } in.close(); try { - // To prevent accessing online dbs in too quick succession. - Thread.sleep( 20 ); + // To prevent accessing online dbs in too quick succession. + Thread.sleep( SLEEP ); } catch ( final InterruptedException e ) { e.printStackTrace(); @@ -305,23 +260,30 @@ public final class SequenceDbWsTools { return result; } - public static List queryEmblDb( final Identifier id, final int max_lines_to_return ) throws IOException { + public static List queryEmblDb( final Accession acc, final int max_lines_to_return ) throws IOException { final StringBuilder url_sb = new StringBuilder(); - url_sb.append( BASE_EMBL_DB_URL ); - if ( ForesterUtil.isEmpty( id.getProvider() ) || id.getProvider().equalsIgnoreCase( Identifier.NCBI ) ) { - url_sb.append( SequenceDbWsTools.EMBL_DBS_EMBL ); - url_sb.append( '/' ); + // url_sb.append( BASE_EMBL_DB_URL ); + System.out.println( "source: " + acc.getSource() ); + if ( acc.getSource().equals( Source.NCBI.toString() ) ) { + url_sb.append( EMBL_GENBANK ); + //url_sb.append( '/' ); } - else if ( id.getProvider().equalsIgnoreCase( Identifier.REFSEQ ) ) { - if ( id.getValue().toUpperCase().indexOf( 'P' ) == 1 ) { - url_sb.append( SequenceDbWsTools.EMBL_DBS_REFSEQ_P ); - url_sb.append( '/' ); - } - else { - url_sb.append( SequenceDbWsTools.EMBL_DBS_REFSEQ_N ); - url_sb.append( '/' ); - } + else if ( acc.getSource().equals( Source.REFSEQ.toString() ) ) { + url_sb.append( EMBL_REFSEQ ); + } + else if ( acc.getSource().equals( Source.EMBL.toString() ) ) { + url_sb.append( EMBL_EMBL ); } + else { + throw new IllegalArgumentException( "unable to handle source: " + acc.getSource() ); + } + return queryDb( acc.getValue(), max_lines_to_return, url_sb.toString() ); + } + + public static List queryEmblDbForRefSeqEntry( final Accession id, final int max_lines_to_return ) + throws IOException { + final StringBuilder url_sb = new StringBuilder(); + url_sb.append( EMBL_REFSEQ ); return queryDb( id.getValue(), max_lines_to_return, url_sb.toString() ); } @@ -329,10 +291,178 @@ public final class SequenceDbWsTools { return queryDb( query, max_lines_to_return, BASE_UNIPROT_URL ); } + final static String extractFrom( final String target, final String a ) { + final int i_a = target.indexOf( a ); + return target.substring( i_a + a.length() ).trim(); + } + + final static String extractFromTo( final String target, final String a, final String b ) { + final int i_a = target.indexOf( a ); + final int i_b = target.indexOf( b ); + if ( ( i_a < 0 ) || ( i_b < i_a ) ) { + throw new IllegalArgumentException( "attempt to extract from \"" + target + "\" between \"" + a + + "\" and \"" + b + "\"" ); + } + return target.substring( i_a + a.length(), i_b ).trim(); + } + + final static String extractTo( final String target, final String b ) { + final int i_b = target.indexOf( b ); + return target.substring( 0, i_b ).trim(); + } + + private static void addDataFromDbToNode( final boolean allow_to_set_taxonomic_data, + final int lines_to_return, + final SortedSet not_found, + final PhylogenyNode node, + final Accession acc ) throws IOException { + SequenceDatabaseEntry db_entry = null; + final String query = acc.getValue(); + if ( acc.getSource().equals( Source.UNIPROT.toString() ) ) { + if ( DEBUG ) { + System.out.println( "uniprot: " + query ); + } + try { + db_entry = obtainUniProtEntry( query, lines_to_return ); + } + catch ( final FileNotFoundException e ) { + // Eat this, and move to next. + } + } + else if ( acc.getSource().equals( Source.REFSEQ.toString() ) ) { + if ( DEBUG ) { + System.out.println( "refseq: " + query ); + } + try { + db_entry = obtainRefSeqEntryFromEmbl( new Accession( query ), lines_to_return ); + } + catch ( final FileNotFoundException e ) { + // Eat this, and move to next. + } + } + else if ( acc.getSource().equals( Source.EMBL.toString() ) || acc.getSource().equals( Source.NCBI.toString() ) + || acc.getSource().equals( Source.EMBL.toString() ) ) { + if ( DEBUG ) { + System.out.println( acc.toString() ); + } + try { + db_entry = obtainEmblEntry( acc, lines_to_return ); + } + catch ( final FileNotFoundException e ) { + // Eat this, and move to next. + } + } + else if ( acc.getSource().equals( Source.GI.toString() ) ) { + if ( DEBUG ) { + System.out.println( "gi: " + query ); + } + try { + db_entry = obtainRefSeqEntryFromEmbl( new Accession( query ), lines_to_return ); + } + catch ( final FileNotFoundException e ) { + // Eat this, and move to next. + } + } + if ( ( db_entry != null ) && !db_entry.isEmpty() ) { + final Sequence seq = node.getNodeData().isHasSequence() ? node.getNodeData().getSequence() : new Sequence(); + if ( !ForesterUtil.isEmpty( db_entry.getAccession() ) ) { + seq.setAccession( new Accession( db_entry.getAccession(), acc.getSource() ) ); + } + if ( !ForesterUtil.isEmpty( db_entry.getSequenceName() ) ) { + seq.setName( db_entry.getSequenceName() ); + } + if ( !ForesterUtil.isEmpty( db_entry.getGeneName() ) ) { + seq.setGeneName( db_entry.getGeneName() ); + } + if ( !ForesterUtil.isEmpty( db_entry.getSequenceSymbol() ) ) { + try { + seq.setSymbol( db_entry.getSequenceSymbol() ); + } + catch ( final PhyloXmlDataFormatException e ) { + // Eat this exception. + } + } + if ( ( db_entry.getMolecularSequence() != null ) + && !ForesterUtil.isEmpty( db_entry.getMolecularSequence().getMolecularSequenceAsString() ) + && ( ALLOW_TO_OVERWRITE_MOL_SEQ || seq.getMolecularSequence().isEmpty() ) ) { + seq.setMolecularSequence( db_entry.getMolecularSequence().getMolecularSequenceAsString() ); + seq.setMolecularSequenceAligned( false ); + if ( db_entry.getMolecularSequence().getType() == TYPE.AA ) { + seq.setType( "protein" ); + } + else if ( db_entry.getMolecularSequence().getType() == TYPE.DNA ) { + seq.setType( "dna" ); + } + else if ( db_entry.getMolecularSequence().getType() == TYPE.RNA ) { + seq.setType( "rna" ); + } + } + if ( ( db_entry.getGoTerms() != null ) && !db_entry.getGoTerms().isEmpty() ) { + for( final GoTerm go : db_entry.getGoTerms() ) { + final Annotation ann = new Annotation( go.getGoId().getId() ); + ann.setDesc( go.getName() ); + seq.addAnnotation( ann ); + } + } + if ( ( db_entry.getCrossReferences() != null ) && !db_entry.getCrossReferences().isEmpty() ) { + for( final Accession x : db_entry.getCrossReferences() ) { + seq.addCrossReference( x ); + } + } + if ( !ForesterUtil.isEmpty( db_entry.getChromosome() ) && !ForesterUtil.isEmpty( db_entry.getMap() ) ) { + seq.setLocation( "chr " + db_entry.getChromosome() + ", " + db_entry.getMap() ); + } + else if ( !ForesterUtil.isEmpty( db_entry.getChromosome() ) ) { + seq.setLocation( "chr " + db_entry.getChromosome() ); + } + else if ( !ForesterUtil.isEmpty( db_entry.getMap() ) ) { + seq.setLocation( db_entry.getMap() ); + } + final Taxonomy tax = node.getNodeData().isHasTaxonomy() ? node.getNodeData().getTaxonomy() : new Taxonomy(); + if ( !ForesterUtil.isEmpty( db_entry.getTaxonomyScientificName() ) ) { + tax.setScientificName( db_entry.getTaxonomyScientificName() ); + } + if ( allow_to_set_taxonomic_data && !ForesterUtil.isEmpty( db_entry.getTaxonomyIdentifier() ) ) { + tax.setIdentifier( new Identifier( db_entry.getTaxonomyIdentifier(), "uniprot" ) ); + } + node.getNodeData().setTaxonomy( tax ); + node.getNodeData().setSequence( seq ); + } + else { + if ( node.isExternal() || !node.isEmpty() ) { + not_found.add( node.toString() ); + } + } + try { + Thread.sleep( SLEEP ); + } + catch ( final InterruptedException ie ) { + } + } + private static String encode( final String str ) throws UnsupportedEncodingException { return URLEncoder.encode( str.trim(), URL_ENC ); } + private static List getTaxonomiesFromCommonName( final String cn, final int max_taxonomies_return ) + throws IOException { + final List result = getTaxonomyStringFromCommonName( cn, max_taxonomies_return ); + if ( result.size() > 0 ) { + return parseUniProtTaxonomy( result ); + } + return null; + } + + private static List getTaxonomiesFromScientificName( final String sn, + final int max_taxonomies_return ) + throws IOException { + final List result = getTaxonomyStringFromScientificName( sn, max_taxonomies_return ); + if ( result.size() > 0 ) { + return parseUniProtTaxonomy( result ); + } + return null; + } + private static List getTaxonomyStringFromCommonName( final String cn, final int max_lines_to_return ) throws IOException { return queryUniprot( "taxonomy/?query=common%3a%22" + encode( cn ) + "%22&format=tab", max_lines_to_return ); @@ -353,6 +483,13 @@ public final class SequenceDbWsTools { return queryUniprot( "taxonomy/?query=mnemonic%3a%22" + encode( code ) + "%22&format=tab", max_lines_to_return ); } + private final static boolean isAccessionAcceptable( final Accession acc ) { + return ( !( ( acc == null ) || ForesterUtil.isEmpty( acc.getSource() ) || ForesterUtil.isEmpty( acc.getValue() ) || ( ( acc + .getSource().equals( Source.UNIPROT.toString() ) ) + && ( acc.getSource().toString().equals( Source.EMBL.toString() ) ) && ( acc.getSource().toString() + .equals( Source.REFSEQ.toString() ) ) ) ) ); + } + private static List parseUniProtTaxonomy( final List result ) throws IOException { final List taxonomies = new ArrayList(); for( final String line : result ) { @@ -376,8 +513,4 @@ public final class SequenceDbWsTools { } return taxonomies; } - - public enum Db { - UNIPROT, EMBL, NCBI, NONE, REFSEQ; - } }