Mercurial > repos > mvdbeek > bam_readtagger
changeset 41:29734dbad6de draft
planemo upload for repository https://github.com/bardin-lab/readtagger/tree/master/galaxy commit 2e991a52c828b0234020e89748eafe1028e44a7d
| author | mvdbeek |
|---|---|
| date | Sun, 01 Apr 2018 03:38:16 -0400 |
| parents | bb09f4ee2ff8 |
| children | fbfd1082de71 |
| files | add_matesequence.xml allow_dovetailing.xml bam_readtagger.xml deploy.sh filter_insertions.xml findcluster.xml plot_coverage.xml test-data/all_control_insertions.gff test-data/all_treatment_insertions.gff test-data/confirmed_insertions.gff test-data/dm6.bam test-data/putative_insertions.gff test-data/three_cluster_out.bam test-data/three_cluster_out.gff test-data/three_cluster_out.vcf update_mapq.xml write_supplementary_fastq.xml |
| diffstat | 16 files changed, 2053 insertions(+), 24 deletions(-) [+] |
line wrap: on
line diff
--- a/add_matesequence.xml Sat Jan 06 07:21:07 2018 -0500 +++ b/add_matesequence.xml Sun Apr 01 03:38:16 2018 -0400 @@ -1,7 +1,7 @@ -<tool id="add_matesequence" name="Add matesequence" version="0.4.6"> +<tool id="add_matesequence" name="Add matesequence" version="0.4.10"> <description>into tag field</description> <requirements> - <requirement type="package" version="0.4.6">readtagger</requirement> + <requirement type="package" version="0.4.10">readtagger</requirement> </requirements> <version_command>add_matesequence --version</version_command> <command detect_errors="aggressive"><![CDATA[
--- a/allow_dovetailing.xml Sat Jan 06 07:21:07 2018 -0500 +++ b/allow_dovetailing.xml Sun Apr 01 03:38:16 2018 -0400 @@ -1,7 +1,7 @@ -<tool id="allow_dovetailing" name="Allow dovetailing" version="0.4.6"> +<tool id="allow_dovetailing" name="Allow dovetailing" version="0.4.10"> <description>modifies proper_pair flag in bam files</description> <requirements> - <requirement type="package" version="0.4.6">readtagger</requirement> + <requirement type="package" version="0.4.10">readtagger</requirement> </requirements> <command detect_errors="aggressive"><![CDATA[ allow_dovetailing -i '$input' -o '$output'
--- a/bam_readtagger.xml Sat Jan 06 07:21:07 2018 -0500 +++ b/bam_readtagger.xml Sun Apr 01 03:38:16 2018 -0400 @@ -1,10 +1,10 @@ -<tool id="bam_readtagger" name="Tag alignment files" version="0.4.6"> +<tool id="bam_readtagger" name="Tag alignment files" version="0.4.10"> <description>from multiple bam files</description> <macros> <import>macros.xml</import> </macros> <requirements> - <requirement type="package" version="0.4.6">readtagger</requirement> + <requirement type="package" version="0.4.10">readtagger</requirement> </requirements> <command detect_errors="aggressive"><![CDATA[ readtagger -t '$tag_file' -s
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/filter_insertions.xml Sun Apr 01 03:38:16 2018 -0400 @@ -0,0 +1,59 @@ +<tool id="confirm_insertions" name="Confirm or discard insertions" version="0.4.10"> + <description>based on softclips in control files</description> + <requirements> + <requirement type="package" version="0.4.10">readtagger</requirement> + </requirements> + <version_command>confirm_insertions --version</version_command> + <command detect_errors="aggressive"><![CDATA[ + ln -fs $putative_insertions putative_insertions.gff && + ln -fs $all_insertions all_insertions.gff && + ln -fs $all_controls all_controls_path.gff && + confirm_insertions + --putative_insertions_path putative_insertions.gff + --all_treatments_path all_insertions.gff + --all_controls_path all_controls_path.gff + --output_path '$output' + $output_discard + ]]></command> + <inputs> + <param name="putative_insertions" label="Path to filtered putative insertinos" argument="--putative_insertions_path" type="data" format="gff3"/> + <param name="all_insertions" label="Path to all unfiltered insertions" argument="--all_treatments_path" type="data" format="gff3"/> + <param name="all_controls" label="Path to all unfiltered insertions from control sample" argument="--all_controls_path" type="data" format="gff3"/> + <param name="output_discard" label="Output discarded records? Discarded records will be marked as failed." truevalue="--output_discarded_records" falsevalue="--no_output_discarded_records" checked="True" type="boolean"/> + </inputs> + <outputs> + <data name="output" format="gff3"/> + </outputs> + <tests> + <test> + <param name="putative_insertions" value="putative_insertions.gff" ftype="gff3"/> + <param name="all_insertions" value="all_treatment_insertions.gff" ftype="gff3"/> + <param name="all_controls" value="all_control_insertions.gff" ftype="gff3"/> + <output name="output" file="confirmed_insertions.gff" ftype="gff3"/> + </test> + </tests> + <help><![CDATA[ +.. code-block:: + + + Usage: confirm_insertions [OPTIONS] + + Confirm insertions by checking that control file does not contain the same + clipping pattern. + + Options: + -p, --putative_insertions_path PATH + Path to file containing putative somatic + insertions [required] + -t, --all_treatments_path PATH Path to file containing all treatment + insertions [required] + -c, --all_controls_path PATH Path to file containing all control + insertions [required] + -o, --output_path TEXT Write annotated output to this Path. + --output_discarded_records / --no_output_discarded_records + Discard an alternative flag if the current + read is in a proper pair. + --version Show the version and exit. + --help Show this message and exit. + ]]></help> +</tool>
--- a/findcluster.xml Sat Jan 06 07:21:07 2018 -0500 +++ b/findcluster.xml Sun Apr 01 03:38:16 2018 -0400 @@ -1,10 +1,10 @@ -<tool id="findcluster" name="Find clusters of reads" version="0.4.6"> +<tool id="findcluster" name="Find clusters of reads" version="0.4.10"> <description>in bam files</description> <macros> <import>macros.xml</import> </macros> <requirements> - <requirement type="package" version="0.4.6">readtagger</requirement> + <requirement type="package" version="0.4.10">readtagger</requirement> </requirements> <version_command>findcluster --version</version_command> <command detect_errors="aggressive"><![CDATA[ @@ -28,33 +28,57 @@ --genome_bwa_index '$genome_source.ref_file.fields.path' #end if #end if - --output_bam '$output_bam' - --output_gff '$output_gff' - --output_fasta '$output_fasta' + #if str($make_bam) == 'True': + --output_bam '$output_bam' + #end if + #if str($make_vcf) == 'True' + --output_vcf '$output_vcf' + #end if + #if str($make_gff) == 'True' + --output_gff '$output_gff' + #end if + #if str($make_fasta) == 'True' + --output_fasta '$output_fasta' + #end if --sample_name '$sample_name' --threads "\${GALAXY_SLOTS:-2}" ]]></command> <inputs> <param name="input" argument="--input_path" type="data" format="bam"/> - + <param name="make_bam" type="boolean" truevalue="True" checked="true" label="Produce an alignment file containing evidence for insertions."/> + <param name="make_vcf" type="boolean" truevalue="True" checked="true" label="Produce a VCF file describing the insertion that have been found."/> + <param name="make_gff" type="boolean" truevalue="True" checked="false" label="produce a GFF file describing the insertions that have been found."/> + <param name="make_fasta" type="boolean" checked="True" truevalue="True" label="Produce a fasta file containing assembled contigs."/> <expand macro="reference_source_conditional" reference_type="transposon"/> <expand macro="reference_source_conditional" reference_type="genome"/> </inputs> <outputs> - <data name="output_bam" format="bam" label="findcluster BAM on $on_string"/> - <data name="output_fasta" format="fasta" label="findcluster contigs on $on_string"/> - <data name="output_gff" format="gff3" label="findcluster GFF on $on_string"/> + <data name="output_bam" format="bam" label="findcluster BAM on $on_string"> + <filter>make_bam</filter> + </data> + <data name="output_fasta" format="fasta" label="findcluster contigs on $on_string"> + <filter>make_fasta</filter> + </data> + <data name="output_vcf" format="vcf" label="findcluster VCF on $on_string"> + <filter>make_vcf</filter> + </data> + <data name="output_gff" format="gff3" label="findcluster GFF on $on_string"> + <filter>make_gff</filter> + </data> </outputs> <tests> <test> <param name="input" value="extended_and_annotated_roi.bam" ftype="bam"/> + <param name="make_gff" value="true"/> <output name="output_bam" file="three_cluster_out.bam" ftype="bam" lines_diff="2"/> - <output name="output_gff" file="three_cluster_out.gff" ftype="gff3" lines_diff="0"/> + <output name="output_gff" file="three_cluster_out.gff" ftype="gff3" compare="sim_size"/> + <output name="output_vcf" file="three_cluster_out.vcf" ftype="vcf" compare="sim_size"/> </test> <test> <param name="input" value="extended_and_annotated_roi.bam" ftype="bam"/> <param name="transposon_source|reference_source_selector" value="history"/> <param name="transposon_source|ref_file" value="reference.fasta" ftype="fasta"/> + <param name="make_gff" value="true"/> <output name="output_bam" file="three_cluster_out.bam" ftype="bam" lines_diff="2"/> <output name="output_gff"> <assert_contents>
--- a/plot_coverage.xml Sat Jan 06 07:21:07 2018 -0500 +++ b/plot_coverage.xml Sun Apr 01 03:38:16 2018 -0400 @@ -1,7 +1,7 @@ -<tool id="plot_coverage" name="Plot coverage" version="0.4.6"> +<tool id="plot_coverage" name="Plot coverage" version="0.4.10"> <description>as area plot between BAM files</description> <requirements> - <requirement type="package" version="0.4.6">readtagger</requirement> + <requirement type="package" version="0.4.10">readtagger</requirement> </requirements> <command detect_errors="aggressive"><