from seq1 import Seq
print("-- Practice 1, Exercise 9 --")
folder = "../session-04/"
FILENAME =folder + "U5.txt"
# -- Create a Null sequence
s = Seq()
# -- Initialize the null seq with the given file in fasta format
s.read_fasta(FILENAME)
print(f"Sequence : (Length: {s.len()}) {s}")
print("\tBases: ", s.count())
print("\tRev: ", s.reverse())
print("\tComp:", s.complement())
from seq1 import Seq
# -- Create a Null sequence
seq_1 = Seq()
# -- Create a valid sequence
seq_2 = Seq("ACTGA")
# -- Create an invalid sequence
seq_3 = Seq("Invalid sequence")
print("Practice 1, Exercise 8")
BASES = ["A", "G", "T", "C"]
print("Sequence 1:(length: ", seq_1.len(), ")", seq_1)
print("Bases:", seq_1.count())
print("Rev:", seq_1.reverse())
print("Comp:", seq_1.complement())
print("\nSequence 2: (length:", seq_2.len(), ")", seq_2)
print("Bases:", seq_2.count())
print("Rev", seq_2.reverse())
print("Comp:", seq_2.complement())
print("\nSequence 3: (length:", seq_3.len(), ")", seq_3)
print("Bases:", seq_3.count())
print("Rev:", seq_3.reverse())
print("comp:", seq_3.complement())
# INFO command
elif "INFO" in msg:
seq_info = Seq(argument_command)
count_bases_string = ""
for base, count in seq_info.count().items():
s_base = str(base) + ": " + str(count) + " (" + str(
round(count / seq_info.len() * 100, 2)) + "%)" + "\n"
count_bases_string += s_base
response = ("Sequence: " + str(seq_info) + "\n" +
"Total length: " + str(seq_info.len()) + "\n" +
count_bases_string)
elif "COMP" in msg:
seq_comp = Seq(argument_command)
response = seq_comp.complement() + "\n"
elif "REV" in msg:
seq_rev = Seq(argument_command)
response = seq_rev.reverse() + "\n"
elif "GENE" in msg:
gene = argument_command
s = Seq()
s.read_fasta("../Session-04/" + gene + ".txt")
response = str(s) + "\n"
# -- The message has to be encoded into bytes
# Server Console
termcolor.cprint(msg[:msg.find(" ")], "green")