Skip to main content
se añadieron 59 caracteres en el cuerpo
Origen Enlace

Escribí el siguiente código para calcular la relación entre transiciones y transversiones, es decir, la división entre el total de transiciones y el total de transversiones que ocurren en dos cadenas (cuando se cambia una letra por otra, según el caso).

transition = {'A':'G', 'G':'A', 'T':'C', 'C':'T'}

transversion = {'A':('T','C'), 'G':('T','C'), 'T':('A','G'), 'C':('A','G')}

def rel_trans(s1, s2):
    x = 0
    y = 0
   for i,j in zip(s1, s2):
       if transition[i] == j:
          x += 1 
       elif j in transversion[i]:
          y += 1
   return (x/y)

s1 = 'GCAACGCACAACGAAAACCCTTAGGGACTGGATTATTTCGTGATCGTTGTAGTTATTGGAAGTACGGGCATCAACCCAGTT'
s2 = 'TTATCTGACAAAGAAAGCCGTCAACGGCTGGATAATTTCGCGATCGTGCTGGTTACTGGCGGTACGAGTGTTCCTTTGGGT'

print(rel_trans(s1, s2))

El problema es que este código solo funciona en cadenas de una línea, pero cuando intento ejecutarlo en una de varias líneas, arroja un error. Es decir, quiero que me funcione en lo siguiente:

s1 = """TTTTGTCATTTTGCCTCGTCGTCAGGAATATTGGGACGAATTGGACTCTACT
TCTGTGGGCTGAGCCGGTTCACCCATTGACGAGTGCTAGTCACAGATGAATTACAATGAA
TTGATTTATG"""

s2 = """CCTTGTAAGTACGCGTAGCCTCCAGGAACTTTTAGGCGAATTGTACGTCACC
TCCGGGAGCCGGGCCTCCGGGAGCCGGGCCGGTTGATCCTTTGACAAGTCCTCACAGATGTGTCACAACCAA
TTGATCTG"""

¿Podría alguien decirme qué debo modificar? parece que zip es el problema.

Escribí el siguiente código para calcular la relación entre transiciones y transversiones, es decir, la división entre el total de transiciones y el total de transversiones que ocurren en dos cadenas (cuando se cambia una letra por otra, según el caso).

transition = {'A':'G', 'G':'A', 'T':'C', 'C':'T'}

transversion = {'A':('T','C'), 'G':('T','C'), 'T':('A','G'), 'C':('A','G')}

def rel_trans(s1, s2):
    x = 0
    y = 0
   for i,j in zip(s1, s2):
       if transition[i] == j:
          x += 1 
       elif j in transversion[i]:
          y += 1
   return (x/y)

s1 = 'GCAACGCACAACGAAAACCCTTAGGGACTGGATTATTTCGTGATCGTTGTAGTTATTGGAAGTACGGGCATCAACCCAGTT'
s2 = 'TTATCTGACAAAGAAAGCCGTCAACGGCTGGATAATTTCGCGATCGTGCTGGTTACTGGCGGTACGAGTGTTCCTTTGGGT'

print(rel_trans(s1, s2))

El problema es que este código solo funciona en cadenas de una línea, pero cuando intento ejecutarlo en una de varias líneas, arroja un error. Es decir, quiero que me funcione en lo siguiente:

s1 = """TTTTGTCATTTTGCCTCGTCGTCAGGAATATTGGGACGAATTGGACTCTACT
TCTGTGGGCTGAGCCGGTTCACCCATTGACGAGTGCTAGTCACAGATGAATTACAATGAA
TTGATTTATG"""

s2 = """CCTTGTAAGTACGCGTAGCCTCCAGGAACTTTTAGGCGAATTGTACGTCACC
TCCGGGAGCCGGGCC

¿Podría alguien decirme qué debo modificar? parece que zip es el problema.

Escribí el siguiente código para calcular la relación entre transiciones y transversiones, es decir, la división entre el total de transiciones y el total de transversiones que ocurren en dos cadenas (cuando se cambia una letra por otra, según el caso).

transition = {'A':'G', 'G':'A', 'T':'C', 'C':'T'}

transversion = {'A':('T','C'), 'G':('T','C'), 'T':('A','G'), 'C':('A','G')}

def rel_trans(s1, s2):
    x = 0
    y = 0
   for i,j in zip(s1, s2):
       if transition[i] == j:
          x += 1 
       elif j in transversion[i]:
          y += 1
   return (x/y)

s1 = 'GCAACGCACAACGAAAACCCTTAGGGACTGGATTATTTCGTGATCGTTGTAGTTATTGGAAGTACGGGCATCAACCCAGTT'
s2 = 'TTATCTGACAAAGAAAGCCGTCAACGGCTGGATAATTTCGCGATCGTGCTGGTTACTGGCGGTACGAGTGTTCCTTTGGGT'

print(rel_trans(s1, s2))

El problema es que este código solo funciona en cadenas de una línea, pero cuando intento ejecutarlo en una de varias líneas, arroja un error. Es decir, quiero que me funcione en lo siguiente:

s1 = """TTTTGTCATTTTGCCTCGTCGTCAGGAATATTGGGACGAATTGGACTCTACT
TCTGTGGGCTGAGCCGGTTCACCCATTGACGAGTGCTAGTCACAGATGAATTACAATGAA
TTGATTTATG"""

s2 = """CCTTGTAAGTACGCGTAGCCTCCAGGAACTTTTAGGCGAATTGTACGTCACC
TCCGGGAGCCGGGCCGGTTGATCCTTTGACAAGTCCTCACAGATGTGTCACAACCAA
TTGATCTG"""

¿Podría alguien decirme qué debo modificar? parece que zip es el problema.

se añadieron 12 caracteres en el cuerpo
Origen Enlace
Roger Torné
  • 7.6k
  • 13
  • 20
  • 33

Escribí el siguiente código para calcular la relación entre transiciones y transversiones, es decir, la división entre el total de transiciones y el total de transversiones que ocurren en dos cadenas (cuando se cambia una letra por otra, según el caso).

transition = {'A':'G', 'G':'A', 'T':'C', 'C':'T'}

transversion = {'A':('T','C'), 'G':('T','C'), 'T':('A','G'), 'C':('A','G')}

transition = {'A':'G', 'G':'A', 'T':'C', 'C':'T'}

transversion = {'A':('T','C'), 'G':('T','C'), 'T':('A','G'), 'C':('A','G')}

def rel_trans(s1, s2):
    x = 0
    y = 0
   for i,j in zip(s1, s2):
       if transition[i] == j:
          x += 1 
       elif j in transversion[i]:
          y += 1
   return (x/y)

s1 = 'GCAACGCACAACGAAAACCCTTAGGGACTGGATTATTTCGTGATCGTTGTAGTTATTGGAAGTACGGGCATCAACCCAGTT'
s2 = 'TTATCTGACAAAGAAAGCCGTCAACGGCTGGATAATTTCGCGATCGTGCTGGTTACTGGCGGTACGAGTGTTCCTTTGGGT'

print(rel_trans(s1, s2))

El problema es que este código solo funciona en cadenas de una línea, pero cuando intento ejecutarlo en una de varias líneas, arroja un error. Es decir, quiero que me funcione en lo siguiente:

s1 = """TTTTGTCATTTTGCCTCGTCGTCAGGAATATTGGGACGAATTGGACTCTACT
TCTGTGGGCTGAGCCGGTTCACCCATTGACGAGTGCTAGTCACAGATGAATTACAATGAA
TTGATTTATG"""

s2 = """CCTTGTAAGTACGCGTAGCCTCCAGGAACTTTTAGGCGAATTGTACGTCACC
TCCGGGAGCCGGGCC

¿Podría alguien decirme qué debo modificar? parece que zip es el problema.

Escribí el siguiente código para calcular la relación entre transiciones y transversiones, es decir, la división entre el total de transiciones y el total de transversiones que ocurren en dos cadenas (cuando se cambia una letra por otra, según el caso).

transition = {'A':'G', 'G':'A', 'T':'C', 'C':'T'}

transversion = {'A':('T','C'), 'G':('T','C'), 'T':('A','G'), 'C':('A','G')}

def rel_trans(s1, s2):
    x = 0
    y = 0
   for i,j in zip(s1, s2):
       if transition[i] == j:
          x += 1 
       elif j in transversion[i]:
          y += 1
   return (x/y)

s1 = 'GCAACGCACAACGAAAACCCTTAGGGACTGGATTATTTCGTGATCGTTGTAGTTATTGGAAGTACGGGCATCAACCCAGTT'
s2 = 'TTATCTGACAAAGAAAGCCGTCAACGGCTGGATAATTTCGCGATCGTGCTGGTTACTGGCGGTACGAGTGTTCCTTTGGGT'

print(rel_trans(s1, s2))

El problema es que este código solo funciona en cadenas de una línea, pero cuando intento ejecutarlo en una de varias líneas, arroja un error. Es decir, quiero que me funcione en lo siguiente:

s1 = """TTTTGTCATTTTGCCTCGTCGTCAGGAATATTGGGACGAATTGGACTCTACT
TCTGTGGGCTGAGCCGGTTCACCCATTGACGAGTGCTAGTCACAGATGAATTACAATGAA
TTGATTTATG"""

s2 = """CCTTGTAAGTACGCGTAGCCTCCAGGAACTTTTAGGCGAATTGTACGTCACC
TCCGGGAGCCGGGCC

¿Podría alguien decirme qué debo modificar? parece que zip es el problema.

Escribí el siguiente código para calcular la relación entre transiciones y transversiones, es decir, la división entre el total de transiciones y el total de transversiones que ocurren en dos cadenas (cuando se cambia una letra por otra, según el caso).

transition = {'A':'G', 'G':'A', 'T':'C', 'C':'T'}

transversion = {'A':('T','C'), 'G':('T','C'), 'T':('A','G'), 'C':('A','G')}

def rel_trans(s1, s2):
    x = 0
    y = 0
   for i,j in zip(s1, s2):
       if transition[i] == j:
          x += 1 
       elif j in transversion[i]:
          y += 1
   return (x/y)

s1 = 'GCAACGCACAACGAAAACCCTTAGGGACTGGATTATTTCGTGATCGTTGTAGTTATTGGAAGTACGGGCATCAACCCAGTT'
s2 = 'TTATCTGACAAAGAAAGCCGTCAACGGCTGGATAATTTCGCGATCGTGCTGGTTACTGGCGGTACGAGTGTTCCTTTGGGT'

print(rel_trans(s1, s2))

El problema es que este código solo funciona en cadenas de una línea, pero cuando intento ejecutarlo en una de varias líneas, arroja un error. Es decir, quiero que me funcione en lo siguiente:

s1 = """TTTTGTCATTTTGCCTCGTCGTCAGGAATATTGGGACGAATTGGACTCTACT
TCTGTGGGCTGAGCCGGTTCACCCATTGACGAGTGCTAGTCACAGATGAATTACAATGAA
TTGATTTATG"""

s2 = """CCTTGTAAGTACGCGTAGCCTCCAGGAACTTTTAGGCGAATTGTACGTCACC
TCCGGGAGCCGGGCC

¿Podría alguien decirme qué debo modificar? parece que zip es el problema.

Origen Enlace

Cálculo entre transición y transversión para dos secuencias con python

Escribí el siguiente código para calcular la relación entre transiciones y transversiones, es decir, la división entre el total de transiciones y el total de transversiones que ocurren en dos cadenas (cuando se cambia una letra por otra, según el caso).

transition = {'A':'G', 'G':'A', 'T':'C', 'C':'T'}

transversion = {'A':('T','C'), 'G':('T','C'), 'T':('A','G'), 'C':('A','G')}

def rel_trans(s1, s2):
    x = 0
    y = 0
   for i,j in zip(s1, s2):
       if transition[i] == j:
          x += 1 
       elif j in transversion[i]:
          y += 1
   return (x/y)

s1 = 'GCAACGCACAACGAAAACCCTTAGGGACTGGATTATTTCGTGATCGTTGTAGTTATTGGAAGTACGGGCATCAACCCAGTT'
s2 = 'TTATCTGACAAAGAAAGCCGTCAACGGCTGGATAATTTCGCGATCGTGCTGGTTACTGGCGGTACGAGTGTTCCTTTGGGT'

print(rel_trans(s1, s2))

El problema es que este código solo funciona en cadenas de una línea, pero cuando intento ejecutarlo en una de varias líneas, arroja un error. Es decir, quiero que me funcione en lo siguiente:

s1 = """TTTTGTCATTTTGCCTCGTCGTCAGGAATATTGGGACGAATTGGACTCTACT
TCTGTGGGCTGAGCCGGTTCACCCATTGACGAGTGCTAGTCACAGATGAATTACAATGAA
TTGATTTATG"""

s2 = """CCTTGTAAGTACGCGTAGCCTCCAGGAACTTTTAGGCGAATTGTACGTCACC
TCCGGGAGCCGGGCC

¿Podría alguien decirme qué debo modificar? parece que zip es el problema.