X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=code_perl%2FMaasha%2FFasta.pm;fp=code_perl%2FMaasha%2FFasta.pm;h=c4920a805c8805d50390c2804b22e528c11069df;hb=430634c816f774703a59c0122006daf70ee570f8;hp=cf399cf16ba07854f230ca67f0a8e285d53dcd79;hpb=15502970edc329bf6f14138be809733d8a8f16a2;p=biopieces.git

diff --git a/code_perl/Maasha/Fasta.pm b/code_perl/Maasha/Fasta.pm
index cf399cf..c4920a8 100644
--- a/code_perl/Maasha/Fasta.pm
+++ b/code_perl/Maasha/Fasta.pm
@@ -488,4 +488,30 @@ sub index_retrieve
 }
 
 
+sub biopiece2fasta
+{
+    # Martin A. Hansen, July 2008.
+
+    # Given a biopiece record converts it to a FASTA record.
+    # If no generic SEQ or SEQ_NAME is found, the Q_* and S_* are
+    # tried in that order.
+
+    my ( $record,    # record
+       ) = @_;
+
+    # Returns a tuple.
+
+    my ( $seq_name, $seq );
+
+    $seq_name = $record->{ "SEQ_NAME" } || $record->{ "Q_ID" }  || $record->{ "S_ID" };
+    $seq      = $record->{ "SEQ" }      || $record->{ "Q_SEQ" } || $record->{ "S_SEQ" };
+
+    if ( defined $seq_name and defined $seq ) {
+        return wantarray ? ( $seq_name, $seq ) : [ $seq_name, $seq ];
+    } else {
+        return;
+    }
+}
+
+
 # >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<