Help for this page

Select Code to Download


  1. or download this
    sub scan_seq {
        local our %sub_seqs;
        $_[0] =~ /(.{7})(?{ ++$sub_seqs{$1} })(?!)/;
        return \%sub_seqs;
    }
    
  2. or download this
    use strict;
    use warnings;
    ...
        printf($fmt_b, $seq, map $_||0, @$row{@cols});
    }
    printf($fmt_b, 'Total', map $_||0, @cols{@cols});
    
  3. 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
    
  4. or download this
    my @cols = sort keys %cols;
    
  5. or download this
    my @cols = glob('{A,C,G,T}' x $sub_seq_len);