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)