2007-11-16 04:30:29 +03:00
|
|
|
#!/usr/bin/ruby
|
|
|
|
|
# The Great Computer Language Shootout
|
|
|
|
|
# http://shootout.alioth.debian.org/
|
|
|
|
|
#
|
|
|
|
|
# Contributed by Peter Bjarke Olsen
|
|
|
|
|
# Modified by Doug King
|
|
|
|
|
|
|
|
|
|
seq=Array.new
|
|
|
|
|
|
|
|
|
|
def revcomp(seq)
|
|
|
|
|
seq.reverse!.tr!('wsatugcyrkmbdhvnATUGCYRKMBDHVN','WSTAACGRYMKVHDBNTAACGRYMKVHDBN')
|
|
|
|
|
stringlen=seq.length
|
|
|
|
|
0.step(stringlen-1,60) {|x| print seq.slice(x,60) , "\n"}
|
|
|
|
|
end
|
|
|
|
|
|
|
|
|
|
input = open(File.join(File.dirname($0), 'fasta.output.2500000'), 'rb')
|
|
|
|
|
|
|
|
|
|
while input.gets
|
|
|
|
|
if $_ =~ />/
|
|
|
|
|
if seq.length != 0
|
|
|
|
|
revcomp(seq.join)
|
|
|
|
|
seq=Array.new
|
|
|
|
|
end
|
|
|
|
|
puts $_
|
|
|
|
|
else
|
|
|
|
|
$_.sub(/\n/,'')
|
|
|
|
|
seq.push $_
|
|
|
|
|
end
|
|
|
|
|
end
|
|
|
|
|
revcomp(seq.join)
|
