#! /usr/bin/env python

# driverAltschulEriksonForHPmodel.py
# P.Clote

from altschulEriksonDinuclShuffle import dinuclShuffle



import sys,os

def HP2AU(seq):
  s = ""
  for ch in seq:
    if ch == 'H':
      s += 'A'
    elif ch == 'P':
      s += 'U'
  return s

def AU2HP(seq):
  s = ""
  for ch in seq:
    if ch == 'A':
      s += 'H'
    elif ch == 'U':
      s += 'P'
  return s


def main(filename):
  file  = open(filename)
  line  = file.readline()
  while line:
    FASTA = line.strip()
    assert( FASTA[0]=='>')
    rna  = HP2AU(file.readline().strip())
    n    = len(rna)
    for i in range(10):
      seq = AU2HP(dinuclShuffle(rna))
      print "%s_random%d of length %d" % (FASTA,i,n)
      print "%s" % seq
    line = file.readline()
  file.close()

if __name__ == '__main__':
  if len(sys.argv) < 2:
    print "Usage: %s filename" % sys.argv[0]
    sys.exit(1)
  filename = sys.argv[1]
  main(filename)