- or download this
sub scan_seq {
local our %sub_seqs;
$_[0] =~ /(.{7})(?{ ++$sub_seqs{$1} })(?!)/;
return \%sub_seqs;
}
- or download this
use strict;
use warnings;
...
printf($fmt_b, $seq, map $_||0, @$row{@cols});
}
printf($fmt_b, 'Total', map $_||0, @cols{@cols});
- or download this
ATGATCT ATGATGA ATGCGTA ATGCTAG ATGCTGT CGTATGC CTA
+GACT CTGTACT GATGATC GCGTATG GCTAGAC GCTGTAC GTATGCA TAGACTG
+TGATCTG TGATGAT TGCGTAT TGCTAGA TGCTGTA TGTACTG
ATGCTGTACTG 0 0 0 0 1 0
+ 0 1 0 0 0 1 0 0
+ 0 0 0 0 1 1
...
ATGCGTATGCA 0 0 1 0 0 1
+ 0 0 0 1 0 0 1 0
+ 0 0 1 0 0 0
ATGCTAGACTG 0 0 0 1 0 0
+ 1 0 0 0 1 0 0 1
+ 0 0 0 1 0 0
Total 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1
- or download this
my @cols = sort keys %cols;
- or download this
my @cols = glob('{A,C,G,T}' x $sub_seq_len);