annotate htseq-count.xml @ 2:f5d08224af89

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