def test_revcompl_dna():
# validate test data
test_data_dna.inspect()
mapper = CharsMapper.create_mapper_compl_dna()
revcompl = mapper.transcode(test_data_dna['+']).decode('ascii')[::-1]
assert revcompl == test_data_dna['-']
def render(cls, args, outfile):
compl = CharsMapper.create_mapper_compl_dna()
for fn in args.filenames:
for rec in fasta.read(fn):
rec.cmt += ".RC"
rec.seq = compl.transcode(rec.seq)[::-1]
fasta.write(outfile, rec)
