0
|
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: “<chrom>:<start>-<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 <expand macro="citations" />
|
|
45 </tool>
|