Mercurial > repos > iuc > bedtools

comparison getfastaBed.xml @ 0:4fb5ea02b441 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author iuc
date Tue, 15 Jul 2014 14:42:23 -0400
parents
children 7511823bdea1
comparison
equal deleted inserted replaced
-1:000000000000 0:4fb5ea02b441
1 <tool id="bedtools_getfastabed" name="GetFastaBed" version="@WRAPPER_VERSION@.0">
2 <description></description>
3 <macros>
4 <import>macros.xml</import>
5 </macros>
6 <expand macro="requirements" />
7 <expand macro="stdio" />
8 <command>
9 bedtools getfasta
10 $name
11 $tab
12 $strand
13 $split
14 -fi $fasta
15 -bed $inputA
16 -fo $output
17 </command>
18 <inputs>
19 <param format="bed,vcf,gff,gff3" name="inputA" type="data" label="BED/VCF/GFF file" />
20 <param format="fasta" name="fasta" type="data" label="Fasta file" />
21 <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" />
22 <param name="tab" type="boolean" checked="false" truevalue="-tab" falsevalue="" label="Report extract sequences in a tab-delimited format instead of in FASTA format" />
23 <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." />
24 <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)" />
25 </inputs>
26 <outputs>
27 <data format="fasta" name="output" />
28 </outputs>
29 <help>
30 **What it does**
31
32 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: “&lt;chrom>:&lt;start>-&lt;end>”.
33
34 .. image:: $PATH_TO_IMAGES/getfasta-glyph.png
35
36 .. class:: warningmark
37
38 1. The headers in the input FASTA file must exactly match the chromosome column in the BED file.
39
40 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.
41
42 @REFERENCES@
43 </help>
44 </tool>