Mercurial > repos > malex > bayesase
diff create_sam_file_with_features.xml @ 0:e979cb57a5d5 draft default tip
"planemo upload for repository https://github.com/McIntyre-Lab/BayesASE/tree/main/galaxy commit 9b70598ef46a73632d9e0fa0c6ce6776fb5e9d6a"
author | malex |
---|---|
date | Thu, 14 Jan 2021 21:51:36 +0000 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/create_sam_file_with_features.xml Thu Jan 14 21:51:36 2021 +0000 @@ -0,0 +1,69 @@ +<tool id="base_create_new_sam_file" name="Create new SAM file" version="21.1.13"> + <description> with features of interest in RNAME field</description> + <macros> + <import>macros.xml</import> + </macros> + <expand macro="requirements"/> + <command> +<![CDATA[ + awk -v OFS=' ' 'FNR==NR{a[$1,$2]=$7; next} {$3=a[$3,$4]}1' '$BED3' '$SAMFILE' | awk -F'\t' '$3!=""' > '$NEWSAM' +]]> + </command> + <inputs> + <param name="BED3" type="data" format="tabular" label="BED file generated from the Remove Reads That Do Not Overlap tool" /> + <param name="SAMFILE" type="data" format="sam" label="Input SAM file containing Uniquely Mapped Reads" + help=" SAM file to convert chromosome names to feature IDs in RNAME field" /> + </inputs> + <outputs> + <data name="NEWSAM" format="sam" label="${tool.name} on ${on_string}: New SAM file with features instead of chrom" /> + </outputs> + <tests> + <test> + <param name="BED3" value="align_and_counts_test_data/drop_nonintersecting_reads_for_BASE.tabular" ftype="tabular"/> + <param name="SAMFILE" value="align_and_counts_test_data/W1118_G1_unique_sam_for_BASE.sam" ftype="sam"/> + <output name="NEWSAM" file="align_and_counts_test_data/W1118_G1_create_new_SAM_file_with_features_BASE_test_data.sam" ftype="sam"/> + </test> + </tests> + <help><![CDATA[ + **Tool Description** + +This tool reformats the SAM file containing uniquely mapping reads using the multi-column BED file created by the *Remove reads that do not overlap with a feature for BASE* tool. +The resulting SAM file contains the name of genic features in the 'chrom' field. + + +---------------------------------------------------------------------------------------------- + +**Inputs** + + -The tool has two required inputs: + +**BED File** + The tool requires a BED file of reads overlapping genic features of interest. + The *Remove reads that do not overlap with a feature for BayesASE* tool creates this BED file. + +Example input BED file + + +---------------+---------+---------+---------+---------+------------+------------+ + | Chrom | Start | End | Name | Score | Strand |ThickStart | + +===============+=========+=========+=========+=========+============+============+ + | X | 2190 | 2245 | X | 1 | 2300 |l(1)G0196 | + +---------------+---------+---------+---------+---------+------------+------------+ + | 2R | 1502 | 1834 | 2R | 50 |19000 |Mapmodulin | + +---------------+---------+---------+---------+---------+------------+------------+ + +**Unique SAM File** + Input the SAM file containing alignments for all uniquely mapping reads. This file can be created by the **BWASplitSAM** tool. + +**Outputs** + The tool outputs a new SAM file with the genic feature (ThickStart column in BED file) in the RNAME field instead of the chromosome name. + + ]]></help> + <citations> + <citation type="bibtex">@ARTICLE{Miller20BASE, + author = {Brecca Miller, Alison M. Morse, Elyse Borgert, Zihao Liu, Kelsey Sinclair, Gavin Gamble, Fei Zou, Jeremy Newman, Luis Leon Novello, Fabio Marroni, Lauren M. McIntyre}, + title = {Testcrosses are an efficient strategy for identifying cis regulatory variation: Bayesian analysis of allele imbalance among conditions (BASE)}, + journal = {????}, + year = {submitted for publication} + }</citation> + </citations> +</tool>