Quantcast
Channel: Post Feed
Viewing all articles
Browse latest Browse all 41826

How To Get Number Of Hit For A Query??

$
0
0

hi , i'm using biopython to parse my blast xml output i make condition print on the number of for query. i print my query and sbject if the query have only one hit. but it's not working here is my script please help. thanks

result_handle = open("result.xml")
blast_records = NCBIXML.parse(result_handle)
 for blast_record in blast_records:
    for alignment in blast_record.alignments:
           for hsp in alignment.hsps:
              if(hsp.num_alignments == 1):
                   print "Query:",blast_record.query
                   print "subjet:",alignment.hit_def
                   print "length:", alignment.length
                   print "nombre de gaps:",hsp.gaps

Viewing all articles
Browse latest Browse all 41826

Trending Articles