Python скрипт с использованием ункций Emboss

На Python был написан скрипт, который находит динуклеотид, реальная частота которого наиболее отличается от ожидаемой.
Ожидаемая частота была рассчитана как произведение частот нуклеотидов, составляющих динуклеотид. При написании программа испытывалась на Python 3.6, поэтому рекоммендуется использование этой версии.
Скрипт
Формат ввода: 'python3 emboss.py ***filename***', где ***filename*** - имя файла с последовательностью. Программа выводит динуклеотид, для которого разность реальной и ожидаемой частот была максимальна и саму разность.
Тестовый запуск был произведен на геноме Escherichia coli файл . Наибольшая разница реальной и ожаидаемой частот была найдена для динуклеотида CG и она составила 0.03261496217813238