reads_information

def reads_readfile():
    reads_file = open("reads7.fa","r")
    reads = reads_file.readlines()
    reads_file.close()
    return reads

def reads_check(c_reads):
    reads_amount,reads_length_all, reads_all = 0, 0, []
    for line in c_reads:
        line.replace("\n", "")
        if line[0].isalpha():
            reads_amount = reads_amount+1
            reads_length_all = reads_length_all+len(line)
            reads_all.append(line)
    reads_average_length = reads_length_all/reads_amount
    return reads_amount, reads_average_length, reads_all


def longest_read(l_reads_all):
    fasta = ">unknown\n"+max(l_reads_all, key=len)
    fasta_file = open("fasta_sequence.txt", "w")
    fasta_file.write(fasta)
    fasta_file.close()


def main():
    reads = reads_readfile()
    reads_amount, reads_average_length, reads_all = reads_check(reads)
    print("Amount of reads: "+str(reads_amount)+"\nAverage length of "\
          "reads: "+str(reads_average_length))
    longest_read(reads_all)

main()

print("\nFinished.")