def testReverseComplementAA(self):
"""
The reverseComplement function must raise a C{ValueError} when called
on an amino acid sequence.
"""
read = Read('id', 'atcg', type='aa')
error = 'Cannot reverse complement an amino acid sequence'
with self.assertRaisesRegexp(ValueError, error):
read.reverseComplement()
def testReverseComplementAmbiguousRNA(self):
"""
The reverseComplement function must work for RNA that includes
ambiguous bases.
"""
read = Read('id', 'aucgmrwsykvhxn', type='rna')
self.assertEqual('NXDBMRSWYKCGAU', read.reverseComplement().sequence)
def testReverseComplementRNA(self):
"""
The reverseComplement function must work for RNA
"""
read = Read('id', 'aucg', type='rna')
self.assertEqual('CGAU', read.reverseComplement().sequence)
def testReverseComplementDNA(self):
"""
The reverseComplement function must work for DNA
"""
read = Read('id', 'atcg', quality='!@#$', type='dna')
self.assertEqual('CGAT', read.reverseComplement().sequence)
def testReverseComplementReversesQuality(self):
"""
The reverseComplement function must return a reversed quality string.
"""
read = Read('id', 'atcg', quality='!@#$')
self.assertEqual('$#@!', read.reverseComplement().quality)
