#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Mon Oct 26 11:52:56 2020

@author: lewis
"""

from sys import argv
from Bio import SeqIO
from statistics import mean

filename1 = argv[1]
filename2 = argv[2]

rec1 = SeqIO.read(filename1, 'abi')
rec2 = SeqIO.read(filename2, 'abi')

signal1 = rec1.annotations['abif_raw']['S/N%1']
noise1 = rec1.annotations['abif_raw']['NOIS1']
signal2 = rec2.annotations['abif_raw']['S/N%1']
noise2 = rec2.annotations['abif_raw']['NOIS1']

avrg_signal1 = mean(signal1)
avrg_noise1 = mean(noise1)
avrg_signal2 = mean(signal2)
avrg_noise2 = mean(noise2)

ratio1 = avrg_signal1 / avrg_noise1
ratio2 = avrg_signal2 / avrg_noise2

print('Среднее соотношение сигнала и шума по каждой из двух последовательностей равно: {:.3f} и {:.3f}'.format(ratio1, ratio2))