Beefy Boxes and Bandwidth Generously Provided by pair Networks
Perl: the Markov chain saw
 
PerlMonks  

Re^5: mismatching characters in dna sequence

by choroba (Cardinal)
on Dec 30, 2011 at 02:36 UTC ( [id://945588]=note: print w/replies, xml ) Need Help??


in reply to Re^4: mismatching characters in dna sequence
in thread mismatching characters in dna sequence

You can get the information back from the strings if needed:
use warnings; use strict; my $source="ATTCCGGG"; print "source: $source\n"; for my $str ( "ATTGCGGG", "ATACCGGC", "ACTGCGGT" , "ATCGGGGG" ) { print "string: $str\n"; my @res = unpack "c*", $source ^ $str; for my $i (grep $res[$_], 0 .. $#res) { print "$i: ", substr($source, $i, 1), "->", substr($str, $i, 1), "\n"; } }
I am not sure how fast this is compared to the PDL solution.

Log In?
Username:
Password:

What's my password?
Create A New User
Domain Nodelet?
Node Status?
node history
Node Type: note [id://945588]
help
Chatterbox?
and the web crawler heard nothing...

How do I use this?Last hourOther CB clients
Other Users?
Others surveying the Monastery: (4)
As of 2024-03-28 20:00 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    No recent polls found