from Bio import SeqIO

 
#recolherfuncoes
#traduzir os genes na sequencia


record = SeqIO.read("sequence.gb", "gb")
func=[]
for feature in record.features:
    if feature.type == "CDS":
        if 'function' in feature.qualifiers:
            func.append((feature.qualifiers["locus_tag"][0], feature.qualifiers["function"][0]))
trad=[]
for feature in record.features:
    if feature.type=="CDS":
        trad.append((feature.qualifiers["locus_tag"][0], feature.qualifiers["translation"][0]))

#auxilio grupo 14

#print de INF importante

print "Gene e Funcao:" + "\n" 
print "Genes com funcao definida: " + str(len(func))
print "Genes com funcao desconhecida: " + str(214-len(func))+ "\n" 
print "Genes e as respetivas funcoes: " + "\n"
for function in func:
    print function
print "\n" + "Gene e Traducao:" + "\n"
for translation in trad:
    print translation