# Path: dark/reads.py
# class DNARead(_NucleotideRead):
#     """
#     Hold information and methods to work with DNA reads.
#     """
#     ALPHABET = set('ATCG')
# 
#     COMPLEMENT_TABLE = _makeComplementTable(ambiguous_dna_complement)


import sys
from numpy.random import choice
from dark.fasta import FastaReads
from dark.reads import DNARead
#!/usr/bin/env python

"""
Read DNA FASTA from stdin and print FASTA to stdout, with the sequences
randomized. The read ids and lengths are preserved.

Note: This produces DNA sequences. If you have AA reads and you need this
      functionality, we can add it.
"""

from __future__ import print_function




if __name__ == '__main__':
    for read in FastaReads(sys.stdin):
        seq = ''.join(choice(['A', 'C', 'G', 'T'], len(read), replace=True))

# Path: dark/reads.py
# class DNARead(_NucleotideRead):
#     """
#     Hold information and methods to work with DNA reads.
#     """
#     ALPHABET = set('ATCG')
# 
#     COMPLEMENT_TABLE = _makeComplementTable(ambiguous_dna_complement)


# Path: bin/randomize-fasta.py
import sys
from numpy.random import choice
from dark.fasta import FastaReads
from dark.reads import DNARead

#!/usr/bin/env python

"""
Read DNA FASTA from stdin and print FASTA to stdout, with the sequences
randomized. The read ids and lengths are preserved.

Note: This produces DNA sequences. If you have AA reads and you need this
      functionality, we can add it.
"""

from __future__ import print_function




if __name__ == '__main__':
    for read in FastaReads(sys.stdin):
        seq = ''.join(choice(['A', 'C', 'G', 'T'], len(read), replace=True))

        print(DNARead(read.id, seq).toString('fasta'))