44 lines
1.1 KiB
Factor
44 lines
1.1 KiB
Factor
|
USING: io io.files io.streams.duplex kernel sequences
|
||
|
|
sequences.private strings vectors words memoize splitting
|
||
|
|
hints ;
|
||
|
|
IN: benchmark.reverse-complement
|
||
|
|
|
||
|
|
MEMO: trans-map ( -- str )
|
||
|
|
256 >string
|
||
|
|
"TGCAAKYRMBDHV" "ACGTUMRYKVHDB"
|
||
|
|
[ pick set-nth ] 2each ;
|
||
|
|
|
||
|
|
: do-trans-map ( str -- )
|
||
|
|
[ ch>upper trans-map nth ] change-each ;
|
||
|
|
|
||
|
|
HINTS: do-trans-map string ;
|
||
|
|
|
||
|
|
: translate-seq ( seq -- str )
|
||
|
|
concat dup reverse-here dup do-trans-map ;
|
||
|
|
|
||
|
|
: show-seq ( seq -- )
|
||
|
|
translate-seq 60 <groups> [ print ] each ;
|
||
|
|
|
||
|
|
: do-line ( seq line -- seq )
|
||
|
|
dup first ">;" memq? [
|
||
|
|
over show-seq print dup delete-all
|
||
|
|
] [
|
||
|
|
over push
|
||
|
|
] if ;
|
||
|
|
|
||
|
|
: (reverse-complement) ( seq -- )
|
||
|
|
readln [ do-line (reverse-complement) ] [ show-seq ] if* ;
|
||
|
|
|
||
|
|
: reverse-complement ( infile outfile -- )
|
||
|
|
<file-writer> >r <file-reader> r> <duplex-stream> [
|
||
|
|
500000 <vector> (reverse-complement)
|
||
|
|
] with-stream ;
|
||
|
|
|
||
|
|
: reverse-complement-main ( -- )
|
||
|
|
"reverse-complement-in.txt"
|
||
|
|
"reverse-complement-out.txt"
|
||
|
|
[ home swap path+ ] 2apply
|
||
|
|
reverse-complement ;
|
||
|
|
|
||
|
|
MAIN: reverse-complement-main
|