read-sequences-v1.py
485 Bytes
if __name__ == "__main__":
# Abrir archivo con secuencias
listLines = []
length_seq = 0
count_seq = 0
with open("hga-sequences-1000.txt", mode="r", encoding="utf8") as iFile:
listLines = iFile.readlines()
for line in listLines:
if not line.startswith("seqname"):
listElements = line.split("\t")
length_seq += len(listElements[4])
count_seq += 1
print("Longitud promedio: {}".format(length_seq/count_seq))