push @args, "-v $options->{ 'mismatches' }"; # DANGER: using seed mismatches as alignment mismatches - may work, may not!
push @args, "-f" if $type eq "FASTA";
push @args, "-p $options->{ 'cpus' }";
+push @args, "--phred64-quals";
if ( defined $options->{ 'max_hits' } ) {
push @args, "-k $options->{ 'max_hits' }";
$record->{ 'S_LEN' } = length $entry->[ 4 ];
$record->{ 'SEQ_LEN' } = length $entry->[ 4 ];
$record->{ 'S_END' } = $record->{ 'S_BEG' } + $record->{ 'SEQ_LEN' } - 1;
- $record->{ 'SCORES' } =~ s/(.)/ord( $1 ) - 33 . ";"/ge; # http://maq.sourceforge.net/fastq.shtml
- $record->{ 'SCORE_MEAN' } = sprintf( "%.2f", Maasha::Calc::mean( [ split /;/, $record->{ 'SCORES' } ] ) );
+ $record->{ 'SCORES' } =~ s/(.)/chr( ( ord( $1 ) - 33 ) + 64 )/ge; # convert phred scores to illumina scores
$record->{ 'HITS' } = '.';
$record->{ 'REC_TYPE' } = "BOWTIE";