1
|
1 <tool id="misa2gff_1" name="MISA to GFF3">
|
|
2 <description>convert MISA to gvf/gff3</description>
|
|
3 <command >cat $inputFile | awk '(NR > 1) {OFS="\t";print $1,"MISA","microsatellite",$6,$7,".",".",".","ID="$1":MISA:microsatellite:"$6";SSR="$4,";Type="$3}' > $outputfile</command>
|
|
4 <inputs>
|
|
5 <param format="tabular" name="inputFile" type="data" label="Input MISA File" help="MISA file from MISA.pl" />
|
|
6 </inputs>
|
|
7 <outputs>
|
|
8 <data format="gff3" name="outputfile" />
|
|
9
|
|
10 </outputs>
|
|
11
|
|
12 <help>
|
|
13
|
|
14 .. class:: infomark
|
|
15
|
|
16 **TIP**
|
|
17
|
|
18 This tool parses MISA columnar output to valid GFF3 format
|
|
19
|
|
20 -------------------
|
|
21
|
|
22 **Example**
|
|
23
|
|
24 *input*
|
|
25
|
|
26 ::
|
|
27
|
|
28 ID SSR nr. SSR type SSR size start end
|
|
29 FW9DLND01C1Y2S 1 p2 (AC)11 22 89 110
|
|
30 FW9DLND01CSIIJ 1 p2 (AC)7 14 226 239
|
|
31 FW9DLND01DOGBK 1 p1 (A)12 12 162 173
|
|
32 FW9DLND01B7MCF 1 p2 (AT)7 14 22 35
|
|
33 FW9DLND01DHCGH 1 p2 (AT)7 14 108 121
|
|
34
|
|
35 *output*
|
|
36
|
|
37 ::
|
|
38
|
|
39 FW9DLND01C1Y2S MISA microsatellite 89 110 . . . ID=FW9DLND01C1Y2S:MISA:microsatellite:89;SSR=(AC)11 ;Type=p2
|
|
40 FW9DLND01CSIIJ MISA microsatellite 226 239 . . . ID=FW9DLND01CSIIJ:MISA:microsatellite:226;SSR=(AC)7 ;Type=p2
|
|
41 FW9DLND01DOGBK MISA microsatellite 162 173 . . . ID=FW9DLND01DOGBK:MISA:microsatellite:162;SSR=(A)12 ;Type=p1
|
|
42 FW9DLND01B7MCF MISA microsatellite 22 35 . . . ID=FW9DLND01B7MCF:MISA:microsatellite:22;SSR=(AT)7 ;Type=p2
|
|
43 FW9DLND01DHCGH MISA microsatellite 108 121 . . . ID=FW9DLND01DHCGH:MISA:microsatellite:108;SSR=(AT)7 ;Type=p2
|
|
44
|
|
45
|
|
46 </help>
|
|
47 </tool>
|
|
48
|