# -*- coding:utf-8 -*-
from collections import defaultdict
tiqu_gi = open('tiqu_gi.txt','r')
gi2taxid = open('cutted_nucl_gb.accession2taxid','r')
taxid2name = open('Taxonomy/names.dmp','r')
get_name = open('scientific_name.txt','w')
taxid_name_dict={}
for lines in taxid2name:
if "scientific name" in lines:
line = lines.strip().split("|")
taxid = line[0].strip()
name = line[1].strip()
taxid_name_dict[taxid]=name
tiqu_dict=defaultdict(list)
for lines in tiqu_gi:
line = lines.strip().split("\ ")
gi = line[1].split("|")[1]
tiqu_dict[gi].append("\ ".join(line))
gi_taxid_dict={}
for lines in gi2taxid:
line = lines.strip().split("\ ")