in reply to Counting overlapping dimers for multiple sequences
Code snippet below should give you some idea
Update: Fixed typo. Anyhow post seems similar to of flexvault :)while (my $genome= <FH>) { for (0..length($genome)-2) { my $dimer = substr($genome, $_, 2); $dimers[$idx]{$dimer}++ ; #seq wise length of dimer } $idx++; }
|
---|