def test_iupac_consensus(self): filename = "GFF/multi_1.fna" format = "fasta" alignment = AlignIO.read(filename, format, alphabet=unambiguous_dna) summary = SummaryInfo(alignment) c = summary.iupac_consensus(ambiguous="N") self.assertEqual(str(c), 'TAWSRWYSNBATCG') self.assertNotEqual(c.alphabet, unambiguous_dna) self.assertTrue(isinstance(c.alphabet, DNAAlphabet)) c = summary.gap_iupac_consensus(ambiguous="N") self.assertEqual(str(c), 'TAWSRWYSNBATCG') self.assertNotEqual(c.alphabet, unambiguous_dna) self.assertTrue(isinstance(c.alphabet, DNAAlphabet))