Python get_taxid示例

示例#1

文件： test_BlastTools.py 项目： hobrien/Python

def test_get_taxid():
  assert get_taxid('gi|455475748|gb|KC437649.1|') == P_venosa_id 
  assert get_taxid('gi|455475748|gb|KC437649.1') == P_venosa_id 
  assert get_taxid('KC437649') == P_venosa_id
  assert get_taxid('KC437649.1') == P_venosa_id
  assert get_taxid('gi|455475748') == P_venosa_id 
  assert get_taxid('455475748') == P_venosa_id 

示例#2

文件： GetTaxID.py 项目： hobrien/Python

#!/usr/local/bin/python

import fileinput
from BlastTools import get_taxid

  
for line in fileinput.input():
  if "Hit_ID" in line: continue
  fields = line.split()
  taxid =  get_taxid(fields[2])
  if taxid > 0:
    print '\t'.join(fields[0], taxid)