from Seq1 import Seq
PROJECT_PATH = "./PROJECT/"
gene_list = ["U5", "ADA", "FRAT1", "FXN"]
print('----| Exercise 10 |----')
s1 = Seq()
for gene in gene_list:
sequence = Seq()
sequence.read_fasta(PROJECT_PATH + gene + '.txt')
print('Gene', gene, ':', 'Most frequent Base:', sequence.frequent_base())
from Seq1 import Seq
def print_result(i, sequence):
print('Sequence ' + str(1) + ': (Length ' + str(sequence.len()) + ' ) ' +
str(sequence))
print('Bases:', sequence.count())
print('Rev:', sequence.reverse())
print('Comp: ', sequence.complement())
PROJECT_PATH = '../P0/sequences/'
print('-----| Practice 1, Exercise 10|------')
s1 = Seq()
s1.read_fasta(PROJECT_PATH + 'ADA.txt')
print('Gene', 'ADA.txt', ': Most frequent Base:', Seq.frequent_base(s1))
