Mercurial > repos > mvdbeek > mismatch_frequencies
diff mismatch_frequencies.xml @ 22:942464ea4211
remove heads
author | Marius van den Beek <m.vandenbeek@gmail.com> |
---|---|
date | Sun, 24 May 2015 17:29:41 +0200 |
parents | 11a1dae1ab8e |
children | ca7b7890ed20 |
line wrap: on
line diff
--- a/mismatch_frequencies.xml Sun May 24 17:26:57 2015 +0200 +++ b/mismatch_frequencies.xml Sun May 24 17:29:41 2015 +0200 @@ -1,11 +1,11 @@ -<tool id="mismatch_frequencies" name="Mismatch Frequencies" version="0.0.9" hidden="false" > - <description>Analyze mismatch frequencies in BAM/SAM alignments</description> - <requirements> - <requirement type="package" version="0.7.7">pysam</requirement> - <requirement type="package" version="0.14.1">pandas</requirement> - <requirement type="package" version="1.2.1">matplotlib</requirement> - </requirements> - <command interpreter="python">mismatch_frequencies.py --input +<tool id="mismatch_frequencies" name="Mismatch Frequencies" version="0.0.3" hidden="false" > + <description>Analyze mismatch frequencies in BAM/SAM alignments</description> + <requirements> + <requirement type="package" version="0.7.7">pysam</requirement> + <requirement type="package" version="0.14">pandas</requirement> + <requirement type="package" version="1.4">matplotlib</requirement> + </requirements> + <command interpreter="python">mismatch_frequencies.py --input #for i in $rep "$i.input_file" #end for @@ -17,73 +17,20 @@ --n_mm $number_of_mismatches --five_p $five_p --three_p $three_p - --expanded_output_tab $expanded_tab - --possible_mismatches $possible_mismatches - </command> - <inputs> - <repeat name="rep" title="alignment files"> - <param name="input_file" type="data" format="bam,sam" label="Alignment file" help="The input alignment file(s) for which to analyze the mismatches."/> - </repeat> - <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"/> - <param name="possible_mismatches" label="Specify mismatches that should be counted" help="Ignores mismatches that are not listed" type="text" value="AC AG AT CA CG CT GA GC GT TA TC TG"> - <validator type="expression" message="Allowed values are AGCTN, seperated by space.">len([False for char in value if not char in " AGCTN"]) == 0</validator> - </param> - <param name="min_length" label="Minumum read length to analyse" type="integer" value="21"/> - <param name="max_length" label="Maximum read length to analyse" type="integer" value="21"/> - <param name="five_p" label="Ignore mismatches in the first N nucleotides of a read" type="integer" value="0"/> - <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"/> - <param help="Output expanded tabular format" label="Nucleotide mismatches per reference sequence" name="expanded" type="select"> - <option select="true" value="false">No</option> - <option value="expanded">Yes</option> - </param> - </inputs> - <outputs> - <data format="tabular" name="output_tab" /> - <data format="fasta" name="expanded_tab"> - <filter> expanded == "expanded"</filter> - </data> - <data format="pdf" name="output_pdf" /> - </outputs> - <tests> - <test> - <param name="rep_0|input_file" value="3mismatches_ago2ip_s2.bam" ftype="bam" /> - <param name="rep_1|input_file" value="3mismatches_ago2ip_ovary.bam" ftype="bam" /> - <param name="number_of_mismatches" value="1" /> - <param name="min_length" value="21" /> - <param name="max_length" value="21" /> - <param name="three_p" value="0" /> - <param name="five_p" value="0" /> - <output name="tabular" file="mismatch.tab" ftype="tabular"/> - <!-- - <output name="pdf" file="mismatch.pdf" ftype="pdf"/> - --> - </test> - </tests> - <help> + </command> + <inputs> + <repeat name="rep" title="alignment files"> + <param name="input_file" type="data" format="bam,sam" label="Alignment file" help="The input alignment file(s) for which to analyze the mismatches."/> + </repeat> + <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"/> + <param name="min_length" label="Minumum read length to analyse" type="integer" value="21"/> + <param name="max_length" label="Maximum read length to analyse" type="integer" value="21"/> + <param name="five_p" label="Ignore mismatches in the first N nucleotides of a read" type="integer" value="0"/> + <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"/> + </inputs> + <outputs> + <data format="pdf" name="output_pdf" /> + <data format="tabular" name="output_tab" /> + </outputs> -.. class:: infomark - - -***What it does*** - -This tool reconstitues for each aligned read of an alignment file in SAM/BAM format whether -a mismatch is annotated in the MD tag, and if that is the case counts the identity of the -mismatch relative to the reference sequence. The output is a PDF document with the calculated -frequency for each mismatch that occured relative to the total number of valid reads and a table -with the corresponding values. Read length can be limited to a specific read length, and 5 prime and -3 prime-most nucleotides of a read can be ignored. - ----- - -.. class:: warningmark - -***Warning*** - -This tool skips all read that have insertions and has been tested only with bowtie and bowtie2 -generated alignment files. - -Written by Marius van den Beek, m.vandenbeek at gmail . com - </help> - <citations> - </citations> </tool>