# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
+use Maasha::Calc;
+use Data::Dumper;
+
use vars qw( @ISA @EXPORT_OK );
require Exporter;
# Returns a list.
- my ( $seq_header, $seq, $score_head, $score, @scores );
+ my ( $seq_header, $seq, $score_head, $score );
$seq_header = <$fh>;
$seq = <$fh>;
chomp $score_header;
chomp $score;
- @scores = split //, $score;
-
- map { $_ = int( score_oct2dec( $_ ) ) } @scores;
-
$seq_header =~ s/^@//;
-
- $score = join( ";", @scores );
+ $score =~ s/(.)/int( score_oct2dec( $1 ) ) . ";"/ge;
return wantarray ? ( $seq_header, $seq, $score ) : [ $seq_header, $seq, $score ];
}
# Returns a float
- return 10 * log(1 + 10 ** ((ord( $char ) - 64) / 10.0)) / log(10);
+ return 10 * log( 1 + 10 ** ( ( ord( $char ) - 64 ) / 10.0 ) ) / log( 10 );
}
# Returns a list.
- my ( $seq_header, $seq, $score_head, $score, @scores );
+ my ( $seq_header, $seq, $score_head, $score );
$seq_header = <$fh>;
$seq = <$fh>;
chomp $score_header;
chomp $score;
- @scores = split / /, $score;
-
$seq_header =~ s/^@//;
-
- $score = join( ";", @scores );
+ $score =~ s/ /;/g;
return wantarray ? ( $seq_header, $seq, $score ) : [ $seq_header, $seq, $score ];
}
# Converts a Solexa entry to a Biopiece record.
- my ( $entry, # Solexa entry
- $qualilty, # Quality cutoff
+ my ( $entry, # Solexa entry
+ $quality, # Quality cutoff
) = @_;
# Returns a hashref.