def test_complement_ambiguous_rna_values(self):
for ambig_char, values in sorted(ambiguous_rna_values.items()):
# Will default to DNA if neither T nor U found...
compl_values = str(Seq.Seq(values).complement().transcribe())
ambig_values = ambiguous_rna_values[
ambiguous_rna_complement[ambig_char]]
self.assertEqual(set(compl_values), set(ambig_values))
def test_complement_ambiguous_rna_values(self):
for ambig_char, values in sorted(ambiguous_rna_values.items()):
compl_values = str(
Seq.Seq(values, alphabet=IUPAC.ambiguous_rna).complement())
self.assertEqual(
set(compl_values),
set(ambiguous_rna_values[
ambiguous_rna_complement[ambig_char]]))
def test_complement_ambiguous_rna_values(self):
for ambig_char, values in sorted(ambiguous_rna_values.items()):
compl_values = str(
# Will default to DNA if neither T or U found...
# Turn black code style off
# fmt: off
Seq.Seq(values).complement().transcribe()
# fmt: on
)
ambig_values = ambiguous_rna_values[
ambiguous_rna_complement[ambig_char]]
self.assertEqual(set(compl_values), set(ambig_values))
print("")
print("DNA Ambiguity mapping: %s" % sorted_dict(ambiguous_dna_values))
print("DNA Complement mapping: %s" % sorted_dict(ambiguous_dna_complement))
for ambig_char, values in sorted(ambiguous_dna_values.items()):
compl_values = complement(values)
print("%s={%s} --> {%s}=%s" % (ambig_char, values, compl_values,
ambiguous_dna_complement[ambig_char]))
assert set(compl_values) == set(
ambiguous_dna_values[ambiguous_dna_complement[ambig_char]])
print("")
print("RNA Ambiguity mapping: %s" % sorted_dict(ambiguous_rna_values))
print("RNA Complement mapping: %s" % sorted_dict(ambiguous_rna_complement))
for ambig_char, values in sorted(ambiguous_rna_values.items()):
compl_values = complement(values).replace(
"T", "U") # need to help as no alphabet
print("%s={%s} --> {%s}=%s" % (ambig_char, values, compl_values,
ambiguous_rna_complement[ambig_char]))
assert set(compl_values) == set(
ambiguous_rna_values[ambiguous_rna_complement[ambig_char]])
print("")
print("Reverse complements:")
for sequence in [
Seq.Seq("".join(sorted(ambiguous_rna_values))),
Seq.Seq("".join(sorted(ambiguous_dna_values))),
Seq.Seq("".join(sorted(ambiguous_rna_values)), Alphabet.generic_rna),
Seq.Seq("".join(sorted(ambiguous_dna_values)), Alphabet.generic_dna),
Seq.Seq("".join(sorted(ambiguous_rna_values)).replace("X", ""),
return "{%s}" % ", ".join("%s: %s" % (repr(k), repr(v))
for k, v in sorted(d.items()))
print("")
print("DNA Ambiguity mapping: %s" % sorted_dict(ambiguous_dna_values))
print("DNA Complement mapping: %s" % sorted_dict(ambiguous_dna_complement))
for ambig_char, values in sorted(ambiguous_dna_values.items()):
compl_values = complement(values)
print("%s={%s} --> {%s}=%s" %
(ambig_char, values, compl_values, ambiguous_dna_complement[ambig_char]))
assert set(compl_values) == set(ambiguous_dna_values[ambiguous_dna_complement[ambig_char]])
print("")
print("RNA Ambiguity mapping: %s" % sorted_dict(ambiguous_rna_values))
print("RNA Complement mapping: %s" % sorted_dict(ambiguous_rna_complement))
for ambig_char, values in sorted(ambiguous_rna_values.items()):
compl_values = complement(values).replace("T", "U") # need to help as no alphabet
print("%s={%s} --> {%s}=%s" %
(ambig_char, values, compl_values, ambiguous_rna_complement[ambig_char]))
assert set(compl_values) == set(ambiguous_rna_values[ambiguous_rna_complement[ambig_char]])
print("")
print("Reverse complements:")
for sequence in [Seq.Seq("".join(sorted(ambiguous_rna_values))),
Seq.Seq("".join(sorted(ambiguous_dna_values))),
Seq.Seq("".join(sorted(ambiguous_rna_values)), Alphabet.generic_rna),
Seq.Seq("".join(sorted(ambiguous_dna_values)), Alphabet.generic_dna),
Seq.Seq("".join(sorted(ambiguous_rna_values)).replace("X", ""), IUPAC.IUPACAmbiguousRNA()),
Seq.Seq("".join(sorted(ambiguous_dna_values)).replace("X", ""), IUPAC.IUPACAmbiguousDNA()),
Seq.Seq("AWGAARCKG")]: # Note no U or T
print("%s -> %s"
def test_complement_ambiguous_rna_values(self):
for ambig_char, values in sorted(ambiguous_rna_values.items()):
compl_values = str(Seq.Seq(values, alphabet=IUPAC.ambiguous_rna).complement())
self.assertEqual(set(compl_values),
set(ambiguous_rna_values[ambiguous_rna_complement[ambig_char]]))