- or download this
my @CONV = glob( "{T,C,A,G}" x 4 );
my $dna = join "", @CONV[ unpack "C*", $raw ];
- or download this
## takes 16 bits (= 8 bases = unsigned short) at a time
my @CONV = glob( "{T,C,A,G}" x 8 );
my $dna = join "", @CONV[ unpack "S*", $raw ];
- or download this
substr($dna, $_, $mblock{$_}, lc(substr($dna, $_, $mblock{$_})))
- or download this
substr($dna, $_, $mblock{$_}) ^= (" " x $mblock{$_});