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)