import_genbank.py 文件源码

def create_profile_hits(cur, feature, gene_id):
    '''Create profile hit entries for a feature'''
    detected_domains = parse_domains_detected(feature)
    for domain in detected_domains:
        domain['gene_id'] = gene_id
        cur.execute("""
SELECT gene_id FROM antismash.profile_hits WHERE
    gene_id = %(gene_id)s AND
    name = %(name)s AND
    evalue = %(evalue)s AND
    bitscore = %(bitscore)s""", domain)
        ret = cur.fetchone()
        if ret is None:
            try:
                cur.execute("""
INSERT INTO antismash.profile_hits (gene_id, name, evalue, bitscore, seeds)
    VALUES (%(gene_id)s, %(name)s, %(evalue)s, %(bitscore)s, %(seeds)s)""", domain)
            except psycopg2.IntegrityError:
                print(feature)
                print(domain)
                raise