Mercurial > repos > lparsons > bam_to_bigwig
annotate bam_to_bigwig.xml @ 5:5b40b93ebae3 default tip
Added advanced settings
| author | Lance Parsons <lparsons@princeton.edu> |
|---|---|
| date | Wed, 24 Aug 2011 16:44:28 -0400 |
| parents | b8525303b945 |
| children |
| rev | line source |
|---|---|
| 0 | 1 <tool id="bam_to_bigwig" name="BAM to BigWig" version="0.1"> |
| 2 <description>Calculates coverage from a BAM alignment file</description> | |
| 3 <requirements> | |
| 4 <requirement type="binary">genomeCoverageBed</requirement> | |
| 5 <requirement type="binary">bedGraphToBigWig</requirement> | |
| 6 </requirements> | |
| 7 <command> | |
| 8 genomeCoverageBed -ibam $align -bg $split -g | |
| 9 `awk '($0 !~ /^#.*/) && ($2 == "$align.dbkey") {print $3}' ${GALAXY_DATA_INDEX_DIR}/sam_fa_indices.loc`.fai | |
| 10 > tmp.bedgraph; | |
|
5
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
11 bedGraphToBigWig |
|
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
12 #if $settings.settingsType == "full": |
|
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
13 -blockSize=${settings.blockSize} -itemsPerSlot=${settings.itemsPerSlot} ${settings.unc} |
|
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
14 #end if |
|
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
15 tmp.bedgraph |
| 0 | 16 `awk '($0 !~ /^#.*/) && ($2 == "$align.dbkey") {print $3}' ${GALAXY_DATA_INDEX_DIR}/sam_fa_indices.loc`.fai |
| 17 $out; | |
| 18 </command> | |
| 19 <inputs> | |
| 20 <param name="align" type="data" format="bam" label="Select the BAM file to generate the pileup file for"> | |
| 21 <validator type="unspecified_build" /> | |
| 22 <validator type="dataset_metadata_in_file" filename="sam_fa_indices.loc" metadata_name="dbkey" metadata_column="1" message="Sequences are not currently available for the specified build." line_startswith="index" /> | |
| 23 </param> | |
| 24 <param name="split" type="boolean" value="false" label="Separate split reads" help="Do not use gaps when computing coverage. This is recommended for RNA-Seq where reads are mapped across long splice junctions. Uses the CIGAR "N" and "D" operations to infer the blocks for computing coverage." truevalue="-split" falsevalue="" /> | |
|
5
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
25 <conditional name="settings"> |
|
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
26 <param name="settingsType" type="select" label="Converter settings to use" help="Default settings should usually be used."> |
|
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
27 <option value="preset">Default</option> |
|
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
28 <option value="full">Full parameter list</option> |
|
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
29 </param> |
|
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
30 <when value="preset" /> |
|
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
31 <when value="full"> |
|
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
32 <param name="blockSize" size="4" type="integer" value="256" label="Items to bundle in r-tree" help="Default is 256 (blockSize)" /> |
|
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
33 <param name="itemsPerSlot" size="4" type="integer" value="1024" label="Data points bundled at lowest level" help="Default is 1024 (itemsPerSlot)" /> |
|
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
34 <param name="unc" type="boolean" truevalue="-unc" falsevalue="" checked="False" label="Do not use compression" help="(unc)"/> |
|
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
35 </when> |
|
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
36 </conditional> |
|
5b40b93ebae3
Added advanced settings
Lance Parsons <lparsons@princeton.edu>
parents:
3
diff
changeset
|
37 |
| 0 | 38 </inputs> |
| 39 <outputs> | |
| 40 <data format="bigwig" name="out" /> | |
| 41 </outputs> | |
| 42 | |
| 43 <help> | |
| 44 **What it does** | |
| 45 | |
| 46 Creates a coverage file in BigWig format, given a BAM alignment file. | |
| 47 | |
| 48 **Input** | |
| 49 | |
| 50 A BAM alignment file. This needs to have the genome database build used in alignment annotated. If your file has '?' for the database build, click on the pencil icon to edit the alignment attributes, and specify the organism used to align against. | |
| 51 | |
| 52 The BAM file must also be sorted. If you generated the BAM file in Galaxy, this will be true, however, if you uploaded the file, you must make sure it was sorted. | |
| 53 | |
| 54 **Output** | |
| 55 | |
| 56 BigWig files can be loaded directly from Galaxy into the UCSC browser. They can be loaded incrementally by UCSC, so a single file can be used to represent the entire genome without having to upload the entire thing as a custom track. | |
| 57 </help> | |
| 58 | |
| 59 </tool> |