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)))
Exemple #2
0
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))
Exemple #3
0
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)
Exemple #4
0
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)))
Exemple #6
0
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])