I am trying to parse blast-xml file from blastp search using searchIO in Biopython. However, I am not sure which one I should use for extracting % identity.
In blast-tab file, I can use hsp.pident, does anyone know the equivalent attribute name for blast-xml. Or I should derive myself, hsp.ident_num/hsp.aln_span*100. Thanks
•
written
by
helenhvalask •
1 answer
The BLAST XML output format does not contain the percentage identify as an explicit field, so yes, you must calculate it from the number of identities and the alignment length.
See for example my BLAST XML to tabular conversion script: https://github.com/peterjc/galaxy_blast/blob/master/tools/ncbi_blast_plus/blastxml_to_tabular.py#L230
(Note if you are using Python 2, beware of integer division if you do the calculation as currently written!)
Log in to answer this question.
-
ERROR When Using Ensembl Plants Biomart Wget API
written by Dora •
-
parsing blast xml to MSA/PSSM
written by patc •
-
AttributeError: 'Hit' object has no attribute 'sacc_ver'
written by bioinfoSeeker •
-
Blast parse script: ValueError: Required query and/or hit ID field not found.
written by bioinfoSeeker •
-
TypeError: 'generator' object is unsubscriptable in biopython
written by bio90029 •
-
parsing hammer 2.3.2 results using biopython
written by frcamacho
-
Processing Blast XML output generated by GNU Parallel with Biopython SearchIO
written by mossmatters •
-
Attribute Error Using Biopython Bio.Searchio.Convert
written by charles.bridges
-
Parsing Uniprot Blast Output
written by jcastrofigueroa
-
Extracting Fasta Alignments From Parsed Blastxml File
written by Eric •
