usage = """
%prog $Filename $Taxonomy
"""
# argv
fi = sys.argv[1] # Filename, fi = 'entrezgene.txt'
ti = sys.argv[2] # Taxonomy, ti = '10116'
# load
Entrez.email = "*****@*****.**"
handle = Entrez.efetch(db="Taxonomy", id=ti, retmode="xml")
records = Entrez.read(handle)
si = records[0]["ScientificName"] # Specics_Id, si = 'Rattus norvegicus'
mg = mygene.MyGeneInfo()
kg = KEGG()
kt = str(kg.lookfor_organism(si)).replace("u'",
"").split(" ")[1] # KEGG_Taxonomy
kge = kg.list(kt)
# init
lq = [] # List_Query
if fi == "-":
for ln in sys.stdin:
lq.append(ln.strip('\n'))
else:
f = open(fi, 'r')
for ln in f:
lq.append(ln.strip('\n'))
lf = [
'entrezgene', 'ensembl.gene', 'symbol', 'name', 'alias', 'summary',
'refseq', 'unigene', 'ensembl.transcript', 'ensembl.protein', 'uniprot',
