%new_record = %{ $record };
$new_record{ "SEQ" } = Maasha::Seq::translate( $record->{ "SEQ" }, $frame );
- $new_record{ "SEQ_LEN" } = length $record->{ "SEQ" };
+ $new_record{ "SEQ_LEN" } = length $new_record{ "SEQ" };
$new_record{ "FRAME" } = $frame;
put_record( \%new_record, $out );
}
}
+ $record->{ "SEQ_LEN" } = length $record->{ "SEQ" };
+
put_record( $record, $out );
}
}
if ( $options->{ "uniq" } and not $lookup{ $subseq } )
{
- $new_record->{ "REC_TYPE" } = "SPLIT";
$new_record->{ "SEQ_NAME" } = $record->{ "SEQ_NAME" } . "[" . ( $i + 1 ) . "-" . ( $i + $options->{ "word_size" } ) . "]";
$new_record->{ "SEQ" } = $subseq;
}
else
{
- $new_record->{ "REC_TYPE" } = "SPLIT";
$new_record->{ "SEQ_NAME" } = $record->{ "SEQ_NAME" } . "[" . ( $i + 1 ) . "-" . ( $i + $options->{ "word_size" } ) . "]";
$new_record->{ "SEQ" } = $subseq;
{
if ( $record->{ "SEQ_NAME" } and $record->{ "SEQ" } )
{
- $seq_cs = Maasha::Solid::seq2color_space( $record->{ "SEQ" } );
+ $seq_cs = Maasha::Solid::seq2color_space( uc $record->{ "SEQ" } );
Maasha::Fasta::put_entry( [ $record->{ "SEQ_NAME" }, $seq_cs ], $fh, $options->{ "wrap" } );
}
foreach $key ( keys %{ $analysis } )
{
$analysis->{ $key }->{ "MEAN" } = sprintf "%.2f", $analysis->{ $key }->{ "SUM" } / $analysis->{ $key }->{ "COUNT" };
- $analysis->{ $key }->{ "SUM" } = sprintf "%.2f", $analysis->{ $key }->{ "SUN" };
+ $analysis->{ $key }->{ "SUM" } = sprintf "%.2f", $analysis->{ $key }->{ "SUM" };
}
my ( $keys, $types, $counts, $mins, $maxs, $sums, $means );