def test_unambig_dna(self):
unambig_dna_seq = Seq("AGCTTCGAAATGCT")
seqs = list(disambiguate_iupac_dna(unambig_dna_seq))
assert len(seqs) == 1
assert seqs[0] == unambig_dna_seq
def test_ambig_dna(self):
ambig_dna_seq = Seq("NGCTT")
# map to str here bc of annoying biopython warning when hashing a Seq
seqs = set(map(str, disambiguate_iupac_dna(ambig_dna_seq)))
assert len(seqs) == 4
assert seqs == {"AGCTT", "CGCTT", "GGCTT", "TGCTT"}
def test_ambig_dna(self):
ambig_dna_seq = Seq("NGCTT", ambiguous_dna)
seqs = set(disambiguate_iupac_dna(ambig_dna_seq))
assert len(seqs) == 4
assert seqs == {"AGCTT", "CGCTT", "GGCTT", "TGCTT"}