from Bio import Entrez from Bio import SeqIO #Aceder ao ncbi e guardar o ficheiro correspondente a zona do genoma do grupo 4 Entrez.email = "a68367@alunos.uminho.pt" handle = Entrez.efetch(db="nucleotide", rettype="gb", retmode="text", id="52840256", seq_start="693001", seq_stop="934905") seq_record = SeqIO.read(handle, "gb") #Criar ficheiro SeqIO.write(seq_record, "sequence.gb", "genbank") handle.close() #Passar os dados para a consola record = SeqIO.read("sequence.gb", "genbank") print "NCBI" + "\n" print record