|
126
|
1 Blast2html: Blast XML to HTML conversion tool
|
|
|
2 =============================================
|
|
|
3
|
|
|
4 This tool accepts Blast XML as input, and creates an HTML page with a human readable version of the result. The output includes graphical displays of where a sequence matches a target and metrics on the quality of the match.
|
|
|
5
|
|
|
6 Example output:
|
|
|
7
|
|
|
8 .. image:: images/blast2html_screenshot1.png
|
|
|
9 :alt: Screenshot of sample output
|
|
|
10 :align: center
|
|
|
11
|
|
|
12 The output format is based on the graphical summary pages on the `NCBI Blast website`_. But as the code for generating such pages is not available we implemented our own clone.
|
|
|
13
|
|
|
14 .. _`NCBI Blast website`: http://blast.ncbi.nlm.nih.gov/Blast.cgi
|
|
|
15
|
|
|
16
|
|
|
17 Galaxy configuration
|
|
|
18 --------------------
|
|
|
19
|
|
|
20 This tool generates an HTML page from Blast XML input. It can be used as a stand-alone command line tool, but it is meant to be used in workflows in Galaxy_. To use this tool in Galaxy some configuration may be required.
|
|
|
21
|
|
|
22 By default Galaxy strips a lot of style information from HTML files before showing them on the screen. The results of this tool will still be somewhat useful but they will be a lot uglier without the styles. You can still download the original generated HTML of course. If you want to view the full result within Galaxy you need to disable this stripping of style information by adding the following line to your Galaxy's ``galaxy.ini``::
|
|
|
23
|
|
|
24 sanitize_all_html = False
|
|
|
25
|
|
|
26 In order to generate links to a gene bank in the result, you will need to tell blast2html what gene bank to use. If you do not configure a gene bank the result will contain links to the NCBI gene bank, and as a name for the gene bank the generic "Gene Bank" is used. These links will only work if the database you are using uses the same accession codes as the NCBI gene bank. To show links to a different gene bank this needs to be configured in Galaxy.
|
|
|
27
|
|
|
28 As at the moment this tool was built Galaxy did not have a good way to specify such configuration directives to tools, we use a slightly ugly solution. The gene bank configuration are added to the NCBI BLAST+ database definition files in Galaxy's ``tool-data`` directory. This directory should contain some files named ``blastdb.loc``, ``blastdb_p.loc``, ``blastdb_d.loc`` etc. that define the databases that the NCBI BLAST+ tool knows about. For BLAST, these files contain three columns of tab-separated data that define the databases. The normal format is::
|
|
|
29
|
|
|
30 <unique_id> <database_caption> <base_name_path>
|
|
|
31
|
|
|
32 with tabs separating the three components. For blast2html we extend this format with two columns that contain a human-readable gene bank name, and a link template. So each row then contains these tab-separated items::
|
|
|
33
|
|
|
34 <unique_id> <database_caption> <base_name_path> <genebank_name> <genebank_link_template>
|
|
|
35
|
|
|
36 NCBI Blast+ will ignore the extra fields in the file.
|
|
|
37
|
|
|
38 So, for example, for a database named ``nt`` with path ``/depot/data2/galaxy/blastdb/nt/nt.chunk`` that uses NCBI nucleontide database accession codes you can use the following definition line::
|
|
|
39
|
|
|
40 nt_02_Dec_2009 nt 02 Dec 2009 /depot/data2/galaxy/blastdb/nt/nt.chunk NCBI Gene Bank http://www.ncbi.nlm.nih.gov/nucleotide/{accession}?report=genbank&log$=nuclalign
|
|
|
41
|
|
|
42 The syntax for the link template is the same as that used by the ``--genelink-template`` command line option (see next section): It can contain the following replacement elements: ``{id[N]}``, ``{fullid}``, ``{defline[N]}``, ``{fulldefline}``, ``{accession}``, where ``N`` is a number. ``id[N]`` and ``defline[N]`` will be replaced by the Nth element of the id or defline, where '``|``' is the field separator.
|
|
|
43
|
|
|
44
|
|
|
45 .. _Galaxy: https://www.galaxyproject.org/
|
|
|
46
|
|
|
47
|
|
|
48 Command line usage
|
|
|
49 ------------------
|
|
|
50
|
|
|
51 ::
|
|
|
52
|
|
|
53 usage: ./blast2html.py [-i] INPUT [-o OUTPUT] [--genelink-template URL_TEMPLATE] [--dbname DBNAME]
|
|
|
54
|
|
|
55 Convert a BLAST XML result into a nicely readable html page
|
|
|
56
|
|
|
57 positional arguments:
|
|
|
58 INPUT The input Blast XML file, same as -i/--input
|
|
|
59
|
|
|
60 optional arguments:
|
|
|
61 -h, --help show this help message and exit
|
|
|
62 -i INPUT, --input INPUT
|
|
|
63 The input Blast XML file
|
|
|
64 -o OUTPUT, --output OUTPUT
|
|
|
65 The output html file
|
|
|
66 --template TEMPLATE The template file to use. Defaults to
|
|
|
67 blast_html.html.jinja
|
|
|
68 --dbname DBNAME The link text to use for external links to a gene bank
|
|
|
69 database. Defaults to 'Gene Bank'
|
|
|
70 --genelink-template URL_TEMPLATE
|
|
|
71 A link template to link hits to a gene bank webpage.
|
|
|
72 The template string is a Python format string. It can
|
|
|
73 contain the following replacement elements: {id[N]},
|
|
|
74 {fullid}, {defline[N]}, {fulldefline}, {accession},
|
|
|
75 where N is a number. id[N] and defline[N] will be
|
|
|
76 replaced by the Nth element of the id or defline,
|
|
|
77 where '|' is the field separator. The default is 'http
|
|
|
78 ://www.ncbi.nlm.nih.gov/nucleotide/{accession}?report=
|
|
|
79 genbank&log$=nuclalign', which is a link to the NCBI
|
|
|
80 nucleotide database.
|
|
|
81 --db-config-dir DB_CONFIG_DIR
|
|
|
82 The directory where databases are configured in
|
|
|
83 blastdb*.loc files. These files are consulted for
|
|
|
84 creating a gene bank link. The files should conform to
|
|
|
85 the format that Galaxy's BLAST expect, i.e. tab-
|
|
|
86 separated tables (with lines starting with '#'
|
|
|
87 ignored), with two extra fields, for a total of five
|
|
|
88 fields per line instead of three.. The third field of
|
|
|
89 each line should be a database path as used by BLAST.
|
|
|
90 The fourth field is the human readable database name,
|
|
|
91 and the fifth a template link to the gene bank
|
|
|
92 conforming to the syntax for the --genelink-template
|
|
|
93 option. Entries in these config files override links
|
|
|
94 specified using --genelink-template and --dbname.
|
|
|
95
|
|
|
96
|
|
|
97 Authorship
|
|
|
98 ----------
|
|
|
99
|
|
|
100 This tool was created and published by `The Hyve B.V.`_ open source bioinformatics solutions.
|
|
|
101
|
|
|
102 .. image:: images/thehyve_logo.png
|
|
|
103 :alt: The Hyve
|
|
|
104 :align: center
|
|
|
105 :target: http://thehyve.nl
|
|
|
106
|
|
|
107 .. _`The Hyve B.V.`: http://thehyve.nl/
|
|
|
108
|
|
|
109
|
|
|
110 Licensing information
|
|
|
111 ---------------------
|
|
|
112
|
|
|
113 Blast2html is free software: you can redistribute it and/or modify
|
|
|
114 it under the terms of the GNU General Public License as published by
|
|
|
115 the Free Software Foundation, either version 3 of the License, or
|
|
|
116 (at your option) any later version.
|
|
|
117
|
|
|
118 This program is distributed in the hope that it will be useful,
|
|
|
119 but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
|
120 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
|
121 GNU General Public License for more details.
|
|
|
122
|
|
|
123 You should have received a copy of the GNU General Public License
|
|
|
124 along with this program. If not, see <http://www.gnu.org/licenses/>.
|
