chomp $line;
@fields = split /\t/, $line;
-
$record = bowtie2biopiece( \@fields );
Maasha::Biopieces::put_record( $record, $out );
$record->{ 'SCORES' } = $entry->[ 5 ];
$record->{ 'SCORE' } = $entry->[ 6 ] + 1;
$record->{ 'DESCRIPTOR' } = $entry->[ 7 ];
-
$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