biopython_parse_ncbiblast.py

#!/usr/bin/env python3

from Bio.Blast import NCBIXML
result_handle = open("../files/UTKBKAM5014-Alignment.xml")
blast_records = NCBIXML.parse(result_handle)
for blast_record in blast_records:
   query_id = blast_record.query_id
   for alignment in blast_record.alignments:
     for hsp in alignment.hsps:
        if hsp.expect < 1e-10:
           print('qid:' , query_id , 'hit_id:' , alignment.title, 'E:' , hsp.expect )