|
0
|
1 <tool id="megablast_xml_parser" name="Parse blast XML output" version="1.0.0">
|
|
|
2 <description></description>
|
|
|
3 <command interpreter="python">megablast_xml_parser.py $input1 $output1</command>
|
|
|
4 <inputs>
|
|
|
5 <param name="input1" type="data" format="blastxml" label="Megablast XML output" />
|
|
|
6 </inputs>
|
|
|
7 <outputs>
|
|
|
8 <data name="output1" format="tabular"/>
|
|
|
9 </outputs>
|
|
|
10 <tests>
|
|
|
11 <test>
|
|
|
12 <param name="input1" value="megablast_xml_parser_test1.gz" ftype="blastxml" />
|
|
|
13 <output name="output1" file="megablast_xml_parser_test1_out.tabular" ftype="tabular" />
|
|
|
14 </test>
|
|
|
15 </tests>
|
|
|
16 <help>
|
|
|
17
|
|
|
18 **What it does**
|
|
|
19
|
|
|
20 This tool processes the XML output of any NCBI blast tool (if you run your own blast jobs, the XML output can be generated with **-m 7** option).
|
|
|
21
|
|
|
22 -----
|
|
|
23
|
|
|
24 **Output fields**
|
|
|
25
|
|
|
26 This tools returns tab-delimited output with the following fields::
|
|
|
27
|
|
|
28 Description Example
|
|
|
29 ----------------------------------------- -----------------
|
|
|
30
|
|
|
31 1. Name of the query sequence Seq1
|
|
|
32 2. Length of the query sequence 30
|
|
|
33 3. Name of target sequence gnl|BL_ORD_ID|0
|
|
|
34 4. Length of target sequence 5528445
|
|
|
35 5. Alignment bit score 59.96
|
|
|
36 6. E-value 8.38112e-11
|
|
|
37 7. Start of alignment within query 1
|
|
|
38 8. End of alignment within query 30
|
|
|
39 9. Start of alignment within target 5436010
|
|
|
40 10. End of alignment within target 5436039
|
|
|
41 11. Query frame 1
|
|
|
42 12. Target frame 1
|
|
|
43 13. Number of identical bases within 29
|
|
|
44 the alignment
|
|
|
45 14. Alignment length 30
|
|
|
46 15. Aligned portion (sequence) of query CGGACAGCGCCGCCACCAACAAAGCCACCA
|
|
|
47 16. Aligned portion (sequence) of target CGGACAGCGCCGCCACCAACAAAGCCATCA
|
|
|
48 17. Midline indicating positions of ||||||||||||||||||||||||||| ||
|
|
|
49 matches within the alignment
|
|
|
50
|
|
|
51 ------
|
|
|
52
|
|
|
53 .. class:: infomark
|
|
|
54
|
|
|
55 Note that this form of output does not contain alignment identify value. However, it can be computed by dividing the number of identical bases within the alignment (Field 13) by the alignment length (Field 14) using *Text Manipulation->Compute* tool
|
|
|
56
|
|
|
57
|
|
|
58
|
|
|
59 </help>
|
|
|
60 </tool>
|
