Repository 'htseq_count'
hg clone https://toolshed.g2.bx.psu.edu/repos/lparsons/htseq_count

Changeset 0:3fdeebd7e710 (2012-09-04)
Next changeset 1:14e18dc9ed13 (2012-09-04)
Commit message:
Initial commit
added:
htseq-count.xml
sam_fa_indices.loc.sample
static/images/count_modes.png
test-data/htseq-test.bam
test-data/htseq-test.gff
test-data/htseq-test.sam
test-data/htseq-test_counts.tsv
tool_data_table_conf.xml.sample
tool_dependencies.xml
b
diff -r 000000000000 -r 3fdeebd7e710 htseq-count.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/htseq-count.xml Tue Sep 04 13:45:36 2012 -0400
[
b'@@ -0,0 +1,225 @@\n+<tool id="htseq_count" name="htseq-count" version="0.1">\n+    <description> - Count aligned reads in a BAM file that overlap features in a GFF file</description>\n+    <version_command>htseq-count -h | grep version | sed \'s/^\\(.*\\)*\\(version .*\\)\\./\\2/\'</version_command>\n+    <requirements>\n+        <requirement type="package" version="0.5.3p9">htseq</requirement>\n+        <requirement type="package" version="0.1.18">samtools</requirement>\n+    </requirements>\n+    <command>\n+    ##set up input files\n+    #set $reference_fasta_filename = "localref.fa"\n+    #if $samout_conditional.samout:\n+        #if str( $samout_conditional.reference_source.reference_source_selector ) == "history":\n+            ln -s "${samout_conditional.reference_source.ref_file}" "${reference_fasta_filename}" &amp;&amp;\n+            samtools faidx "${reference_fasta_filename}" 2&gt;&amp;1 || echo "Error running samtools faidx for htseq-count" &gt;&amp;2 &amp;&amp;\n+        #else:\n+            #set $reference_fasta_filename = str( $samout_conditional.reference_source.ref_file.fields.path )\n+        #end if\n+    #end if\n+\n+    #if $samfile.extension == "bam":\n+        samtools view $samfile | \n+    #end if\n+    htseq-count \n+    --mode=$mode \n+    --stranded=$stranded \n+    --minaqual=$minaqual \n+    --type=$type \n+    --idattr=$idattr \n+    #if $samout_conditional.samout:\n+        --samout=$__new_file_path__/${samoutfile.id}_tmp\n+    #end if\n+    #if $samfile.extension == "bam":\n+        - \n+    #else\n+        $samfile \n+    #end if\n+    $gfffile \n+    &gt; $counts\n+    #if $samout_conditional.samout:\n+        &amp;&amp; samtools view -Su -t ${reference_fasta_filename}.fai $__new_file_path__/${samoutfile.id}_tmp | samtools sort -o - sorted > $samoutfile\n+    #end if</command>\n+    <inputs>\n+        <param format="sam, bam" name="samfile" type="data" label="Aligned SAM File">\n+            <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>\n+        </param>\n+        <param format="gff" name="gfffile" type="data" label="GFF File"/>\n+        <param name="mode" type="select" label="Mode">\n+            <help>Mode to handle reads overlapping more than one feature.</help>\n+            <option value="union" selected="true">Union</option>\n+            <option value="intersection-strict">Intersection (strict)</option>\n+            <option value="intersection-nonempty">Intersection (nonempty)</option>\n+        </param>\n+        <param name="stranded" type="select" label="Stranded">\n+            <help>Specify whether the data is from a strand-specific assay. \'Reverse\' means yes with reversed strand interpretation.</help>\n+            <option value="yes" selected="true">Yes</option>\n+            <option value="no">No</option>\n+            <option value="reverse">Reverse</option>\n+        </param>\n+        <param name="minaqual" type="integer" value="0" label="Minimum alignment quality">\n+            <help>Skip all reads with alignment quality lower than the given minimum value</help>\n+        </param>\n+        <param name="type" type="text" value="exon" label="Feature type">\n+            <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>\n+        </param>\n+        <param name="idattr" type="text" value="gene_id" label="ID Attribute">\n+            <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>\n+        </param>\n+        <conditional name="samout_conditional">\n+            <param name="samout" type="boolean" value="False" truevalue="True" falsevalue="False" label="Additional BAM Output">\n+                <help>'..b"on as a feature, e.g., in order to check for alternative splicing. For \n+comparative ChIP-Seq, the features might be binding regions from a pre-determined \n+list.\n+\n+**Paired-end Data MUST be sorted by QUERY NAME first**\n+\n+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.\n+\n+\n+Overlap Modes\n+-------------\n+\n+Special care must be taken to decide how to deal with reads that overlap more than one feature. \n+\n+The htseq-count script allows to choose between three modes: *union*, *intersection-strict*, and *intersection-nonempty*.\n+\n+The following figure illustrates the effect of these three modes:\n+\n+.. image:: /static/images/count_modes.png\n+    :width: 500\n+\n+Strandedness\n+------------\n+\n+**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!\n+\n+Output\n+------\n+\n+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\n+\n+- *no_feature*: reads which could not be assigned to any feature (set S as described above was empty).\n+\n+- *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).\n+\n+- *too_low_aQual*: reads which were not counted due to the -a option, see below\n+\n+- *not_aligned*: reads in the SAM file without alignment\n+\n+- *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.)\n+\n+\n+Options Summary\n+---------------\n+\n+Usage: htseq-count [options] sam_file gff_file\n+\n+This script takes an alignment file in SAM format and a feature file in GFF\n+format and calculates for each feature the number of reads mapping to it. See\n+http://www-huber.embl.de/users/anders/HTSeq/doc/count.html for details.\n+\n+Options:\n+  -h, --help            show this help message and exit\n+  -m MODE, --mode=MODE  mode to handle reads overlapping more than one\n+                        feature(choices: union, intersection-strict,\n+                        intersection-nonempty; default: union)\n+  -s STRANDED, --stranded=STRANDED\n+                        whether the data is from a strand-specific assay.\n+                        Specify 'yes', 'no', or 'reverse' (default: yes).\n+                        'reverse' means 'yes' with reversed strand\n+                        interpretation\n+  -a MINAQUAL, --minaqual=MINAQUAL\n+                        skip all reads with alignment quality lower than the\n+                        given minimum value (default: 0)\n+  -t FEATURETYPE, --type=FEATURETYPE\n+                        feature type (3rd column in GFF file) to be used, all\n+                        features of other type are ignored (default, suitable\n+                        for Ensembl GTF files: exon)\n+  -i IDATTR, --idattr=IDATTR\n+                        GFF attribute to be used as feature ID (default,\n+                        suitable for Ensembl GTF files: gene_id)\n+  -o SAMOUT, --samout=SAMOUT\n+                        write out all SAM alignment records into an output SAM\n+                        file called SAMOUT, annotating each line with its\n+                        feature assignment (as an optional field with tag\n+                        'XF')\n+  -q, --quiet           suppress progress report and warnings\n+\n+Written by Simon Anders (sanders@fs.tum.de), European Molecular Biology\n+Laboratory (EMBL). (c) 2010. Released under the terms of the GNU General\n+Public License v3. Part of the 'HTSeq' framework.\n+    </help>\n+</tool>\n"
b
diff -r 000000000000 -r 3fdeebd7e710 sam_fa_indices.loc.sample
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/sam_fa_indices.loc.sample Tue Sep 04 13:45:36 2012 -0400
b
@@ -0,0 +1,28 @@
+#This is a sample file distributed with Galaxy that enables tools
+#to use a directory of Samtools indexed sequences data files.  You will need
+#to create these data files and then create a sam_fa_indices.loc file 
+#similar to this one (store it in this directory) that points to 
+#the directories in which those files are stored. The sam_fa_indices.loc 
+#file has this format (white space characters are TAB characters):
+#
+#index <seq> <location>
+#
+#So, for example, if you had hg18 indexed stored in 
+#/depot/data2/galaxy/sam/, 
+#then the sam_fa_indices.loc entry would look like this:
+#
+#index hg18 /depot/data2/galaxy/sam/hg18.fa
+#
+#and your /depot/data2/galaxy/sam/ directory
+#would contain hg18.fa and hg18.fa.fai files:
+#
+#-rw-r--r--  1 james    universe 830134 2005-09-13 10:12 hg18.fa
+#-rw-r--r--  1 james    universe 527388 2005-09-13 10:12 hg18.fa.fai
+#
+#Your sam_fa_indices.loc file should include an entry per line for 
+#each index set you have stored.  The file in the path does actually
+#exist, but it should never be directly used. Instead, the name serves
+#as a prefix for the index file.  For example:
+#
+#index hg18 /depot/data2/galaxy/sam/hg18.fa
+#index hg19 /depot/data2/galaxy/sam/hg19.fa
b
diff -r 000000000000 -r 3fdeebd7e710 static/images/count_modes.png
b
Binary file static/images/count_modes.png has changed
b
diff -r 000000000000 -r 3fdeebd7e710 test-data/htseq-test.bam
b
Binary file test-data/htseq-test.bam has changed
b
diff -r 000000000000 -r 3fdeebd7e710 test-data/htseq-test.gff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/htseq-test.gff Tue Sep 04 13:45:36 2012 -0400
b
b'@@ -0,0 +1,492 @@\n+I\tprotein_coding\texon\t335\t649\t.\t+\t.\t gene_id "YAL069W"; transcript_id "YAL069W"; exon_number "1"; gene_name "YAL069W"; transcript_name "YAL069W";\n+I\tprotein_coding\tCDS\t335\t646\t.\t+\t0\t gene_id "YAL069W"; transcript_id "YAL069W"; exon_number "1"; gene_name "YAL069W"; transcript_name "YAL069W"; protein_id "YAL069W";\n+I\tprotein_coding\tstart_codon\t335\t337\t.\t+\t0\t gene_id "YAL069W"; transcript_id "YAL069W"; exon_number "1"; gene_name "YAL069W"; transcript_name "YAL069W";\n+I\tprotein_coding\tstop_codon\t647\t649\t.\t+\t0\t gene_id "YAL069W"; transcript_id "YAL069W"; exon_number "1"; gene_name "YAL069W"; transcript_name "YAL069W";\n+I\tprotein_coding\texon\t538\t792\t.\t+\t.\t gene_id "YAL068W-A"; transcript_id "YAL068W-A"; exon_number "1"; gene_name "YAL068W-A"; transcript_name "YAL068W-A";\n+I\tprotein_coding\tCDS\t538\t789\t.\t+\t0\t gene_id "YAL068W-A"; transcript_id "YAL068W-A"; exon_number "1"; gene_name "YAL068W-A"; transcript_name "YAL068W-A"; protein_id "YAL068W-A";\n+I\tprotein_coding\tstart_codon\t538\t540\t.\t+\t0\t gene_id "YAL068W-A"; transcript_id "YAL068W-A"; exon_number "1"; gene_name "YAL068W-A"; transcript_name "YAL068W-A";\n+I\tprotein_coding\tstop_codon\t790\t792\t.\t+\t0\t gene_id "YAL068W-A"; transcript_id "YAL068W-A"; exon_number "1"; gene_name "YAL068W-A"; transcript_name "YAL068W-A";\n+I\tprotein_coding\texon\t1807\t2169\t.\t-\t.\t gene_id "YAL068C"; transcript_id "YAL068C"; exon_number "1"; gene_name "PAU8"; transcript_name "PAU8";\n+I\tprotein_coding\tCDS\t1810\t2169\t.\t-\t0\t gene_id "YAL068C"; transcript_id "YAL068C"; exon_number "1"; gene_name "PAU8"; transcript_name "PAU8"; protein_id "YAL068C";\n+I\tprotein_coding\tstart_codon\t2167\t2169\t.\t-\t0\t gene_id "YAL068C"; transcript_id "YAL068C"; exon_number "1"; gene_name "PAU8"; transcript_name "PAU8";\n+I\tprotein_coding\tstop_codon\t1807\t1809\t.\t-\t0\t gene_id "YAL068C"; transcript_id "YAL068C"; exon_number "1"; gene_name "PAU8"; transcript_name "PAU8";\n+I\tprotein_coding\texon\t2480\t2707\t.\t+\t.\t gene_id "YAL067W-A"; transcript_id "YAL067W-A"; exon_number "1"; gene_name "YAL067W-A"; transcript_name "YAL067W-A";\n+I\tprotein_coding\tCDS\t2480\t2704\t.\t+\t0\t gene_id "YAL067W-A"; transcript_id "YAL067W-A"; exon_number "1"; gene_name "YAL067W-A"; transcript_name "YAL067W-A"; protein_id "YAL067W-A";\n+I\tprotein_coding\tstart_codon\t2480\t2482\t.\t+\t0\t gene_id "YAL067W-A"; transcript_id "YAL067W-A"; exon_number "1"; gene_name "YAL067W-A"; transcript_name "YAL067W-A";\n+I\tprotein_coding\tstop_codon\t2705\t2707\t.\t+\t0\t gene_id "YAL067W-A"; transcript_id "YAL067W-A"; exon_number "1"; gene_name "YAL067W-A"; transcript_name "YAL067W-A";\n+I\tprotein_coding\texon\t7236\t9017\t.\t-\t.\t gene_id "YAL067C"; transcript_id "YAL067C"; exon_number "1"; gene_name "SEO1"; transcript_name "SEO1";\n+I\tprotein_coding\tCDS\t7239\t9017\t.\t-\t0\t gene_id "YAL067C"; transcript_id "YAL067C"; exon_number "1"; gene_name "SEO1"; transcript_name "SEO1"; protein_id "YAL067C";\n+I\tprotein_coding\tstart_codon\t9015\t9017\t.\t-\t0\t gene_id "YAL067C"; transcript_id "YAL067C"; exon_number "1"; gene_name "SEO1"; transcript_name "SEO1";\n+I\tprotein_coding\tstop_codon\t7236\t7238\t.\t-\t0\t gene_id "YAL067C"; transcript_id "YAL067C"; exon_number "1"; gene_name "SEO1"; transcript_name "SEO1";\n+I\tprotein_coding\texon\t10092\t10400\t.\t+\t.\t gene_id "YAL066W"; transcript_id "YAL066W"; exon_number "1"; gene_name "YAL066W"; transcript_name "YAL066W";\n+I\tprotein_coding\tCDS\t10092\t10397\t.\t+\t0\t gene_id "YAL066W"; transcript_id "YAL066W"; exon_number "1"; gene_name "YAL066W"; transcript_name "YAL066W"; protein_id "YAL066W";\n+I\tprotein_coding\tstart_codon\t10092\t10094\t.\t+\t0\t gene_id "YAL066W"; transcript_id "YAL066W"; exon_number "1"; gene_name "YAL066W"; transcript_name "YAL066W";\n+I\tprotein_coding\tstop_codon\t10398\t10400\t.\t+\t0\t gene_id "YAL066W"; transcript_id "YAL066W"; exon_number "1"; gene_name "YAL066W"; transcript_name "YAL066W";\n+I\tprotein_coding\texon\t11566\t11952\t.\t-\t.\t gene_id "YAL065C"; transcript_id "YAL065C"; exon_number "1"; gene_name "YAL065C"; transcript_name "YAL065C";\n+I\tprotein_coding\tCDS\t11569\t11952\t.\t-\t'..b'W";\n+I\tprotein_coding\tstop_codon\t221649\t221651\t.\t+\t0\t gene_id "YAR066W"; transcript_id "YAR066W"; exon_number "1"; gene_name "YAR066W"; transcript_name "YAR066W";\n+I\tprotein_coding\texon\t222397\t222882\t.\t+\t.\t gene_id "YAR068W"; transcript_id "YAR068W"; exon_number "1"; gene_name "YAR068W"; transcript_name "YAR068W";\n+I\tprotein_coding\tCDS\t222397\t222879\t.\t+\t0\t gene_id "YAR068W"; transcript_id "YAR068W"; exon_number "1"; gene_name "YAR068W"; transcript_name "YAR068W"; protein_id "YAR068W";\n+I\tprotein_coding\tstart_codon\t222397\t222399\t.\t+\t0\t gene_id "YAR068W"; transcript_id "YAR068W"; exon_number "1"; gene_name "YAR068W"; transcript_name "YAR068W";\n+I\tprotein_coding\tstop_codon\t222880\t222882\t.\t+\t0\t gene_id "YAR068W"; transcript_id "YAR068W"; exon_number "1"; gene_name "YAR068W"; transcript_name "YAR068W";\n+I\tprotein_coding\texon\t224002\t224295\t.\t-\t.\t gene_id "YAR069C"; transcript_id "YAR069C"; exon_number "1"; gene_name "YAR069C"; transcript_name "YAR069C";\n+I\tprotein_coding\tCDS\t224005\t224295\t.\t-\t0\t gene_id "YAR069C"; transcript_id "YAR069C"; exon_number "1"; gene_name "YAR069C"; transcript_name "YAR069C"; protein_id "YAR069C";\n+I\tprotein_coding\tstart_codon\t224293\t224295\t.\t-\t0\t gene_id "YAR069C"; transcript_id "YAR069C"; exon_number "1"; gene_name "YAR069C"; transcript_name "YAR069C";\n+I\tprotein_coding\tstop_codon\t224002\t224004\t.\t-\t0\t gene_id "YAR069C"; transcript_id "YAR069C"; exon_number "1"; gene_name "YAR069C"; transcript_name "YAR069C";\n+I\tprotein_coding\texon\t224554\t224853\t.\t-\t.\t gene_id "YAR070C"; transcript_id "YAR070C"; exon_number "1"; gene_name "YAR070C"; transcript_name "YAR070C";\n+I\tprotein_coding\tCDS\t224557\t224853\t.\t-\t0\t gene_id "YAR070C"; transcript_id "YAR070C"; exon_number "1"; gene_name "YAR070C"; transcript_name "YAR070C"; protein_id "YAR070C";\n+I\tprotein_coding\tstart_codon\t224851\t224853\t.\t-\t0\t gene_id "YAR070C"; transcript_id "YAR070C"; exon_number "1"; gene_name "YAR070C"; transcript_name "YAR070C";\n+I\tprotein_coding\tstop_codon\t224554\t224556\t.\t-\t0\t gene_id "YAR070C"; transcript_id "YAR070C"; exon_number "1"; gene_name "YAR070C"; transcript_name "YAR070C";\n+I\tprotein_coding\texon\t225451\t226854\t.\t+\t.\t gene_id "YAR071W"; transcript_id "YAR071W"; exon_number "1"; gene_name "PHO11"; transcript_name "PHO11";\n+I\tprotein_coding\tCDS\t225451\t226851\t.\t+\t0\t gene_id "YAR071W"; transcript_id "YAR071W"; exon_number "1"; gene_name "PHO11"; transcript_name "PHO11"; protein_id "YAR071W";\n+I\tprotein_coding\tstart_codon\t225451\t225453\t.\t+\t0\t gene_id "YAR071W"; transcript_id "YAR071W"; exon_number "1"; gene_name "PHO11"; transcript_name "PHO11";\n+I\tprotein_coding\tstop_codon\t226852\t226854\t.\t+\t0\t gene_id "YAR071W"; transcript_id "YAR071W"; exon_number "1"; gene_name "PHO11"; transcript_name "PHO11";\n+I\tprotein_coding\texon\t227733\t228944\t.\t+\t.\t gene_id "YAR073W"; transcript_id "YAR073W"; exon_number "1"; gene_name "IMD1"; transcript_name "IMD1";\n+I\tprotein_coding\tCDS\t227733\t228941\t.\t+\t0\t gene_id "YAR073W"; transcript_id "YAR073W"; exon_number "1"; gene_name "IMD1"; transcript_name "IMD1"; protein_id "YAR073W";\n+I\tprotein_coding\tstart_codon\t227733\t227735\t.\t+\t0\t gene_id "YAR073W"; transcript_id "YAR073W"; exon_number "1"; gene_name "IMD1"; transcript_name "IMD1";\n+I\tprotein_coding\tstop_codon\t228942\t228944\t.\t+\t0\t gene_id "YAR073W"; transcript_id "YAR073W"; exon_number "1"; gene_name "IMD1"; transcript_name "IMD1";\n+I\tprotein_coding\texon\t228835\t229308\t.\t+\t.\t gene_id "YAR075W"; transcript_id "YAR075W"; exon_number "1"; gene_name "YAR075W"; transcript_name "YAR075W";\n+I\tprotein_coding\tCDS\t228835\t229305\t.\t+\t0\t gene_id "YAR075W"; transcript_id "YAR075W"; exon_number "1"; gene_name "YAR075W"; transcript_name "YAR075W"; protein_id "YAR075W";\n+I\tprotein_coding\tstart_codon\t228835\t228837\t.\t+\t0\t gene_id "YAR075W"; transcript_id "YAR075W"; exon_number "1"; gene_name "YAR075W"; transcript_name "YAR075W";\n+I\tprotein_coding\tstop_codon\t229306\t229308\t.\t+\t0\t gene_id "YAR075W"; transcript_id "YAR075W"; exon_number "1"; gene_name "YAR075W"; transcript_name "YAR075W";\n'
b
diff -r 000000000000 -r 3fdeebd7e710 test-data/htseq-test.sam
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/htseq-test.sam Tue Sep 04 13:45:36 2012 -0400
b
b'@@ -0,0 +1,1717 @@\n+@HD\tVN:1.0\tSO:coordinate\n+@SQ\tSN:2-micron\tLN:6318\n+@SQ\tSN:MT\tLN:85779\n+@SQ\tSN:I\tLN:230208\n+@SQ\tSN:VI\tLN:270148\n+@SQ\tSN:III\tLN:316617\n+@SQ\tSN:IX\tLN:439885\n+@SQ\tSN:VIII\tLN:562643\n+@SQ\tSN:V\tLN:576869\n+@SQ\tSN:XI\tLN:666454\n+@SQ\tSN:X\tLN:745745\n+@SQ\tSN:XIV\tLN:784333\n+@SQ\tSN:II\tLN:813178\n+@SQ\tSN:XIII\tLN:924429\n+@SQ\tSN:XVI\tLN:948062\n+@SQ\tSN:XII\tLN:1078175\n+@SQ\tSN:VII\tLN:1090946\n+@SQ\tSN:XV\tLN:1091289\n+@SQ\tSN:IV\tLN:1531918\n+HWI-EAS225:1:10:1473:444#0/1\t0\tI\t12855\t255\t36M\t*\t0\t0\tTTGGTGATTCAGAAGAGATAGAAGAGGAAGAAGAGA\tABA@9>2<@@=?529.::=:?99?::99=:99=:<7\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:10:1426:363#0/1\t0\tI\t34474\t255\t36M\t*\t0\t0\tACATATGATAACGGGCAGAGTCAATATTGAGGACGG\tBA675>?6<>>>?37>59:=8;03?38#########\tXA:i:1\tMD:Z:24C11\tNM:i:1\n+HWI-EAS225:1:10:1661:622#0/1\t0\tI\t34541\t255\t36M\t*\t0\t0\tGAGAAGGAGTCTACAATCAAGATTATTCTGACTCCA\tBBA@;8=:>;;>;999;;;8878;:::;8788999:\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:10:1296:245#0/1\t0\tI\t34561\t255\t36M\t*\t0\t0\tGATTATTCTGACTCCAAACAATCACGGAGAGTTGAA\tB?@A>AA???;;??>>;;;;9;999>>9:;;9;;;;\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:11:7:1134#0/1\t0\tI\t35512\t255\t36M\t*\t0\t0\tATGTGATGCTTGTCAGAGGGGCAGTGAAAATCTATG\tB7A@@ABAAAAA?@AA@A@9@;8=7?;>=74;3684\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:10:1500:217#0/1\t0\tI\t36126\t255\t36M\t*\t0\t0\tAAGAAGTTGTTCGTGCCATCCACAACGGAGACATCA\tB54+48874286418588888882/42436263677\tXA:i:0\tMD:Z:35G0\tNM:i:1\n+HWI-EAS225:1:10:1400:1623#0/1\t0\tI\t36164\t255\t36M\t*\t0\t0\tATGGAAGATTGTAAGCAACTAATCACTGGTAAGCAA\tBB@??>;<?A><?9279?:>>>>>::;88::99999\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:10:1488:698#0/1\t0\tI\t36165\t255\t36M\t*\t0\t0\tTGGAAGATTGTAAGCAACTAATCACTGGTAAGCAAA\tB@>:7=8;?;9>;9<<<<<;;<<69<<<<858<:76\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:10:1740:1132#0/1\t0\tI\t36647\t255\t36M\t*\t0\t0\tAAAGGATTTAAGTGAACTACCTAGAGCTGAAAAAAC\tBA?@;(=9>>>=6(0?3999677459:766557:9#\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:10:1525:567#0/1\t0\tI\t36679\t255\t36M\t*\t0\t0\tAAACCGATCTGACTAATATTTTGATTAGAACAGCAG\tB@;-?B@AAA@?8B>@A>BAAA;;;A@?>>?>>>>>\tXA:i:1\tMD:Z:8T27\tNM:i:1\n+HWI-EAS225:1:11:47:610#0/1\t0\tI\t38579\t255\t36M\t*\t0\t0\tCGGAAATTGAAAACCCGCTATACTACGAAGAGCAGC\tBAAA@@A@>?>>>>?>?8;<=>9=>:::886/697:\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:10:1313:1915#0/1\t0\tI\t47769\t255\t36M\t*\t0\t0\tCGTGTGTCATTCATTCCTAAGATTCAAAATGATTTC\t<;BBBB9BBBBCBBBBBBABBBAAAAABB@A>B@@5\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:10:1542:1844#0/1\t0\tI\t47901\t255\t36M\t*\t0\t0\tACTTTCGGCAAGAACCTGAATGCAAACACAAATACA\tBBA<@@???@@?8<<;@?@@>??@>@?A;>?>;>>>\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:10:1690:1175#0/1\t0\tI\t47939\t255\t36M\t*\t0\t0\tAAGACTCTTTGATGATGAGACTAGTTCATCCTCTTT\tBBA?ABBAB@?@=<9?@===<?8?;@@>>?9=>??@\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:10:1629:1866#0/1\t0\tI\t48153\t255\t36M\t*\t0\t0\tGAATCAAGTCTCTACCTGGGGAATTCCAATAAATCA\tBBBBBA>@>AAA@A@A?>>@?=<?A?@@>@>>>@??\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:10:1411:1429#0/1\t0\tI\t51101\t255\t36M\t*\t0\t0\tGAAAGGGTGAATTATGGCACGGCATATTTACGATGA\t@==<A@?5>=;;>@>>>?>?29>:4>?=6<9@83;:\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:10:1586:129#0/1\t16\tI\t52586\t255\t36M\t*\t0\t0\tTTAGATGCCATTTGATTGGAATAATTTTGTTCTGCT\t<>44533:;97.7>@7><7<@>@A@AABAB@A@B@B\tXA:i:1\tMD:Z:18A17\tNM:i:1\n+HWI-EAS225:1:10:1504:470#0/1\t16\tI\t55187\t255\t36M\t*\t0\t0\tGATGTGGCTGCTCTGGGGTTTCCATGGCAATCAGTG\t77678876776788877775457777666752577B\tXA:i:1\tMD:Z:27A8\tNM:i:1\n+HWI-EAS225:1:11:58:66#0/1\t16\tI\t56517\t255\t36M\t*\t0\t0\tAAGATTTAGACCCTCGATAATCGAATTTCTAAAATC\t#####63779999:::<<<:::<<>:::>>@AA??B\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:10:1460:1769#0/1\t16\tI\t57422\t255\t36M\t*\t0\t0\tTCTCAGTGACCAAATAAAGAAACTTATATATGTTCC\t=99711<55(55<??;<>>??<@??@BBA?BBB@?B\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:10:1517:636#0/1\t0\tI\t57485\t255\t36M\t*\t0\t0\tATCGTTAGCAGCAAACGAAAGTGGAAAGAGAATAAT\tBB@=9A@@@9:>8<<?A<;=<6>@9<79:=7====@\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:10:1414:1008#0/1\t0\tI\t57518\t255\t36M\t*\t0\t0\tAATGTTCAAGAGAGCAATGAGCACAGATGGCCCCGT\tB?83/:=;=876=86;7<:86:=3(38(/5)8=856\tXA:i:0\tMD:Z:30T5\tNM:i:1\n+HWI-EAS225:1:10:1645:950#0/1\t16\tI\t58310\t255\t36M\t*\t0\t0\tCATGTTGGGAAGTGTACCTCACGGCTTGAGGTCCTT\t8898988889889988786686898=688=8<<@@B\tXA:i:0\tMD:Z:36\tNM:i:0\n+HWI-EAS225:1:10:1490:1708#0/1\t16\tI\t58435\t255\t36M\t*\t0\t0\tGTGGTCCATAGTCTAGTAGTGCGTAACATTGTCGAT'..b"1320#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tAAGAAGAGCTACGGCACTTCTGCCGTTACCGGCGGT\t>'258665############################\tXM:i:0\n+HWI-EAS225:1:11:65:1738#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tGAGACAGTATATTCGAAATGGCCGGTGTTAAAGCTT\tAA@??@@4>>>:@>?9>9>::>>?>:9/5;9:::::\tXM:i:0\n+HWI-EAS225:1:11:65:659#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tAATTGAACCGATACGAAAATAAAAAAAAAAAAAAAA\tBBBCB@B@5@>@?@A?A@>@@B@?*A>@>@6>?=??\tXM:i:0\n+HWI-EAS225:1:11:65:762#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tAGTACCAACGTCAGGACCAATGCTAGTACTAATGCG\tBB@B@7>?6?:?<;;@:;89@@??9>6:88698786\tXM:i:0\n+HWI-EAS225:1:11:66:1105#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tAAACCTCTGACATCTCGTATGCCGTCTTCTGCTTGA\tBB?:>99<87:>:>>9>8:878798998768884/4\tXM:i:0\n+HWI-EAS225:1:11:66:1289#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tAAAGCTGTTTAGTGAGGATGAAGGACAATGGGAGAC\t7@;?9@A?@?(<:A7<<117:>=68:8;<88;3687\tXM:i:0\n+HWI-EAS225:1:11:66:1547#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tAAGGGATATACTGGCAAACATATGATAAAAACAACA\tBCB@@?ABBB?A>@>@@?<?@?A@>@>;><:<:<::\tXM:i:0\n+HWI-EAS225:1:11:66:2018#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tTGAATACACCATTAACTTGCACAAAAGATTGCACGG\t>>A@>@=:=;>;;;?91(-6;>>;;;;;8:99;;96\tXM:i:0\n+HWI-EAS225:1:11:67:1064#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tATTAAATGAACATGGTATTAGTATACTCCCAGCTAA\tBCB@@;=::<<<::9;<;;:::<<<87999334:8:\tXM:i:0\n+HWI-EAS225:1:11:67:1411#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tAGGGATCCTGTACTTCAACTTCCATACTTGAAAGAT\tB@@:=??;@?=??A???@;?=<<@><78;579688;\tXM:i:0\n+HWI-EAS225:1:11:67:742#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tAGATTCGCTGTCAAACTTTCAACAACGGATCTCTTG\tBB@A?>:?>84>=68:::>9857718747::98:##\tXM:i:0\n+HWI-EAS225:1:11:68:1280#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tAAAGTCCAGAGAGACGGATACGGAAGAGACTGTGTA\tBB?7899>98:887C579:9988988899:888989\tXM:i:0\n+HWI-EAS225:1:11:68:656#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tAAAATATGGTAGACTTGACAATTAAGGGTGATGATG\tBCB@B>@=<6>>:>?<;.889??7886949:<>399\tXM:i:0\n+HWI-EAS225:1:11:68:969#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tCACATGATTTAACTCCCCCTTGACAATTGTCATATG\tB99;9;##############################\tXM:i:0\n+HWI-EAS225:1:11:69:1005#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tGGTAAATTCCATCCACAGCTAAATAGTGGCGAGAGC\tB98:96/667634#######################\tXM:i:0\n+HWI-EAS225:1:11:69:1560#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tAGCGTTATGTTTATCCTATGTGCGAGCATGGTATCA\tBBAA?A@<?@@AAA>>A?@?:???>>:>>999:::;\tXM:i:0\n+HWI-EAS225:1:11:70:1046#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tCATAGATGACGGAATGACTGGAAAAAAAAAAAAAAA\tBCCCBBBAAABBAABAB?BAA@AAAA?@@@=@>>?@\tXM:i:0\n+HWI-EAS225:1:11:70:1197#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tATAAATTAGAGGAAAATGGCTAACTTGCGTACTCAA\tBC@@@A@>><=:9998:7:<<<;;9::9949<9979\tXM:i:0\n+HWI-EAS225:1:11:70:1268#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tCAGATAAGACTCCTGATATTTAAAAAAAAAAAAAAA\tBCBBBBBB?ABBAB@BBBBBBB@B@@@??>:??;99\tXM:i:0\n+HWI-EAS225:1:11:70:441#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tAACGTGAAAAGGCTTTGAACGAAGAAGCAAGCTAGA\tBB=?876;::::;>:>>:::>::::626:635:756\tXM:i:0\n+HWI-EAS225:1:11:71:612#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tAGGAGGAGGAGGACTCTGATGTGAGGAGGACTGGTG\tBB????7<766<68=8=55<7677776555652457\tXM:i:0\n+HWI-EAS225:1:11:72:947#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAA\tBBBB@A>?A>>AA@?B@>@=9AA@8?<<<?:::4:2\tXM:i:0\n+HWI-EAS225:1:11:73:1341#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tATAGATAATCAAACAAAAAAATGGCTAAGGGAACAA\tBCBA@;?9<=?<3'4:=>:=6/5.670/434.3,59\tXM:i:0\n+HWI-EAS225:1:11:74:1238#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAA\tBBBA:8:?@@?=:<;?B?@A9@AA>:<:>;7:<8##\tXM:i:0\n+HWI-EAS225:1:11:74:1895#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAC\t;@B@@;4@>A@?;39??:@=8B??9>3;########\tXM:i:0\n+HWI-EAS225:1:11:74:646#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAC\tBCBCA=?BBBAB@A@BB@BA=BBB?@=>@?>95###\tXM:i:0\n+HWI-EAS225:1:11:75:1484#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tGCGCATTGCCAAGATCAGCTCAGCAAAGTTAAGACC\t:A=9==<;89==887988:::<888:9989788886\tXM:i:0\n+HWI-EAS225:1:11:75:1494#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tATGGACGATGGACGGAAAAAAGAATTGCATGACTTG\t>B??ABB>B>BA<?;AB@BBB@5??;2:@A>>C;9=\tXM:i:0\n+HWI-EAS225:1:11:75:58#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tAATAGGACTAAAAATGGCTGGCGTTAAAGCTTACGA\tBBA@@>2;@?@@==><389505758;;788877755\tXM:i:0\n+HWI-EAS225:1:11:76:1414#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAA\t?BACCBBCCCBBC?BCCBCBBBBA9B@@BA=@@###\tXM:i:0\n+HWI-EAS225:1:11:76:762#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tCAGTATATTCGAAATGGCCGGTGTTAAAGCTTACGA\tBBB?@A?AA?@???@?==>?=9?=?=:9=9999888\tXM:i:0\n+HWI-EAS225:1:11:7:1482#0/1\t4\t*\t0\t0\t*\t*\t0\t0\tAAGGTAGCCAAATGCCTCGTCATCTAATTAGAGGGT\tBA9<;@;??8??9>8?8<68::9698:<7888:684\tXM:i:0\n"
b
diff -r 000000000000 -r 3fdeebd7e710 test-data/htseq-test_counts.tsv
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/htseq-test_counts.tsv Tue Sep 04 13:45:36 2012 -0400
b
b'@@ -0,0 +1,7129 @@\n+15S_rRNA\t0\n+21S_rRNA\t0\n+HRA1\t0\n+LSR1\t0\n+NME1\t0\n+Q0010\t0\n+Q0017\t0\n+Q0032\t0\n+Q0045\t0\n+Q0050\t0\n+Q0055\t0\n+Q0060\t0\n+Q0065\t0\n+Q0070\t0\n+Q0075\t0\n+Q0080\t0\n+Q0085\t0\n+Q0092\t0\n+Q0105\t0\n+Q0110\t0\n+Q0115\t0\n+Q0120\t0\n+Q0130\t0\n+Q0140\t0\n+Q0142\t0\n+Q0143\t0\n+Q0144\t0\n+Q0160\t0\n+Q0182\t0\n+Q0250\t0\n+Q0255\t0\n+Q0275\t0\n+Q0297\t0\n+R0010W\t3\n+R0020C\t14\n+R0030W\t1\n+R0040C\t20\n+RDN18-1\t205\n+RDN18-2\t2\n+RDN25-1\t1408\n+RDN25-2\t59\n+RDN5-1\t0\n+RDN5-2\t1\n+RDN5-3\t0\n+RDN5-4\t0\n+RDN5-5\t3\n+RDN5-6\t0\n+RDN58-1\t85\n+RDN58-2\t1\n+RNA170\t0\n+RPM1\t0\n+RPR1\t0\n+RUF5-1\t0\n+RUF5-2\t0\n+SCR1\t0\n+SRG1\t1\n+TLC1\t0\n+YAL001C\t3\n+YAL002W\t2\n+YAL003W\t31\n+YAL004W\t0\n+YAL005C\t31\n+YAL007C\t8\n+YAL008W\t0\n+YAL009W\t0\n+YAL010C\t0\n+YAL011W\t1\n+YAL012W\t6\n+YAL013W\t1\n+YAL014C\t0\n+YAL015C\t1\n+YAL016C-A\t0\n+YAL016C-B\t0\n+YAL016W\t4\n+YAL017W\t3\n+YAL018C\t0\n+YAL019W\t5\n+YAL019W-A\t0\n+YAL020C\t0\n+YAL021C\t1\n+YAL022C\t0\n+YAL023C\t7\n+YAL024C\t0\n+YAL025C\t3\n+YAL026C\t4\n+YAL026C-A\t0\n+YAL027W\t0\n+YAL028W\t0\n+YAL029C\t7\n+YAL030W\t2\n+YAL031C\t2\n+YAL031W-A\t0\n+YAL032C\t0\n+YAL033W\t3\n+YAL034C\t0\n+YAL034C-B\t0\n+YAL034W-A\t1\n+YAL035W\t23\n+YAL036C\t7\n+YAL037C-A\t0\n+YAL037C-B\t0\n+YAL037W\t0\n+YAL038W\t105\n+YAL039C\t3\n+YAL040C\t5\n+YAL041W\t3\n+YAL042C-A\t0\n+YAL042W\t6\n+YAL043C\t3\n+YAL044C\t2\n+YAL044W-A\t2\n+YAL045C\t0\n+YAL046C\t0\n+YAL047C\t2\n+YAL047W-A\t0\n+YAL048C\t0\n+YAL049C\t1\n+YAL051W\t1\n+YAL053W\t4\n+YAL054C\t0\n+YAL055W\t0\n+YAL056C-A\t0\n+YAL056W\t0\n+YAL058W\t1\n+YAL059C-A\t0\n+YAL059W\t2\n+YAL060W\t4\n+YAL061W\t3\n+YAL062W\t0\n+YAL063C\t0\n+YAL063C-A\t0\n+YAL064C-A\t0\n+YAL064W\t0\n+YAL064W-B\t0\n+YAL065C\t0\n+YAL066W\t0\n+YAL067C\t0\n+YAL067W-A\t0\n+YAL068C\t0\n+YAL068W-A\t0\n+YAL069W\t0\n+YAR002C-A\t2\n+YAR002W\t2\n+YAR003W\t3\n+YAR007C\t4\n+YAR008W\t1\n+YAR009C\t0\n+YAR010C\t0\n+YAR014C\t5\n+YAR015W\t5\n+YAR018C\t0\n+YAR019C\t1\n+YAR019W-A\t0\n+YAR020C\t0\n+YAR023C\t0\n+YAR027W\t2\n+YAR028W\t0\n+YAR029W\t0\n+YAR030C\t0\n+YAR031W\t0\n+YAR033W\t0\n+YAR035C-A\t0\n+YAR035W\t0\n+YAR042W\t2\n+YAR047C\t0\n+YAR050W\t0\n+YAR053W\t0\n+YAR060C\t0\n+YAR061W\t0\n+YAR062W\t0\n+YAR064W\t0\n+YAR066W\t0\n+YAR068W\t11\n+YAR069C\t0\n+YAR070C\t0\n+YAR071W\t1\n+YAR073W\t0\n+YAR075W\t0\n+YBL001C\t2\n+YBL002W\t6\n+YBL003C\t20\n+YBL004W\t8\n+YBL005W\t2\n+YBL005W-A\t0\n+YBL005W-B\t3\n+YBL006C\t1\n+YBL006W-A\t0\n+YBL007C\t4\n+YBL008W\t2\n+YBL008W-A\t0\n+YBL009W\t3\n+YBL010C\t1\n+YBL011W\t4\n+YBL012C\t0\n+YBL013W\t0\n+YBL014C\t0\n+YBL015W\t0\n+YBL016W\t0\n+YBL017C\t2\n+YBL018C\t2\n+YBL019W\t0\n+YBL020W\t0\n+YBL021C\t0\n+YBL022C\t11\n+YBL023C\t0\n+YBL024W\t11\n+YBL025W\t0\n+YBL026W\t0\n+YBL027W\t92\n+YBL028C\t5\n+YBL029C-A\t0\n+YBL029W\t0\n+YBL030C\t6\n+YBL031W\t1\n+YBL032W\t7\n+YBL033C\t1\n+YBL034C\t2\n+YBL035C\t2\n+YBL036C\t1\n+YBL037W\t3\n+YBL038W\t3\n+YBL039C\t6\n+YBL039C-A\t0\n+YBL039W-B\t0\n+YBL040C\t0\n+YBL041W\t9\n+YBL042C\t4\n+YBL043W\t0\n+YBL044W\t0\n+YBL045C\t3\n+YBL046W\t2\n+YBL047C\t6\n+YBL048W\t0\n+YBL049W\t0\n+YBL050W\t1\n+YBL051C\t6\n+YBL052C\t2\n+YBL053W\t0\n+YBL054W\t0\n+YBL055C\t0\n+YBL056W\t7\n+YBL057C\t5\n+YBL058W\t5\n+YBL059C-A\t4\n+YBL059W\t0\n+YBL060W\t1\n+YBL061C\t2\n+YBL062W\t0\n+YBL063W\t2\n+YBL064C\t1\n+YBL065W\t0\n+YBL066C\t0\n+YBL067C\t4\n+YBL068W\t1\n+YBL068W-A\t0\n+YBL069W\t1\n+YBL070C\t0\n+YBL071C\t0\n+YBL071C-B\t0\n+YBL071W-A\t4\n+YBL072C\t24\n+YBL073W\t0\n+YBL074C\t0\n+YBL075C\t0\n+YBL076C\t18\n+YBL077W\t0\n+YBL078C\t1\n+YBL079W\t2\n+YBL080C\t0\n+YBL081W\t3\n+YBL082C\t0\n+YBL083C\t0\n+YBL084C\t0\n+YBL085W\t2\n+YBL086C\t4\n+YBL087C\t6\n+YBL088C\t2\n+YBL089W\t1\n+YBL090W\t2\n+YBL091C\t2\n+YBL091C-A\t0\n+YBL092W\t30\n+YBL093C\t1\n+YBL094C\t0\n+YBL095W\t2\n+YBL096C\t0\n+YBL097W\t2\n+YBL098W\t1\n+YBL099W\t3\n+YBL100C\t0\n+YBL100W-A\t0\n+YBL100W-B\t1\n+YBL100W-C\t0\n+YBL101C\t2\n+YBL102W\t1\n+YBL103C\t0\n+YBL104C\t3\n+YBL105C\t5\n+YBL106C\t1\n+YBL107C\t4\n+YBL107W-A\t0\n+YBL108C-A\t0\n+YBL108W\t0\n+YBL109W\t0\n+YBL111C\t0\n+YBL112C\t0\n+YBL113C\t2\n+YBL113W-A\t0\n+YBR001C\t3\n+YBR002C\t0\n+YBR003W\t3\n+YBR004C\t1\n+YBR005W\t0\n+YBR006W\t0\n+YBR007C\t3\n+YBR008C\t0\n+YBR009C\t23\n+YBR010W\t16\n+YBR011C\t17\n+YBR012C\t0\n+YBR012W-A\t0\n+YBR012W-B\t7\n+YBR013C\t0\n+YBR014C\t1\n+YBR015C\t5\n+YBR016W\t0\n+YBR017C\t4\n+YBR018C\t0\n+YBR019C\t0\n+YBR020W\t0\n+YBR021W\t0\n+YBR022W\t3\n+YBR023C\t2\n+YBR024W\t1\n+YBR025C\t17\n+YBR026C\t0\n+YBR027C\t0\n+YBR028C\t3\n+YBR029C\t3\n+YBR030W\t0\n+YBR031W\t65\n+YBR032W\t0\n+YBR033W\t0\n+YBR034C\t4\n+YBR035C\t2\n+YBR036C\t2\n+YBR037C\t1\n+YBR038W\t2\n+YBR039W\t10\n+YBR040W\t0\n+YBR041W\t1\n+YBR042C\t0\n+YBR043C\t1\n+YBR044C\t0\n+YBR045C\t0\n+YBR046C\t1\n+YBR047W\t0\n+YBR048W\t35\n+YBR049C\t9\n+YBR050C\t'..b'0\n+snR7-S\t0\n+snR70\t0\n+snR71\t0\n+snR72\t0\n+snR73\t0\n+snR74\t0\n+snR75\t0\n+snR76\t1\n+snR77\t1\n+snR78\t0\n+snR79\t0\n+snR8\t0\n+snR80\t0\n+snR81\t0\n+snR82\t0\n+snR83\t0\n+snR84\t0\n+snR85\t0\n+snR86\t1\n+snR87\t0\n+snR9\t1\n+tA(AGC)D\t0\n+tA(AGC)F\t0\n+tA(AGC)G\t0\n+tA(AGC)H\t0\n+tA(AGC)J\t0\n+tA(AGC)K1\t0\n+tA(AGC)K2\t0\n+tA(AGC)L\t0\n+tA(AGC)M1\t0\n+tA(AGC)M2\t0\n+tA(AGC)P\t0\n+tA(UGC)A\t0\n+tA(UGC)E\t0\n+tA(UGC)G\t0\n+tA(UGC)L\t0\n+tA(UGC)O\t0\n+tA(UGC)Q\t0\n+tC(GCA)B\t0\n+tC(GCA)G\t0\n+tC(GCA)P1\t0\n+tC(GCA)P2\t0\n+tC(GCA)Q\t0\n+tD(GUC)B\t0\n+tD(GUC)D\t0\n+tD(GUC)G1\t0\n+tD(GUC)G2\t0\n+tD(GUC)I1\t0\n+tD(GUC)I2\t0\n+tD(GUC)J1\t0\n+tD(GUC)J2\t0\n+tD(GUC)J3\t0\n+tD(GUC)J4\t0\n+tD(GUC)K\t0\n+tD(GUC)L1\t0\n+tD(GUC)L2\t0\n+tD(GUC)M\t0\n+tD(GUC)N\t0\n+tD(GUC)O\t0\n+tD(GUC)Q\t0\n+tE(CUC)D\t0\n+tE(CUC)I\t0\n+tE(UUC)B\t0\n+tE(UUC)C\t0\n+tE(UUC)E1\t0\n+tE(UUC)E2\t0\n+tE(UUC)E3\t0\n+tE(UUC)G1\t0\n+tE(UUC)G2\t0\n+tE(UUC)G3\t0\n+tE(UUC)I\t0\n+tE(UUC)J\t0\n+tE(UUC)K\t0\n+tE(UUC)L\t0\n+tE(UUC)M\t0\n+tE(UUC)P\t0\n+tE(UUC)Q\t0\n+tF(GAA)B\t0\n+tF(GAA)D\t0\n+tF(GAA)F\t0\n+tF(GAA)G\t0\n+tF(GAA)H1\t0\n+tF(GAA)H2\t0\n+tF(GAA)M\t0\n+tF(GAA)N\t0\n+tF(GAA)P1\t0\n+tF(GAA)P2\t0\n+tF(GAA)Q\t0\n+tG(CCC)D\t0\n+tG(CCC)O\t0\n+tG(GCC)B\t0\n+tG(GCC)C\t0\n+tG(GCC)D1\t0\n+tG(GCC)D2\t0\n+tG(GCC)E\t0\n+tG(GCC)F1\t0\n+tG(GCC)F2\t0\n+tG(GCC)G1\t0\n+tG(GCC)G2\t0\n+tG(GCC)J1\t0\n+tG(GCC)J2\t0\n+tG(GCC)M\t0\n+tG(GCC)O1\t0\n+tG(GCC)O2\t0\n+tG(GCC)P1\t0\n+tG(GCC)P2\t0\n+tG(UCC)G\t0\n+tG(UCC)N\t0\n+tG(UCC)O\t0\n+tG(UCC)Q\t0\n+tH(GUG)E1\t0\n+tH(GUG)E2\t0\n+tH(GUG)G1\t0\n+tH(GUG)G2\t0\n+tH(GUG)H\t0\n+tH(GUG)K\t0\n+tH(GUG)M\t0\n+tH(GUG)Q\t0\n+tI(AAU)B\t0\n+tI(AAU)D\t0\n+tI(AAU)E1\t0\n+tI(AAU)E2\t0\n+tI(AAU)G\t0\n+tI(AAU)I1\t0\n+tI(AAU)I2\t0\n+tI(AAU)L1\t0\n+tI(AAU)L2\t0\n+tI(AAU)N1\t0\n+tI(AAU)N2\t0\n+tI(AAU)P1\t0\n+tI(AAU)P2\t0\n+tI(GAU)Q\t0\n+tI(UAU)D\t0\n+tI(UAU)L\t0\n+tK(CUU)C\t0\n+tK(CUU)D1\t0\n+tK(CUU)D2\t0\n+tK(CUU)E1\t0\n+tK(CUU)E2\t0\n+tK(CUU)F\t0\n+tK(CUU)G1\t0\n+tK(CUU)G2\t0\n+tK(CUU)G3\t0\n+tK(CUU)I\t0\n+tK(CUU)J\t0\n+tK(CUU)K\t0\n+tK(CUU)M\t0\n+tK(CUU)P\t0\n+tK(UUU)D\t0\n+tK(UUU)G1\t0\n+tK(UUU)G2\t0\n+tK(UUU)K\t0\n+tK(UUU)L\t0\n+tK(UUU)O\t0\n+tK(UUU)P\t0\n+tK(UUU)Q\t0\n+tL(CAA)A\t0\n+tL(CAA)C\t0\n+tL(CAA)D\t0\n+tL(CAA)G1\t0\n+tL(CAA)G2\t0\n+tL(CAA)G3\t0\n+tL(CAA)K\t0\n+tL(CAA)L\t0\n+tL(CAA)M\t0\n+tL(CAA)N\t0\n+tL(GAG)G\t0\n+tL(UAA)B1\t0\n+tL(UAA)B2\t0\n+tL(UAA)D\t0\n+tL(UAA)J\t0\n+tL(UAA)K\t0\n+tL(UAA)L\t0\n+tL(UAA)N\t0\n+tL(UAA)Q\t0\n+tL(UAG)J\t0\n+tL(UAG)L1\t0\n+tL(UAG)L2\t0\n+tM(CAU)C\t0\n+tM(CAU)D\t0\n+tM(CAU)E\t0\n+tM(CAU)J1\t0\n+tM(CAU)J2\t0\n+tM(CAU)J3\t0\n+tM(CAU)M\t0\n+tM(CAU)O1\t0\n+tM(CAU)O2\t0\n+tM(CAU)P\t0\n+tM(CAU)Q1\t0\n+tM(CAU)Q2\t0\n+tN(GUU)C\t0\n+tN(GUU)F\t0\n+tN(GUU)G\t0\n+tN(GUU)K\t0\n+tN(GUU)L\t0\n+tN(GUU)N1\t0\n+tN(GUU)N2\t0\n+tN(GUU)O1\t0\n+tN(GUU)O2\t0\n+tN(GUU)P\t0\n+tN(GUU)Q\t0\n+tP(AGG)C\t0\n+tP(AGG)N\t0\n+tP(UGG)A\t0\n+tP(UGG)F\t0\n+tP(UGG)H\t0\n+tP(UGG)L\t0\n+tP(UGG)M\t0\n+tP(UGG)N1\t0\n+tP(UGG)N2\t0\n+tP(UGG)O1\t0\n+tP(UGG)O2\t0\n+tP(UGG)O3\t0\n+tP(UGG)Q\t0\n+tQ(CUG)M\t0\n+tQ(UUG)B\t0\n+tQ(UUG)C\t0\n+tQ(UUG)D1\t0\n+tQ(UUG)D2\t0\n+tQ(UUG)D3\t0\n+tQ(UUG)E1\t0\n+tQ(UUG)E2\t0\n+tQ(UUG)H\t0\n+tQ(UUG)L\t0\n+tQ(UUG)Q\t0\n+tR(ACG)D\t0\n+tR(ACG)E\t0\n+tR(ACG)J\t0\n+tR(ACG)K\t0\n+tR(ACG)L\t0\n+tR(ACG)O\t0\n+tR(ACG)Q2\t0\n+tR(CCG)L\t0\n+tR(CCU)J\t0\n+tR(UCU)B\t0\n+tR(UCU)D\t0\n+tR(UCU)E\t0\n+tR(UCU)G1\t0\n+tR(UCU)G2\t0\n+tR(UCU)G3\t0\n+tR(UCU)J1\t0\n+tR(UCU)J2\t0\n+tR(UCU)K\t0\n+tR(UCU)M1\t0\n+tR(UCU)M2\t0\n+tR(UCU)Q1\t0\n+tS(AGA)A\t0\n+tS(AGA)B\t0\n+tS(AGA)D1\t0\n+tS(AGA)D2\t0\n+tS(AGA)D3\t0\n+tS(AGA)E\t0\n+tS(AGA)G\t0\n+tS(AGA)H\t0\n+tS(AGA)J\t0\n+tS(AGA)L\t0\n+tS(AGA)M\t0\n+tS(CGA)C\t0\n+tS(GCU)F\t0\n+tS(GCU)L\t0\n+tS(GCU)O\t0\n+tS(GCU)Q1\t0\n+tS(UGA)E\t0\n+tS(UGA)I\t0\n+tS(UGA)P\t0\n+tS(UGA)Q2\t0\n+tT(AGU)B\t0\n+tT(AGU)C\t0\n+tT(AGU)D\t0\n+tT(AGU)H\t0\n+tT(AGU)I1\t0\n+tT(AGU)I2\t0\n+tT(AGU)J\t0\n+tT(AGU)N1\t0\n+tT(AGU)N2\t0\n+tT(AGU)O1\t0\n+tT(AGU)O2\t0\n+tT(CGU)K\t0\n+tT(UGU)G1\t0\n+tT(UGU)G2\t0\n+tT(UGU)H\t0\n+tT(UGU)P\t0\n+tT(UGU)Q1\t0\n+tT(XXX)Q2\t0\n+tV(AAC)E1\t0\n+tV(AAC)E2\t0\n+tV(AAC)G1\t0\n+tV(AAC)G2\t0\n+tV(AAC)G3\t0\n+tV(AAC)H\t0\n+tV(AAC)J\t0\n+tV(AAC)K1\t0\n+tV(AAC)K2\t0\n+tV(AAC)L\t0\n+tV(AAC)M1\t0\n+tV(AAC)M2\t0\n+tV(AAC)M3\t0\n+tV(AAC)O\t0\n+tV(CAC)D\t0\n+tV(CAC)H\t0\n+tV(UAC)B\t0\n+tV(UAC)D\t0\n+tV(UAC)Q\t0\n+tW(CCA)G1\t0\n+tW(CCA)G2\t0\n+tW(CCA)J\t0\n+tW(CCA)K\t0\n+tW(CCA)M\t0\n+tW(CCA)P\t0\n+tW(UCA)Q\t0\n+tX(XXX)D\t0\n+tY(GUA)D\t0\n+tY(GUA)F1\t0\n+tY(GUA)F2\t0\n+tY(GUA)J1\t0\n+tY(GUA)J2\t0\n+tY(GUA)M1\t0\n+tY(GUA)M2\t0\n+tY(GUA)O\t0\n+tY(GUA)Q\t0\n+no_feature\t1437\n+ambiguous\t67\n+too_low_aQual\t0\n+not_aligned\t1336\n+alignment_not_unique\t0\n'
b
diff -r 000000000000 -r 3fdeebd7e710 tool_data_table_conf.xml.sample
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/tool_data_table_conf.xml.sample Tue Sep 04 13:45:36 2012 -0400
b
@@ -0,0 +1,8 @@
+<!-- Use the file tool_data_table_conf.xml.oldlocstyle if you don't want to update your loc files as changed in revision 4550:535d276c92bc-->
+<tables>
+    <!-- Location of SAMTools indexes and other files -->
+    <table name="sam_fa_indexes" comment_char="#">
+        <columns>line_type, value, path</columns>
+        <file path="tool-data/sam_fa_indices.loc" />
+    </table>
+</tables>
b
diff -r 000000000000 -r 3fdeebd7e710 tool_dependencies.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/tool_dependencies.xml Tue Sep 04 13:45:36 2012 -0400
b
@@ -0,0 +1,44 @@
+<?xml version="1.0"?>
+<tool_dependency>
+    <package name="htseq" version="0.5.3p9">
+        <install version="1.0">
+            <actions>
+                <action type="download_by_url">http://pypi.python.org/packages/source/H/HTSeq/HTSeq-0.5.3p9.tar.gz</action>
+                <action type="shell_command">mkdir -p $INSTALL_DIR/lib/python</action>
+                <action type="set_environment">
+                    <environment_variable name="PYTHONPATH" action="append_to">$INSTALL_DIR/lib/python</environment_variable>
+                </action>
+                <action type="shell_command">export PYTHONPATH=$PYTHONPATH:$INSTALL_DIR/lib/python &amp;&amp; python setup.py install --home $INSTALL_DIR --install-scripts $INSTALL_DIR/bin</action>
+                <action type="set_environment">
+                    <environment_variable name="PATH" action="prepend_to">$INSTALL_DIR/bin</environment_variable>
+                </action>
+            </actions>
+        </install>
+        <readme>
+            Installation of HTSeq requires Python 2.5+ (does not yet work with Python 3), and the Nympy Python package.
+        </readme>
+    </package>
+    <package name="samtools" version="0.1.18">
+        <install version="1.0">
+            <actions>
+                <action type="download_by_url">http://sourceforge.net/projects/samtools/files/samtools/0.1.18/samtools-0.1.18.tar.bz2</action>
+                <action type="shell_command">sed -i.bak -e 's/-lcurses/-lncurses/g' Makefile</action>
+                <action type="shell_command">make</action>
+                <action type="move_file">
+                    <source>samtools</source>
+                    <destination>$INSTALL_DIR/bin</destination>
+                </action>
+                <action type="move_file">
+                    <source>misc/maq2sam-long</source>
+                    <destination>$INSTALL_DIR/bin</destination>
+                </action>
+                <action type="set_environment">
+                    <environment_variable name="PATH" action="prepend_to">$INSTALL_DIR/bin</environment_variable>
+                </action>
+            </actions>
+        </install>
+        <readme>
+            Compiling SAMtools requires the ncurses and zlib development libraries.
+        </readme>
+    </package>
+</tool_dependency>