Mercurial > repos > iuc > bedtools
diff getfastaBed.xml @ 1:82aac94b06c3 draft
Uploaded
author | iuc |
---|---|
date | Thu, 08 Jan 2015 14:25:51 -0500 |
parents | b8348686a0b9 |
children | 607c0576c6ab |
line wrap: on
line diff
--- a/getfastaBed.xml Tue Nov 04 01:45:04 2014 -0500 +++ b/getfastaBed.xml Thu Jan 08 14:25:51 2015 -0500 @@ -6,30 +6,59 @@ <expand macro="requirements" /> <expand macro="stdio" /> <command> +<![CDATA[ bedtools getfasta $name $tab $strand $split -fi $fasta - -bed $inputA + -bed $input -fo $output +]]> </command> <inputs> - <param format="bed,vcf,gff,gff3" name="inputA" type="data" label="BED/VCF/GFF file" /> + <param format="bed,vcf,gff,gff3" name="input" type="data" label="BED/VCF/GFF file" /> <param format="fasta" name="fasta" type="data" label="Fasta file" /> - <param name="name" type="boolean" checked="false" truevalue="-name" falsevalue="" label="Use the “name” column in the BED file for the FASTA headers in the output FASTA file" /> - <param name="tab" type="boolean" checked="false" truevalue="-tab" falsevalue="" label="Report extract sequences in a tab-delimited format instead of in FASTA format" /> - <param name="strand" type="boolean" checked="false" truevalue="-s" falsevalue="" label="Force strandedness" help="If the feature occupies the antisense strand, the sequence will be reverse complemented." /> - <param name="split" type="boolean" checked="false" truevalue="-split" falsevalue="" label="Given BED12 input, extract and concatenate the sequences from the BED 'blocks' (e.g., exons)" /> + <param name="name" type="boolean" checked="false" truevalue="-name" falsevalue="" + label="Use the 'name' column in the BED file for the FASTA headers in the output FASTA file" + help="(-name)" /> + <param name="tab" type="boolean" checked="false" truevalue="-tab" falsevalue="" + label="Report extract sequences in a tab-delimited format instead of in FASTA format" + help="(-tab)" /> + <param name="strand" type="boolean" checked="false" truevalue="-s" falsevalue="" + label="Force strandedness" + help="If the feature occupies the antisense strand, the sequence will be reverse complemented. (-s)" /> + <expand macro="split" /> </inputs> <outputs> - <data format="fasta" name="output" /> + <data format="fasta" name="output"> + <change_format> + <when input="tab" value="-tab" format="tabular" /> + </change_format> + </data> </outputs> + <tests> + <test> + <param name="input" value="nucBed1.bed" ftype="bed" /> + <param name="fasta" value="nucBed1.fasta" ftype="fasta" /> + <param name="tab" value="False" /> + <param name="split" value="False" /> + <output name="output" file="getfastaBed_result1.bed" ftype="fasta" /> + </test> + <test> + <param name="input" value="nucBed1.bed" ftype="bed" /> + <param name="fasta" value="nucBed1.fasta" ftype="fasta" /> + <param name="tab" value="True" /> + <param name="split" value="False" /> + <output name="output" file="getfastaBed_result2.tabular" ftype="tabular" /> + </test> + </tests> <help> +<![CDATA[ **What it does** -bedtools getfasta will extract the sequence defined by the coordinates in a BED interval and create a new FASTA entry in the output file for each extracted sequence. By default, the FASTA header for each extracted sequence will be formatted as follows: “<chrom>:<start>-<end>”. +bedtools getfasta will extract the sequence defined by the coordinates in a BED interval and create a new FASTA entry in the output file for each extracted sequence. By default, the FASTA header for each extracted sequence will be formatted as follows: “>chrom>:<start>-<end>”. .. image:: $PATH_TO_IMAGES/getfasta-glyph.png @@ -40,6 +69,7 @@ 2. You can use the UNIX fold command to set the line width of the FASTA output. For example, fold -w 60 will make each line of the FASTA file have at most 60 nucleotides for easy viewing. @REFERENCES@ +]]> </help> <expand macro="citations" /> </tool>