Re: Comparing a large set of DNA sequences

in reply to Comparing a large set of DNA sequences

This will perform the same job, but run significantly faster (Approx 7 hours).

$| = 1;
my $compare_count = 0;
for (my $i = 0; $i < @kmers; $i++ ) {
    for (my $j = $i + 1; $j < @kmers; $j++ ) {
        my $sims = ( $kmers[ $i ] ^ $kmers[ $j ] ) =~ tr[\0][];
        if(  $sims == 19 ) {
            print "Match found: $kmers[ $i ], $kmers[ $j ]";
            push @varient_kmers, ($kmers[$i], $kmers[$j]);
        }
        ++$compare_count;
    }
}
print "\rFinished: $compare_count comparisions made.\n";
[download]

With the rise and rise of 'Social' network sites: 'Computers are making people easier to use everyday'

Examine what is said, not who speaks -- Silence betokens consent -- Love the truth but pardon error.

"Science is about questioning the status quo. Questioning authority".

In the absence of evidence, opinion is indistinguishable from prejudice.

In Section Seekers of Perl Wisdom