X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=forester%2Fjava%2Fsrc%2Forg%2Fforester%2Fapplication%2Faa.java;h=b835e52cb907670adccfad1b4d1c1c538587fa1e;hb=f47b5de663d2ca59a4d53587ea8d5cd20c594db2;hp=6d7be4f6eac52eed9627736b1b4e899480225bd0;hpb=cfda6a6b1a23304bdcdac75f629e9ce4beb9cd9a;p=jalview.git

diff --git a/forester/java/src/org/forester/application/aa.java b/forester/java/src/org/forester/application/aa.java
index 6d7be4f..b835e52 100644
--- a/forester/java/src/org/forester/application/aa.java
+++ b/forester/java/src/org/forester/application/aa.java
@@ -1,15 +1,17 @@
+//
 
 package org.forester.application;
 
 import java.io.FileInputStream;
 import java.util.ArrayList;
+import java.util.Arrays;
 import java.util.HashSet;
 import java.util.List;
 import java.util.Set;
-import java.util.TreeSet;
 
 import org.forester.io.parsers.FastaParser;
 import org.forester.msa.Msa;
+import org.forester.sequence.BasicSequence;
 import org.forester.sequence.Sequence;
 import org.forester.util.ForesterUtil;
 
@@ -39,12 +41,18 @@ public class aa {
                 final List<Sequence> found_seqs = new ArrayList<Sequence>();
                 for( final Sequence orig_seq : orig ) {
                     final String orig_seq_id = orig_seq.getIdentifier();
-                    if ( orig_seq_id.indexOf( id_ ) >= 0 && orig_seq_id.indexOf( "[" + range + "]" ) >= 0 ) {
+                    if ( ( orig_seq_id.indexOf( id_ ) >= 0 ) && ( orig_seq_id.indexOf( "[" + range + "]" ) >= 0 ) ) {
                         found++;
                         found_seqs.add( orig_seq );
                     }
                 }
                 if ( found > 0 ) {
+                    for( final Sequence found_seq : found_seqs ) {
+                        if ( found_seq.getLength() >= 85 ) {
+                            all_found_seqs.add( BasicSequence.createAaSequence( id, found_seq
+                                    .getMolecularSequenceAsString() ) );
+                        }
+                    }
                     if ( found > 1 ) {
                         System.out.println( i + ": " + id + "=>" + id_ + " " + range );
                         System.out.println( "  found: " + found );
@@ -58,6 +66,17 @@ public class aa {
                     System.exit( -1 );
                 }
             }
+            final String fasta_ary[] = new String[ all_found_seqs.size() ];
+            int i = 0;
+            for( final Sequence sequence : all_found_seqs ) {
+                fasta_ary[ i ] = ">" + sequence.getIdentifier() + "\n" + sequence.getMolecularSequenceAsString();
+                System.out.println( sequence );
+                i++;
+            }
+            Arrays.sort( fasta_ary );
+            for( final String element : fasta_ary ) {
+                System.out.println( element );
+            }
             System.out.println( "DONE." );
         }
         catch ( final Exception e ) {