Estoy trabajando con secuencias de ADN y necesito hacer una clase que tenga un método para calcular el ratio de guanina-citosina, otro para hacer el complementario de una cadena de ADN y otro más que haga el complementario y que además lo invierta, además de pruebas unitarias. El primer método consigo que me funcione, el problema viene en el segundo. Este es código que he hecho:
class SequenceBase:
def __init__ (self, secuencia_nucleotidos):
if len (secuencia_nucleotidos) == None:
raise Exception
else:
self.secuencia_nucleotidos = secuencia_nucleotidos
def __str__ (self):
return self.secuencia_nucleotidos.upper()
def cadenas_iguales (self, otra_secuencia_nucleotidos):
return (self.secuencia_nucleotidos == otra_secuencia_nucleotidos.secuencia_nucleotidos)
def gc_content (self): #Cálculo ratio GC --> ((G+C)/(A+T+G+C))*100
contador_guanina = self.secuencia_nucleotidos.count("G")
contador_citosina = self.secuencia_nucleotidos.count ("C")
contador_adenina = self.secuencia_nucleotidos.count ("A")
contador_timina = self.secuencia_nucleotidos.count ("T")
proporcion_guanina_citosina = float (((contador_guanina + contador_citosina)/(contador_adenina + contador_citosina + contador_timina + contador_guanina))*100)
return round (proporcion_guanina_citosina, 2)
def complement (self):
complementario = self.secuencia_nucleotidos.maketrans ("ATCG","TAGC")
ADN_complementario = self.secuencia_nucleotidos.translate(complementario)
return ADN_complementario
def reverse_complement (self):
ADN_complementario = self.complement ()
cadena_invertida_ADN = reversed(ADN_complementario)
return SequenceBase (cadena_invertida_ADN)
Cuando quiero hacer la prueba unitaria para complement me sale el error AtributeError: 'str' object has no attribute 'complement'. Aquí están mis pruebas unitarias:
def test_complement (self):
self.secuencia_ADN = "AAAACCCACGTTTGCCCG"
self.secuencia_ADN_complementaria = "TTTTGGGTGCAAACGGGC"
self.secuencia_ADN_no_complementaria = "ACCAGT"
self.assertEqual (SequenceBase.complement (self.secuencia_ADN), self.secuencia_ADN_complementaria)
self.assertEqual (SequenceBase.complement (self.secuencia_ADN_complementaria), self.secuencia_ADN)
self.assertNotEqual (SequenceBase.complement (self.secuencia_ADN), self.secuencia_ADN_no_complementaria)
self.assertNotEqual (SequenceBase.complement (self.secuencia_ADN_complementaria), self.secuencia_ADN_no_complementaria)
def test_reverse_complement (self):
self.secuencia_ADN_complementaria = "AAACCAGTTGTGTACACAGT"
self.secuencia_ADN_inversa = "TGACACATGTGTTGACCAAA"
self.assertEqual (SequenceBase.reverse_complement (self.secuencia_ADN_complementaria),self.secuencia_ADN_inversa)
self.assertEqual (SequenceBase.reverse_complement (self.secuencia_ADN_inversa),self.secuencia_ADN_complementaria)
Les agradecería si me pudiesen ayudar