import collections
seq = ''
count = 0

with open('sequence.fasta', 'r') as inp:
  line = inp.readline()
  while line:
    seq += line.strip()
    line = inp.readline()
  for i in range(len(seq)):
    if seq[i:i+2] == 'TA':
      count += 1
  print(count)

nuc_counter = collections.Counter(seq)

cnt_A = nuc_counter['A']/len(seq)
cnt_T = nuc_counter['T']/len(seq)

exp_num = round((cnt_T * cnt_A) * len(seq))   
print(exp_num)
print(count/exp_num)