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))
