Exemplo n.º 1
0
def main_blast_n():
    an_seq=Analise_Seq('NC_000919', start_position=241801, end_position=381900)
    gene_ids=an_seq._get_geneid_hipothetical_P()
    #retorna gene_ids que pertencem a proteinas hipoteticas e nao tem gis
    print(len(gene_ids))
    print(gene_ids)
    blast_n(gene_ids)