fixed STRAND in SAM format

git-svn-id: http://biopieces.googlecode.com/svn/trunk@687 74ccb610-7750-0410-82ae-013aeee3265d

diff --git a/code_perl/Maasha/SAM.pm b/code_perl/Maasha/SAM.pm
index c9008674005e46890794d9b11ce0325d6678ec71..79f697d478ca77494334c38f94cd0b10aa35274b 100644 (file)
--- a/code_perl/Maasha/SAM.pm
+++ b/code_perl/Maasha/SAM.pm
@@ -116,6 +116,12 @@ sub sam2biopiece
     $record->{ 'S_BEG' } -= 1 if $record->{ 'S_BEG' }  != 0;
     $record->{ 'S_BEG' } -= 1 if $record->{ 'S_BEG2' } != 0;
 
+    if ( $record->{ 'FLAG' } & 0x0010 ) {
+        $record->{ 'STRAND' } = '+';
+    } else {
+        $record->{ 'STRAND' } = '-';
+    }
+
     return wantarray ? %{ $record } : $record;
 }