Puedes usar collections.Counter
para contar las bases de cada línea. Para saber si es ARN u ADN basta con que comprobemos si la clave "U" está o no en el diccionario.
También es importante, sobre todo si vas a trabajar con archivos relativamente grandes, no recorrerlos más de una vez.
Hay múltiples formas de organizar el código, una de ellas puede ser:
import collections
import os
def registros(fichero):
registros = collections.OrderedDict()
contador = 0
for linea in fichero:
if linea.startswith('>'):
contador += 1
registros[contador] = collections.Counter()
else:
registros[contador].update(linea)
return registros
def parse_fasta(path):
try:
with open(path, 'r') as f:
if os.path.isfile(path):
print('Se encontró Fichero Fasta:', fichero)
regs = registros(f)
print("El fichero {} contiene {} registros con el siguiente contenido de bases:".format(fichero, len(regs)))
for reg, bases in regs.items():
print(" Registro {}:".format(reg))
if "U" in bases:
print(" Es un registro de tipo ARN:")
print(" Adenina: {}\n Uracilo: {}\n Citosina: {}\n Guanina: {}\n".format(bases["A"],
bases["U"],
bases["C"],
bases["G"]))
else:
print(" Es un registro de tipo ADN:")
print(" Adenina: {}\n Timina: {}\n Citosina: {}\n Guanina: {}\n".format(bases["A"],
bases["T"],
bases["C"],
bases["G"]))
else:
print("La ruta no se corresponde con un fichero.")
except FileNotFoundError:
print('El fichero introducido no se ha encontrado, asegúrese de que se encuentra en ese directorio!')
while True:
fichero = input('Introduzca nombre del fichero FASTA(q para salir):\n')
if fichero == 'q':
break
parse_fasta(fichero)
La salida para tu ejemplo es:
>>> Introduzca nombre del fichero FASTA(q para salir):
2.fasta
Se encontró Fichero Fasta: 2.fasta
El fichero 2.fasta contiene 2 registros con el siguiente contenido de bases:
Registro 1:
Es un registro de tipo ARN:
Adenina: 2
Uracilo: 1
Citosina: 6
Guanina: 3
Registro 2:
Es un registro de tipo ADN:
Adenina: 2
Timina: 2
Citosina: 3
Guanina: 4