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))