Mercurial > repos > lparsons > htseq_count
annotate htseq-count.xml @ 3:c6d73933f601
Fixed issues with tool shed parsing
| author | Lance Parsons <lparsons@princeton.edu> |
|---|---|
| date | Fri, 31 Aug 2012 16:50:42 -0400 |
| parents | 8f0f896ec527 |
| children | 359d40333595 |
| rev | line source |
|---|---|
| 0 | 1 <tool id="htseq_count" name="htseq-count" version="0.1"> |
| 2 <description> - Count aligned reads in a BAM file that overlap features in a GFF file</description> | |
| 3 <version_command>htseq-count -h | grep version | sed 's/^\(.*\)*\(version .*\)\./\2/'</version_command> | |
| 4 <requirements> | |
|
2
8f0f896ec527
Added tool_dependencies.xml
Lance Parsons <lparsons@princeton.edu>
parents:
1
diff
changeset
|
5 <requirement type="package" version="1.6.2">numpy</requirement> |
|
8f0f896ec527
Added tool_dependencies.xml
Lance Parsons <lparsons@princeton.edu>
parents:
1
diff
changeset
|
6 <requirement type="package" version="0.5.3.p9">htseq</requirement> |
|
1
4de3f044aaeb
Added version to requirements
Lance Parsons <lparsons@princeton.edu>
parents:
0
diff
changeset
|
7 <requirement type="package" version="0.1.18">samtools</requirement> |
| 0 | 8 </requirements> |
| 9 <command> | |
| 10 ##set up input files | |
| 11 #set $reference_fasta_filename = "localref.fa" | |
| 12 #if $samout_conditional.samout: | |
| 13 #if str( $samout_conditional.reference_source.reference_source_selector ) == "history": | |
| 14 ln -s "${samout_conditional.reference_source.ref_file}" "${reference_fasta_filename}" && | |
| 15 samtools faidx "${reference_fasta_filename}" 2>&1 || echo "Error running samtools faidx for htseq-count" >&2 && | |
| 16 #else: | |
| 17 #set $reference_fasta_filename = str( $samout_conditional.reference_source.ref_file.fields.path ) | |
| 18 #end if | |
| 19 #end if | |
| 20 | |
| 21 #if $samfile.extension == "bam": | |
| 22 samtools view $samfile | | |
| 23 #end if | |
| 24 htseq-count | |
| 25 --mode=$mode | |
| 26 --stranded=$stranded | |
| 27 --minaqual=$minaqual | |
| 28 --type=$type | |
| 29 --idattr=$idattr | |
| 30 #if $samout_conditional.samout: | |
| 31 --samout=$__new_file_path__/${samoutfile.id}_tmp | |
| 32 #end if | |
| 33 #if $samfile.extension == "bam": | |
| 34 - | |
| 35 #else | |
| 36 $samfile | |
| 37 #end if | |
| 38 $gfffile | |
| 39 > $counts | |
| 40 #if $samout_conditional.samout: | |
| 41 && samtools view -Su -t ${reference_fasta_filename}.fai $__new_file_path__/${samoutfile.id}_tmp | samtools sort -o - sorted > $samoutfile | |
| 42 #end if</command> | |
| 43 <inputs> | |
| 44 <param format="sam, bam" name="samfile" type="data" label="Aligned SAM File"> | |
| 45 <help>Paired-End data must be sorted by QUERY NAME, use Picard Read Mate Fixer and Query name sort order before using this tool on paired data</help> | |
| 46 </param> | |
| 47 <param format="gff" name="gfffile" type="data" label="GFF File"/> | |
| 48 <param name="mode" type="select" label="Mode"> | |
| 49 <help>Mode to handle reads overlapping more than one feature.</help> | |
| 50 <option value="union" selected="true">Union</option> | |
| 51 <option value="intersection-strict">Intersection (strict)</option> | |
| 52 <option value="intersection-nonempty">Intersection (nonempty)</option> | |
| 53 </param> | |
| 54 <param name="stranded" type="select" label="Stranded"> | |
| 55 <help>Specify whether the data is from a strand-specific assay. 'Reverse' means yes with reversed strand interpretation.</help> | |
| 56 <option value="yes" selected="true">Yes</option> | |
| 57 <option value="no">No</option> | |
| 58 <option value="reverse">Reverse</option> | |
| 59 </param> | |
| 60 <param name="minaqual" type="integer" value="0" label="Minimum alignment quality"> | |
| 61 <help>Skip all reads with alignment quality lower than the given minimum value</help> | |
| 62 </param> | |
| 63 <param name="type" type="text" value="exon" label="Feature type"> | |
| 64 <help>Feature type (3rd column in GFF file) to be used. All features of other types are ignored. The default, suitable for RNA-Seq and Ensembl GTF files, is exon.</help> | |
| 65 </param> | |
| 66 <param name="idattr" type="text" value="gene_id" label="ID Attribute"> | |
| 67 <help>GFF attribute to be used as feature ID. Several GFF lines with the same feature ID will be considered as parts of the same feature. The feature ID is used to identity the counts in the output table. The default, suitable for RNA-SEq and Ensembl GTF files, is gene_id.</help> | |
| 68 </param> | |
| 69 <conditional name="samout_conditional"> | |
| 70 <param name="samout" type="boolean" value="False" truevalue="True" falsevalue="False" label="Additional BAM Output"> | |
| 71 <help>Write out all SAM alignment records into an output BAM file, annotating each line with its assignment to a feature or a special counter (as an optional field with tag ‘XF’).</help> | |
| 72 </param> | |
| 73 <when value="True"> | |
| 74 <conditional name="reference_source"> | |
| 75 <param name="reference_source_selector" type="select" label="Choose the source for the reference list"> | |
| 76 <option value="cached">Locally cached</option> | |
| 77 <option value="history">History</option> | |
| 78 </param> | |
| 79 <when value="cached"> | |
| 80 <param name="ref_file" type="select" label="Using reference genome"> | |
| 81 <options from_data_table="sam_fa_indexes"> | |
|
3
c6d73933f601
Fixed issues with tool shed parsing
Lance Parsons <lparsons@princeton.edu>
parents:
2
diff
changeset
|
82 <filter type="data_meta" key="dbkey" ref="samfile" column="3"/> |
| 0 | 83 </options> |
| 84 <validator type="no_options" message="A built-in reference genome is not available for the build associated with the selected input file"/> | |
| 85 </param> | |
| 86 </when> | |
| 87 <when value="history"> <!-- FIX ME!!!! --> | |
| 88 <param name="ref_file" type="data" format="fasta" label="Using reference file" /> | |
| 89 </when> | |
| 90 </conditional> | |
| 91 </when> | |
| 92 </conditional> | |
| 93 </inputs> | |
| 94 | |
| 95 <outputs> | |
| 96 <data format="tabular" name="counts" label="${tool.name} on ${on_string}"/> | |
| 97 <data format="bam" name="samoutfile" label="${tool.name} on ${on_string} (BAM)"> | |
| 98 <filter>samout_conditional['samout']</filter> | |
| 99 </data> | |
| 100 </outputs> | |
| 101 | |
| 102 <stdio> | |
| 103 <exit_code range="1:" level="fatal" description="Unknown error occurred" /> | |
| 104 </stdio> | |
| 105 | |
| 106 <tests> | |
| 107 <test> | |
| 108 <param name="samfile" value="htseq-test.sam" /> | |
| 109 <param name="gfffile" value="htseq-test.gff" /> | |
| 110 <param name="samout" value="False" /> | |
| 111 <output name="counts" file="htseq-test_counts.tsv" /> | |
| 112 </test> | |
| 113 <test> | |
| 114 <param name="samfile" value="htseq-test.bam" /> | |
| 115 <param name="gfffile" value="htseq-test.gff" /> | |
| 116 <param name="samout" value="False" /> | |
| 117 <output name="counts" file="htseq-test_counts.tsv" /> | |
| 118 </test> | |
| 119 <!-- Seems to be an issue setting the $reference_fasta_filename variable during test | |
| 120 <test> | |
| 121 <param name="samfile" value="htseq-test.sam" /> | |
| 122 <param name="gfffile" value="htseq-test.gff" /> | |
| 123 <param name="samout" value="True" /> | |
| 124 <param name="reference_source_selector" value="history" /> | |
| 125 <param name="ref_file" value="htseq-test_reference.fasta" /> | |
| 126 <output name="counts" file="htseq-test_counts.tsv" /> | |
| 127 <output name="samoutfile" file="htseq-test_samout.bam" /> | |
| 128 </test> | |
| 129 --> | |
| 130 | |
| 131 </tests> | |
| 132 | |
| 133 <help> | |
| 134 Overview | |
| 135 -------- | |
| 136 | |
| 137 This tool takes an alignment file in SAM or BAM format and feature file in GFF format | |
| 138 and calculates the number of reads mapping to each feature. It uses the *htseq-count* | |
| 139 script that is part of the HTSeq python module. See | |
| 140 http://www-huber.embl.de/users/anders/HTSeq/doc/count.html for details. | |
| 141 | |
| 142 A feature is an interval (i.e., a range of positions) on a chromosome or a union of | |
| 143 such intervals. In the case of RNA-Seq, the features are typically genes, where | |
| 144 each gene is considered here as the union of all its exons. One may also consider | |
| 145 each exon as a feature, e.g., in order to check for alternative splicing. For | |
| 146 comparative ChIP-Seq, the features might be binding regions from a pre-determined | |
| 147 list. | |
| 148 | |
| 149 **Paired-end Data MUST be sorted by QUERY NAME first** | |
| 150 | |
| 151 This tool requires that paired-end data be sorted by query name, which is NOT the default for Galaxy. Using the Picard Paired Read Mate Fixer with Query name sort FIRST is required for paired end data. | |
| 152 | |
| 153 | |
| 154 Overlap Modes | |
| 155 ------------- | |
| 156 | |
| 157 Special care must be taken to decide how to deal with reads that overlap more than one feature. | |
| 158 | |
| 159 The htseq-count script allows to choose between three modes: *union*, *intersection-strict*, and *intersection-nonempty*. | |
| 160 | |
| 161 The following figure illustrates the effect of these three modes: | |
| 162 | |
| 163 .. image:: /static/images/count_modes.png | |
| 164 :width: 500 | |
| 165 | |
| 166 Strandedness | |
| 167 ------------ | |
| 168 | |
| 169 **Important**: The default for strandedness is yes. If your RNA-Seq data has not been made with a strand-specific protocol, this causes half of the reads to be lost. Hence, make sure to set the option Stranded to 'No' unless you have strand-specific data! | |
| 170 | |
| 171 Output | |
| 172 ------ | |
| 173 | |
| 174 The script outputs a table with counts for each feature, followed by the special counters, which count reads that were not counted for any feature for various reasons, namely | |
| 175 | |
| 176 - *no_feature*: reads which could not be assigned to any feature (set S as described above was empty). | |
| 177 | |
| 178 - *ambiguous*: reads which could have been assigned to more than one feature and hence were not counted for any of these (set S had mroe than one element). | |
| 179 | |
| 180 - *too_low_aQual*: reads which were not counted due to the -a option, see below | |
| 181 | |
| 182 - *not_aligned*: reads in the SAM file without alignment | |
| 183 | |
| 184 - *alignment_not_unique*: reads with more than one reported alignment. These reads are recognized from the NH optional SAM field tag. (If the aligner does not set this field, multiply aligned reads will be counted multiple times.) | |
| 185 | |
| 186 | |
| 187 Options Summary | |
| 188 --------------- | |
| 189 | |
| 190 Usage: htseq-count [options] sam_file gff_file | |
| 191 | |
| 192 This script takes an alignment file in SAM format and a feature file in GFF | |
| 193 format and calculates for each feature the number of reads mapping to it. See | |
| 194 http://www-huber.embl.de/users/anders/HTSeq/doc/count.html for details. | |
| 195 | |
| 196 Options: | |
| 197 -h, --help show this help message and exit | |
| 198 -m MODE, --mode=MODE mode to handle reads overlapping more than one | |
| 199 feature(choices: union, intersection-strict, | |
| 200 intersection-nonempty; default: union) | |
| 201 -s STRANDED, --stranded=STRANDED | |
| 202 whether the data is from a strand-specific assay. | |
| 203 Specify 'yes', 'no', or 'reverse' (default: yes). | |
| 204 'reverse' means 'yes' with reversed strand | |
| 205 interpretation | |
| 206 -a MINAQUAL, --minaqual=MINAQUAL | |
| 207 skip all reads with alignment quality lower than the | |
| 208 given minimum value (default: 0) | |
| 209 -t FEATURETYPE, --type=FEATURETYPE | |
| 210 feature type (3rd column in GFF file) to be used, all | |
| 211 features of other type are ignored (default, suitable | |
| 212 for Ensembl GTF files: exon) | |
| 213 -i IDATTR, --idattr=IDATTR | |
| 214 GFF attribute to be used as feature ID (default, | |
| 215 suitable for Ensembl GTF files: gene_id) | |
| 216 -o SAMOUT, --samout=SAMOUT | |
| 217 write out all SAM alignment records into an output SAM | |
| 218 file called SAMOUT, annotating each line with its | |
| 219 feature assignment (as an optional field with tag | |
| 220 'XF') | |
| 221 -q, --quiet suppress progress report and warnings | |
| 222 | |
| 223 Written by Simon Anders (sanders@fs.tum.de), European Molecular Biology | |
| 224 Laboratory (EMBL). (c) 2010. Released under the terms of the GNU General | |
| 225 Public License v3. Part of the 'HTSeq' framework. | |
| 226 </help> | |
| 227 | |
| 228 </tool> |