Готов подсчитать, но мне нужны сами фрагменты. Пожалуйста, пришлите последовательности фрагментов A, B и C. Укажите формат данных: - являются ли фрагменты двуцепочечными (две последовательности на фрагмент: одна для каждой цепи) или одиночными (одна цепь)? - если двуцепочочные, пришлите пары цепей по каждому фрагменту (например, A1 и A2 для фрагмента A и т. д.). Что нужно посчитать для каждого фрагмента - общее количество нуклеотидов каждого типа: a, t, g, c - количество пар A–T и пар G–C (base pairs) внутри фрагмента - отношения: a/t, g/c и (a+t)/(g+c) Как буду считать 1) Если фрагмент дан как одиночная последовательность (одна цепь): - a = число 'A' (и/или 'a'), t = число 'T', g = число 'G', c = число 'C'. - количество пар A–T и G–C следует трактовать как максимально возможные пары на двойной спирали: - AT-пары ≈ min(a, t) - GC-пары ≈ min(g, c) - далее вычисляются ratios: a/t, g/c, (a+t)/(g+c) (с обработкой деления на ноль). 2) Если фрагменты даны как двуцепочочные (s1 и s2 одинаковой длины): - a = количество A в обеих цепях (sum of 'A' в s1 и в s2), t, g, c аналогично. - AT-пары = количество позиций i, где (s1[i], s2[i]) образуют A–T или T–A - GC-пары = количество позиций i, где (s1[i], s2[i]) образуют G–C или C–G - затем ratios: a/t, g/c, (a+t)/(g+c) Если хотите, могу привести сразу пример формата вывода после получения данных. Пример кода (Python), чтобы вы могли просто прислать данные и получить таблицу результатов - Вариант 1: одиночные фрагменты (одна цепь на фрагмент) def analyze_single(str_seq): s = str_seq.upper() a = s.count('A') t = s.count('T') g = s.count('G') c = s.count('C') at_pairs = min(a, t) # предполагаемая максимальная число А–Т пар gc_pairs = min(g, c) # максимальные Г–Ц пары def safe_ratio(x, y): return x / y if y != 0 else None ratios = { 'a/t': safe_ratio(a, t), 'g/c': safe_ratio(g, c), '(a+t)/(g+c)': safe_ratio(a + t, g + c) } return {'a': a, 't': t, 'g': g, 'c': c, 'AT_pairs': AT_pairs, 'GC_pairs': GC_pairs, 'ratios': ratios} # Вариант 2: двуцепочочные фрагменты (s1 и s2) def analyze_double(s1, s2): a = s1.upper().count('A') + s2.upper().count('A') t = s1.upper().count('T') + s2.upper().count('T') g = s1.upper().count('G') + s2.upper().count('G') c = s1.upper().count('C') + s2.upper().count('C') AT_pairs = sum(1 for a1, a2 in zip(s1.upper(), s2.upper()) if (a1, a2) in [('A','T'), ('T','A')]) GC_pairs = sum(1 for b1, b2 in zip(s1.upper(), s2.upper()) if (b1, b2) in [('G','C'), ('C','G')]) def safe_ratio(x, y): return x / y if y != 0 else None ratios = { 'a/t': safe_ratio(a, t), 'g/c': safe_ratio(g, c), '(a+t)/(g+c)': safe_ratio(a + t, g + c) } return {'a': a, 't': t, 'g': g, 'c': c, 'AT_pairs': AT_pairs, 'GC_pairs': GC_pairs, 'ratios': ratios} Puedo adaptar код под ваш формат ввода и вывести таблицу вида: Fragment: A a = ... t = ... g = ... c = ... AT_pairs = ... GC_pairs = ... a/t = ... g/c = ... (a+t)/(g+c) = ... Прошу пришлите: - фрагменты A, B, C (одна цепь каждая или по две цепи, если двуцепочочные), - и формат, в каком хотите увидеть результат (п просто список чисел или таблица).