Mercurial > repos > yqiancolumbia > ctk_test
annotate parseAlignment.xml @ 0:0475e4175855 draft default tip
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
author | yqiancolumbia |
---|---|
date | Mon, 30 Apr 2018 05:25:11 -0400 |
parents | |
children |
rev | line source |
---|---|
0
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
1 <tool id="parseAlignment" name="Parse alignment"> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
2 <description>in SAM format to BED format</description> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
3 |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
4 <command interpreter="perl"> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
5 /home/galaxy/tools/CTK/parseAlignment.pl -v |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
6 #if $MAPQ.MAPQRequired == "yes": |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
7 --map-qual $MAPQ.MapQual |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
8 #end if |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
9 #if $MinLen.MinLenRequired == "yes": |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
10 --min-len $MinLen.MinimalReadLength |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
11 #end if |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
12 #if $IndeltoEnd.IndeltoEndRequired == "yes": |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
13 --indel-to-end $IndeltoEnd.NucleotidesfromIndeltoEnd |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
14 #end if |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
15 $SplitDel $IndelinScore --mutation-file $outputMutationFile $input $outputBedFile |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
16 </command> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
17 |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
18 <inputs> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
19 <param name="input" format="sam" type="data" label="Input SAM file (.gz file accepted)" /> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
20 <conditional name="MAPQ"> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
21 <param name="MAPQRequired" type="select" label="Min MAPQ score (e.g. to keep only uniquely mapped reads)"> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
22 <option value="yes" selected="True">Yes</option> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
23 <option value="no">No</option> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
24 </param> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
25 <when value="yes"> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
26 <param name="MapQual" type="integer" value="1" label="Input an integer" /> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
27 </when> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
28 <when value="no"> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
29 </when> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
30 </conditional> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
31 |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
32 <conditional name="MinLen"> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
33 <param name="MinLenRequired" type="select" label="Minimal read length to report"> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
34 <option value="yes" selected="True">Yes</option> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
35 <option value="no">No</option> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
36 </param> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
37 <when value="yes"> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
38 <param name="MinimalReadLength" type="integer" value="18" label="Input an integer" /> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
39 </when> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
40 <when value="no"> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
41 </when> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
42 </conditional> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
43 |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
44 <conditional name="IndeltoEnd"> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
45 <param name="IndeltoEndRequired" type="select" label="Nucleotides from indel to end"> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
46 <option value="yes">Yes</option> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
47 <option value="no" selected="True">No</option> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
48 </param> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
49 <when value="yes"> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
50 <param name="NucleotidesfromIndeltoEnd" type="integer" value="" label="Input an integer" /> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
51 </when> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
52 <when value="no"> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
53 </when> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
54 </conditional> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
55 |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
56 <param name="SplitDel" type="boolean" truevalue="--split-del" falsevalue="" checked="no" label="Split oligo deletion into single nucleotides"/> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
57 <param name="IndelinScore" type="boolean" truevalue="--indel-in-score" falsevalue="" checked="no" label="Count indels in as mismatches reported in the score column"/> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
58 </inputs> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
59 |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
60 <outputs> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
61 <data name="outputBedFile" format="bed" label="Parse tag bed file on ${on_string}"></data> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
62 <data name="outputMutationFile" format="tabular" label="Parse tag mutation file on ${on_string}"></data> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
63 </outputs> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
64 |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
65 <help> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
66 |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
67 .. class:: infomark |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
68 |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
69 **What this tool does** |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
70 |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
71 This tool will parse alignment in SAM format to BED format (both tags and mutations). |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
72 |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
73 It will take as input files in SAM format after alignment and output BED files that keeps only unique mappings (with MAPQ >=1, specific for BWA) and a minimal mapping size of 18 nt. It will also generate mutation files. |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
74 |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
75 Note 1: In the tag bed file, the 5′ column records the number of mismatches (substitutions) in each read. |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
76 |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
77 Note 2: The parsing script relies on MD tags, which is an optional field without strict definition in SAM file format specification. Some aligners might have a slightly different format than bwa in how they report mismatches. |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
78 |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
79 Note 3: Other aligners might not use a positive MAPQ as an indication of unique mapping. Another useful option is "nucleotides from indel to end", which specifies the number of nucleotides towards the end from which indels should not be called (default=5 nt). |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
80 |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
81 </help> |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
82 |
0475e4175855
planemo upload commit 81ece2551cea27cbd0e718ef5b7a2fe8d4abd071-dirty
yqiancolumbia
parents:
diff
changeset
|
83 </tool> |