def testGC(self):
"""Verify computation of GC."""
binStats = BinStatistics(threads=1)
seqs = {'S1': 'ACgt', 'S2': 'GGgg', 'S3': 'TTtt', 'S4': 'NNNN'}
seqStats = defaultdict(dict)
meanGC, _ = binStats.calculateGC(seqs, seqStats)
self.assertAlmostEqual(seqStats['S1']['GC'], 0.5)
self.assertAlmostEqual(seqStats['S2']['GC'], 1.0)
self.assertAlmostEqual(seqStats['S3']['GC'], 0.0)
self.assertAlmostEqual(seqStats['S4']['GC'], 0.0)
self.assertAlmostEqual(meanGC, 6.0 / 12.0)
def testGC(self):
"""Verify computation of GC."""
binStats = BinStatistics(threads = 1)
seqs = {'S1':'ACgt', 'S2':'GGgg', 'S3':'TTtt', 'S4':'NNNN'}
seqStats = defaultdict(dict)
meanGC, _ = binStats.calculateGC(seqs, seqStats)
self.assertAlmostEqual(seqStats['S1']['GC'], 0.5)
self.assertAlmostEqual(seqStats['S2']['GC'], 1.0)
self.assertAlmostEqual(seqStats['S3']['GC'], 0.0)
self.assertAlmostEqual(seqStats['S4']['GC'], 0.0)
self.assertAlmostEqual(meanGC, 6.0/12.0)
