Mercurial > repos > jankanis > blast2html
view test-data/blast xml example2.xml @ 120:2729c2326235
Fix for Rikilt issue 13
Hit e-value and identity% should be taken from the hsp with the highest
bit score. Previously each of these values was calculated independently.
Also use arrays for cover calculation instead of python lists and
refactor the hit_info() code a bit.
author | Jan Kanis <jan.code@jankanis.nl> |
---|---|
date | Thu, 31 Jul 2014 16:14:36 +0200 |
parents | 420d0508dcd6 |
children |
line wrap: on
line source
<?xml version="1.0"?> <!DOCTYPE BlastOutput PUBLIC "-//NCBI//NCBI BlastOutput/EN" "http://www.ncbi.nlm.nih.gov/dtd/NCBI_BlastOutput.dtd"> <BlastOutput> <BlastOutput_program>blastn</BlastOutput_program> <BlastOutput_version>BLASTN 2.2.28+</BlastOutput_version> <BlastOutput_reference>Stephen F. Altschul, Thomas L. Madden, Alejandro A. Sch&auml;ffer, Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997), "Gapped BLAST and PSI-BLAST: a new generation of protein database search programs", Nucleic Acids Res. 25:3389-3402.</BlastOutput_reference> <BlastOutput_db>/share/BlastDB/nt</BlastOutput_db> <BlastOutput_query-ID>Query_1</BlastOutput_query-ID> <BlastOutput_query-def>No definition line</BlastOutput_query-def> <BlastOutput_query-len>16</BlastOutput_query-len> <BlastOutput_param> <Parameters> <Parameters_expect>0.001</Parameters_expect> <Parameters_sc-match>2</Parameters_sc-match> <Parameters_sc-mismatch>-3</Parameters_sc-mismatch> <Parameters_gap-open>5</Parameters_gap-open> <Parameters_gap-extend>2</Parameters_gap-extend> <Parameters_filter>L;m;</Parameters_filter> </Parameters> </BlastOutput_param> <BlastOutput_iterations> <Iteration> <Iteration_iter-num>1</Iteration_iter-num> <Iteration_query-ID>Query_1</Iteration_query-ID> <Iteration_query-def>No definition line</Iteration_query-def> <Iteration_query-len>16</Iteration_query-len> <Iteration_hits> </Iteration_hits> <Iteration_stat> <Statistics> <Statistics_db-num>28813188</Statistics_db-num> <Statistics_db-len>53027764389</Statistics_db-len> <Statistics_hsp-len>13</Statistics_hsp-len> <Statistics_eff-space>157959578835</Statistics_eff-space> <Statistics_kappa>0.41</Statistics_kappa> <Statistics_lambda>0.625</Statistics_lambda> <Statistics_entropy>0.78</Statistics_entropy> </Statistics> </Iteration_stat> <Iteration_message>No hits found</Iteration_message> </Iteration> </BlastOutput_iterations> </BlastOutput>