fixed strand bug in SAM.pm

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

diff --git a/code_perl/Maasha/SAM.pm b/code_perl/Maasha/SAM.pm
index be1002aa4d026258e094bd357a4077dfaa6dd6be..197a957619021f11718b923e4940e2dab4bdf657 100644 (file)
--- a/code_perl/Maasha/SAM.pm
+++ b/code_perl/Maasha/SAM.pm
@@ -130,10 +130,10 @@ sub sam2biopiece
     $record->{ 'Q_SEQ_UNMAPPED' } = $record->{ 'FLAG' } & 0x0004;
     $record->{ 'MATE_UNMAPPED' }  = $record->{ 'FLAG' } & 0x0008;
 
-    if ( $record->{ 'FLAG' } & 0x0010 ) {
-        $record->{ 'STRAND' } = '+';
-    } else {
+    if ( $record->{ 'FLAG' } & 0x0010 ) {   # 0 for forward, 1 for reverse
         $record->{ 'STRAND' } = '-';
+    } else {
+        $record->{ 'STRAND' } = '+';
     }
 
     $record->{ 'MATE_STRAND' }    = $record->{ 'FLAG' } & 0x0020;