in reply to Re^2: Create the reverse complement DNA sequence without pattern matching and reverse built-in function?
in thread Create the reverse complement DNA sequence without pattern matching and reverse built-in function?
I tried to alert the OP to this issue in my post:
Sorry. I missed that.
your tr/// method missed reversing the string afterwards, but that won't change the fact that the tr/// approach is the way to go.
Right on both counts. And it make's more of a difference than I was expecting, but still 80x faster than the next best:
tr => q[ ( my $revcmp = reverse $seq ) =~ tr[ACGTacgt][tgcaTGCA]; ], C:\test>junk60 Rate hdb atcroft davido tr hdb 9827/s -- -36% -71% -100% atcroft 15292/s 56% -- -55% -99% davido 33822/s 244% 121% -- -99% tr 2640400/s 26770% 17167% 7707% --
|
|---|