added illumina 1.8 name scheme support to order pairs

[biopieces.git] / bp_bin / align_pair_seq

diff --git a/bp_bin/align_pair_seq b/bp_bin/align_pair_seq
index 24b8c659bc2523a5eee247d677d89b07bced845b..7f622f121d397cf0f124830860d3635935f94c5b 100755 (executable)
--- a/bp_bin/align_pair_seq
+++ b/bp_bin/align_pair_seq
@@ -36,7 +36,7 @@ use Data::Dumper;
 # >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
 
 
-my ( $options, $in, $out, $count, $record, @records );
+my ( $options, $in, $out, $count, $record, @records, $space, $matches );
 
 $options = Maasha::Biopieces::parse_options();
 
@@ -55,13 +55,9 @@ while ( $record = Maasha::Biopieces::get_record( $in ) )
 
         if ( $count == 2 )
         {
-            $records[ 0 ]->{ 'SEQ' } = lc $records[ 0 ]->{ 'SEQ' };
-            $records[ 1 ]->{ 'SEQ' } = lc $records[ 1 ]->{ 'SEQ' };
-
-            my $matches = Maasha::AlignTwoSeq::align_two_seq( { 'Q_SEQ' => \$records[ 0 ]->{ 'SEQ' } , 'S_SEQ' => \$records[ 1 ]->{ 'SEQ' } }, [] );
+            $matches = Maasha::AlignTwoSeq::align_two_seq( { 'Q_SEQ' => \$records[ 0 ]->{ 'SEQ' } , 'S_SEQ' => \$records[ 1 ]->{ 'SEQ' } }, [] );
             
-            print Dumper( $matches );
-
+            Maasha::AlignTwoSeq::shift_case( $matches, \$records[ 0 ]->{ 'SEQ' } ,\$records[ 1 ]->{ 'SEQ' } );
             Maasha::AlignTwoSeq::insert_indels( $matches, \$records[ 0 ]->{ 'SEQ' } ,\$records[ 1 ]->{ 'SEQ' } );
 
             map { Maasha::Biopieces::put_record( $_, $out ) } @records;