comment on

I have a large set (100000+) of short DNA reads 20 characters long. I need to compares all reads against each other and pull out those that vary by just 1 position. Heres the script I came up with.

$| = 1;
my $compare_count = 0;
for (my $i = 0; $i < @kmers; $i++ ) {
    for (my $j = $i + 1; $j < @kmers; $j++ ) {
        print "\rComparing sequence $i to $j";
        my @result = PCCompare::dissimilarity($kmers[$i], $kmers[$j], 
+1);
        if ($result[0] == 1) {
            print "\rMatch found: $kmers[$i], $kmers[$j]\n";
            push @variant_kmers, ($kmers[$i], $kmers[$j]);
        }
        $compare_count++;
    }
}
print "\rFinished: $compare_count comparisions made.\n";
[download]

The problem is that this loop runs very very slow. It takes on the orders of days to process 100000 sequences. Is there a way to make the process faster?

In reply to Comparing a large set of DNA sequences by shamshersingh

Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!

Titles consisting of a single word are discouraged, and in most cases are disallowed outright.

Read Where should I post X? if you're not absolutely sure you're posting in the right place.

Please read these before you post! —

Posts may use any of the Perl Monks Approved HTML tags:

a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr

You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)

	For:		Use:
	&		`&`
	<		`<`
	>		`>`
	[		`[`
	]		`]`

Link using PerlMonks shortcuts! What shortcuts can I use for linking?

See Writeup Formatting Tips and other pages linked from there for more info.


Do you know where your variables are?
	PerlMonks