|
18
|
1 <tool id="mismatch_frequencies" name="Mismatch Frequencies" version="0.0.7" hidden="false" >
|
|
12
|
2 <description>Analyze mismatch frequencies in BAM/SAM alignments</description>
|
|
|
3 <requirements>
|
|
13
|
4 <requirement type="package" version="0.7.7">pysam</requirement>
|
|
12
|
5 <requirement type="package" version="0.14.1">pandas</requirement>
|
|
17
|
6 <requirement type="package" version="1.2.1">matplotlib</requirement>
|
|
12
|
7 </requirements>
|
|
|
8 <command interpreter="python">mismatch_frequencies.py --input
|
|
|
9 #for i in $rep
|
|
|
10 "$i.input_file"
|
|
|
11 #end for
|
|
|
12 --name
|
|
|
13 #for i in $rep
|
|
|
14 "$i.input_file.name"
|
|
|
15 #end for
|
|
|
16 --output_pdf $output_pdf --output_tab $output_tab --min $min_length --max $max_length
|
|
|
17 --n_mm $number_of_mismatches
|
|
|
18 --five_p $five_p
|
|
|
19 --three_p $three_p
|
|
|
20 </command>
|
|
|
21 <inputs>
|
|
16
|
22 <repeat name="rep" title="alignment files">
|
|
|
23 <param name="input_file" type="data" format="bam,sam" label="Alignment file" help="The input alignment file(s) for which to analyze the mismatches."/>
|
|
12
|
24 </repeat>
|
|
|
25 <param name="number_of_mismatches" label="Maximum number of allowed mismatches per read" help="Discard reads with more than the chosen number of mismatches from the frequency calculation" type="integer" value="3"/>
|
|
|
26 <param name="min_length" label="Minumum read length to analyse" type="integer" value="21"/>
|
|
|
27 <param name="max_length" label="Maximum read length to analyse" type="integer" value="21"/>
|
|
|
28 <param name="five_p" label="Ignore mismatches in the first N nucleotides of a read" type="integer" value="0"/>
|
|
|
29 <param name="three_p" label="Ignore mismatches in the last N nucleotides of a read" help="useful to discriminate between tailing events and editing events" type="integer" value="3"/>
|
|
|
30 </inputs>
|
|
|
31 <outputs>
|
|
|
32 <data format="tabular" name="output_tab" />
|
|
|
33 <data format="pdf" name="output_pdf" />
|
|
|
34 </outputs>
|
|
|
35 <tests>
|
|
|
36 <test>
|
|
|
37 <param name="rep_0|input_file" value="3mismatches_ago2ip_s2.bam" ftype="bam" />
|
|
18
|
38 <param name="rep_1|input_file" value="3mismatches_ago2ip_ovary.bam" ftype="bam" />
|
|
12
|
39 <param name="number_of_mismatches" value="1" />
|
|
|
40 <param name="min_length" value="21" />
|
|
|
41 <param name="max_length" value="21" />
|
|
18
|
42 <param name="three_p" value="0" />
|
|
|
43 <param name="five_p" value="0" />
|
|
|
44 <output name="tabular" file="mismatch.tab" ftype="tabular"/>
|
|
|
45 <!--
|
|
|
46 <output name="pdf" file="mismatch.pdf" ftype="pdf"/>
|
|
|
47 -->
|
|
|
48 </test>
|
|
|
49 </tests>
|
|
|
50 <help>
|
|
|
51
|
|
|
52 .. class:: infomark
|
|
|
53
|
|
|
54
|
|
|
55 ***What it does***
|
|
|
56
|
|
|
57 This tool reconstitues for each aligned read of an alignment file in SAM/BAM format whether
|
|
|
58 a mismatch is annotated in the MD tag, and if that is the case counts the identity of the
|
|
|
59 mismatch relative to the reference sequence. The output is a PDF document with the calculated
|
|
|
60 frequency for each mismatch that occured relative to the total number of valid reads and a table
|
|
|
61 with the corresponding values. Read length can be limited to a specific read length, and 5 prime and
|
|
|
62 3 prime-most nucleotides of a read can be ignored.
|
|
|
63
|
|
|
64 ----
|
|
|
65
|
|
|
66 .. class:: warningmark
|
|
|
67
|
|
|
68 ***Warning***
|
|
|
69
|
|
|
70 This tool skips all read that have insertions and has been tested only with bowtie and bowtie2
|
|
|
71 generated alignment files.
|
|
|
72
|
|
|
73 Written by Marius van den Beek, m.vandenbeek at gmail . com
|
|
|
74 </help>
|
|
12
|
75 </tool>
|
