import P0.Seq0 as Seq0 gene_folder = "./sequences/" gene_list = ["U5", "ADA", "FRAT1", "FXN"] print("------|EXERCISE 8|------") for gene in gene_list: seq = Seq0.seq_read_fasta(gene_folder + gene + ".txt") bases_count = Seq0.seq_count(seq) print("Gene " + gene + ": Most frequent base: " + str(Seq0.max_base(bases_count)))
import P0.Seq0 as Seq0 GENE_FOLDER = "./Sequences/" gene_list = ["U5", "ADA", "FRAT1", "FXN"] print("------| Exercise 5 |------") for gene in gene_list: seq = Seq0.seq_read_fasta(GENE_FOLDER + gene + ".txt") print("Gene " + gene + ":", Seq0.seq_count(seq))
import P0.Seq0 as Seq0 gene_folder = "./Sequences/" U5 = "U5" ADA = "ADA" PRAT1 = "PRAT1" FXN = "FXN" base = ["A", "C", "G", "T"] sequence_U5 = Seq0.seq_read_fasta(gene_folder + U5 + ".txt") bases = Seq0.seq_count(sequence_U5) l_U5 = [] for a in bases: l_U5.append(bases[a]) U5_h = l_U5.index(max(l_U5)) sequence_ADA = Seq0.seq_read_fasta(gene_folder + ADA + ".txt") bases2 = Seq0.seq_count(sequence_ADA) l_ADA = [] for b in bases2: l_ADA.append(bases2[b]) ADA_h = l_ADA.index(max(l_ADA)) sequence_PRAT1 = Seq0.seq_read_fasta(gene_folder + PRAT1 + ".txt") bases3 = Seq0.seq_count(sequence_U5) l_PRAT1 = [] for c in bases3: l_PRAT1.append(bases3[c]) PRAT1_h = l_PRAT1.index(max(l_PRAT1)) sequence_FXN = Seq0.seq_read_fasta(gene_folder + FXN + ".txt") bases4 = Seq0.seq_count(sequence_FXN)
import P0.Seq0 as Seq0 gene_folder = "./Sequences/" gene_list = ["U5", "ADA", "PRAT1", "FXN"] print("-----| Exercise 5 |------") for gene in gene_list: sequence = Seq0.seq_read_fasta(gene_folder + gene + ".txt") print("Gene ", gene, ":", Seq0.seq_count(sequence))
import P0.Seq0 as Seq0 gene_folder = "./sequences/" gene_list = ["U5", "ADA", "FRAT1", "FXN"] print("------|EXERCISE 5|------") for gene in gene_list: seq = Seq0.seq_read_fasta(gene_folder + gene + ".txt") print("Gene " + gene + ": " + str(Seq0.seq_count(seq)))
import P0.Seq0 as Seq0 GENE_FOLDER = "./Sequences/" gene_list = ["U5", "ADA", "FRAT1", "FXN"] base_list = ["A", "C", "T", "G"] print("------| Exercise 8 |------") for gene in gene_list: seq = Seq0.seq_read_fasta(GENE_FOLDER + gene + ".txt") sort_freq = sorted(Seq0.seq_count(seq).items(), key=lambda item: item[1], reverse=True) print("Gene " + gene + " most frequent base:", sort_freq[0][0])