"RBMY2YP": "ENSG00000227633",
"FGFR3": "ENSG00000068078",
"KDR": "ENSG00000128052",
"ANK2": "ENSG00000145362"
}
SERVER = "rest.ensembl.org"
ENDPOINT = "/sequence/id/"
PARAMETERS = "?content-type=application/json"
connection = http.client.HTTPConnection(SERVER)
try:
user_gene = input("Enter the gene that you want to analyze:")
id = DICT_GENE[user_gene]
connection.request("GET", ENDPOINT + id + PARAMETERS)
response = connection.getresponse()
if response.status == 200:
response_dict = json.loads(response.read().decode())
sequence = Seq(response_dict["seq"])
s_length = sequence.len()
percentages = sequence.percentage_base(sequence.count_base())
most_frequent_base = sequence.frequent_base(sequence.count())
print("Total length:", s_length)
for value in percentages.values():
print(value)
print("Most frequent base:", most_frequent_base)
except KeyError:
print("The key is not inside our dictionary. Choose one of the following:",
list(DICT_GENE.keys))
from seq import Seq
import operator
print("-------Exercise 10-------")
FASTA_GENE_SEQ = ["U5.txt", "FRAT1.txt", "FXN.txt", "RNU6_269P.txt", "ADA.txt"]
for fasta in FASTA_GENE_SEQ:
s = Seq()
s.read_fasta_format(fasta)
var = s.count()
max_key = max(var, key=var.get)
print(max_key)
