Exemple #1
0
def test_map_features2alignment(key, next_hash, alb_resources, hf):
    alignbuddy = alb_resources.get_one(key)
    if alignbuddy.alpha == IUPAC.protein:
        seqbuddy = Sb.SeqBuddy("%s%sMnemiopsis_pep.gb" % (hf.resource_path, os.path.sep))
    else:
        seqbuddy = Sb.SeqBuddy("%s%sMnemiopsis_cds.gb" % (hf.resource_path, os.path.sep))
    tester = Alb.map_features2alignment(seqbuddy, alignbuddy)
    tester.set_format("genbank")
    assert hf.buddy2hash(tester) == next_hash, tester.write("error_files%s%s" % (next_hash, os.path.sep))