Repository 'qualimap_bamqc'
hg clone https://toolshed.g2.bx.psu.edu/repos/iuc/qualimap_bamqc

Changeset 0:ac607906f10a (2019-10-10)
Next changeset 1:4a89c6f84425 (2019-10-29)
Commit message:
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/qualimap commit b4d43001cc0caa14d760c347fa1c416929f769b2"
added:
qualimap_bamqc.xml
qualimap_macros.xml
test-data/coverage_across_reference.txt
test-data/coverage_histogram.txt
test-data/duplication_rate_histogram.txt
test-data/features.gtf
test-data/genome_fraction_coverage.txt
test-data/genome_results_default.txt
test-data/genome_results_inside_features.txt
test-data/genome_results_outside_features.txt
test-data/mapped_reads_clipping_profile.txt
test-data/mapped_reads_gc-content_distribution.txt
test-data/mapped_reads_nucleotide_content.txt
test-data/mapping_quality_across_reference.txt
test-data/mapping_quality_histogram.txt
test-data/mouse_counts_ensemble_1000_6.tsv
test-data/per_base_coverage_default.txt
test-data/per_base_coverage_inside_features.txt
test-data/per_base_coverage_outside_features.txt
test-data/rnaseq_qc_counts_custom.txt
test-data/rnaseq_qc_results_custom.txt
test-data/rnaseq_qc_results_default.txt
test-data/test_mapped_reads.bam
b
diff -r 000000000000 -r ac607906f10a qualimap_bamqc.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/qualimap_bamqc.xml Thu Oct 10 17:42:04 2019 -0400
[
b'@@ -0,0 +1,474 @@\n+<tool id="qualimap_bamqc" name="QualiMap BamQC" version="@VERSION@">\n+    <macros>\n+        <import>qualimap_macros.xml</import>\n+    </macros>\n+    <expand macro="requirements" />\n+    <expand macro="version_command" />\n+    <command detect_errors="exit_code"><![CDATA[\n+        #import os\n+        @SET_JAVA_OPTS@ &&\n+\n+        ## Set some default file names and paths\n+        ## where we expect tool output to end up.\n+        ## Note that most of these need to be overwritten if the user is\n+        ## interested in regions *outside* those defined in a custom regions\n+        ## file.\n+        #set $out_dir = \'results\'\n+        #set $report_name = \'qualimapReport\'\n+        #set $summary_report = \'genome_results.txt\'\n+        #set $coverage_file = os.path.join($out_dir, \'coverage.txt\')\n+        ## This is the only file path that qualimap does not calculate\n+        ## from $out_dir.\n+        #set $per_base_coverage_target = $coverage_file\n+\n+        qualimap bamqc\n+        -bam \'$input1\' -outdir results -outformat html\n+        --collect-overlap-pairs\n+        #if str($stats_regions.region_select) == \'custom_regions\':\n+            -gff ${stats_regions.regions}\n+            #if $stats_regions.outside_stats:\n+                #set $report_name = \'qualimapReportOutsideRegions\'\n+                #set $summary_report = \'outside_results.txt\'\n+                #set $coverage_file = os.path.join(\n+                    $out_dir, \'outside_coverage.txt\'\n+                )\n+                #if $per_base_coverage:\n+                    #set $per_base_coverage_target = \'/dev/null\'\n+                #end if\n+                ${stats_regions.outside_stats}\n+            #end if\n+        #end if\n+        #if $per_base_coverage:\n+            $per_base_coverage $per_base_coverage_target\n+        #end if\n+        -nw ${plot_specific.n_bins}\n+        ${plot_specific.paint_chromosome_limits}\n+        #if $plot_specific.genome_gc_distr:\n+          --genome-gc-distr ${plot_specific.genome_gc_distr}\n+        #end if\n+        -hm ${plot_specific.homopolymer_size}\n+\n+        #if $duplicate_skipping:\n+          --skip-duplicated\n+          #if str($duplicate_skipping) == \'0,1\':\n+            --skip-dup-mode 2\n+          #else:\n+            --skip-dup-mode ${duplicate_skipping}\n+          #end if\n+        #end if\n+        -nt \\${GALAXY_SLOTS:-1} &&\n+\n+        #if $per_base_coverage:\n+            mv $coverage_file \'$output_per_base_coverage\' &&\n+        #end if\n+        @MASSAGE_OUTPUT@\n+    ]]></command>\n+    <inputs>\n+        <param argument="-bam" name="input1" type="data" format="bam"\n+        label="Mapped reads input dataset" />\n+        <conditional name="stats_regions">\n+            <param name="region_select" type="select" label="Reference genome regions to calculate mapping statistics for">\n+                <option value="all">All (whole genome)</option>\n+                <option value="custom_regions">Select regions</option>\n+            </param>\n+            <when value="all" />\n+            <when value="custom_regions">\n+                <param argument="-gff" name="regions" type="data" format="gff,gtf,bed"\n+                label="Dataset specifying regions" />\n+                <param argument="-os" name="outside_stats" type="boolean" truevalue="--outside-stats" falsevalue="" checked="false"\n+                label="Invert regions"\n+                help="If selected, report read statistics *outside* the regions in the regions file." />\n+            </when>\n+        </conditional>\n+        <param argument="-oc" name="per_base_coverage" type="boolean" truevalue="--output-genome-coverage" falsevalue="" checked="false"\n+        label="Generate per-base coverage output"\n+        help="Produce additional tabular output listing the coverage at every site (omitting only zero-coverage positions) in the selected regions of the genome. Caution: Will generate a huge dataset for anything but small input genomes or restricted regions!" />\n+        <p'..b'age.\n+\n+*Coverage Histogram (0-50X)*\n+\n+Similar to the previous plot, but in this graph genome locations with a\n+coverage greater than 50X are grouped into the last bin.\n+By doing so a higher resolution of the most common values for the coverage rate\n+is obtained.\n+\n+*Genome Fraction Coverage*\n+\n+Provides a visual way of knowing how much reference has been sequenced to at\n+least a given coverage rate.\n+This graph should be interpreted as in this example:\n+If one aims for a coverage rate of at least 25X (x-axis), how much of the\n+reference (y-axis) will be considered?\n+\n+*Duplication Rate Histogram*\n+\n+This plot shows the distribution of duplicated reads.\n+Due to several factors (*e.g.* amount of starting material, sample preparation,\n+*etc.*) it is possible that the same fragments are sequenced several times.\n+For some experiments where enrichment is used (*e.g.* ChIP-seq ) this is\n+expected to some degree.\n+For most experiments, however, a high duplication level of the reads indicates\n+some unwanted bias.\n+\n+*Mapped Reads Nucleotide Content*\n+\n+This plot shows the nucleotide content per position of the mapped reads.\n+\n+*Mapped Reads GC Content Distribution*\n+\n+This graph shows the distribution of GC-content per mapped read.\n+If compared with a precomputed genome distribution, this plot allows to check\n+if there is a shift in the GC content.\n+\n+*Mapped Reads Clipping Profile*\n+\n+Represents the percentage of clipped bases across the reads.\n+Technically, the clipping is detected via SAM format CIGAR codes \xe2\x80\x98H\xe2\x80\x99\n+(hard clipping) and \xe2\x80\x98S\xe2\x80\x99 (soft clipping).\n+In addition, the total number of clipped reads can be found in the report\n+`Summary` section.\n+\n+This plot is not shown if no clipped reads are found.\n+\n+*Homopolymer Indels*\n+\n+This bar plot shows the number of indels that are located within A, C, G and T\n+homopolymers, respectively, as well as the number of indels that are not within\n+any homopolymer. Large numbers of homopolymer indels may indicate a problem in\n+the sequencing process.\n+Technically, Qualimap identifies indels from the CIGAR code of the aligned\n+reads. Indel statistics can also be found in a dedicated section of the report\n+`Summary`.\n+\n+This graph is not shown if the sample doesn\xe2\x80\x99t contain any indels.\n+\n+*Mapping Quality Across Reference*\n+\n+This plot provides the mapping quality distribution across the reference.\n+To construct the plot, the mean mapping quality is computed for each bin.\n+\n+*Mapping Quality Histogram*\n+\n+Histogram of the number of genomic locations having a given mapping quality.\n+To construct the histogram the mean mapping quality is computed at each genome\n+position with non-zero coverage and collected.\n+According to the SAM/BAM format specifications, the range for the mapping\n+quality score is [0-255].\n+\n+*Insert Size Across Reference*\n+\n+This plot provides the insert size distribution across the reference.\n+Technically, the insert size of each pair of aligned reads is collected from\n+the SAM alignment field `TLEN`. Only positive values are taken into account.\n+To construct the plot, the mean insert size is computed for each bin.\n+\n+*Insert Size Histogram*\n+\n+Histogram of insert size distribution.\n+\n+\n+Raw Data\n+--------\n+\n+This is a *Collection* of 10 individual datasets.\n+\n+The *genome_results* dataset provides a plain-text summary of key statistics,\n+most of which can also be found in the *Summary* section of the *HTML Report*.\n+\n+The remaining 9 datasets hold the tabular raw data underlying the plots of the corresponding names in the *HTML Report*.\n+\n+\n+Per-base coverage\n+-----------------\n+\n+Optional. This is a tabular dataset listing the coverage of every base in the\n+reference genome unless that coverage is zero. Since its content is\n+uncompressed text, this dataset can easily become huge, and it is recommended\n+that you generate this dataset only for very small genomes or very limited\n+regions of larger genomes.\n+    ]]>    </help>\n+    <expand macro="citations"/>\n+</tool>\n'
b
diff -r 000000000000 -r ac607906f10a qualimap_macros.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/qualimap_macros.xml Thu Oct 10 17:42:04 2019 -0400
[
@@ -0,0 +1,40 @@
+<macros>
+    <token name="@VERSION@">2.2.2c</token>
+
+    <xml name="requirements">
+        <requirements>
+            <requirement type="package" version="@VERSION@">qualimap</requirement>
+            <yield />
+        </requirements>
+    </xml>
+
+    <xml name="version_command">
+        <version_command>qualimap --help | grep 'QualiMap v.' | cut -d '.' -f 2-</version_command>
+    </xml>
+
+    <xml name="citations">
+        <citations>
+            <citation type="doi">10.1093/bioinformatics/btv566</citation>
+            <citation type="doi">10.1093/bioinformatics/bts503</citation>
+            <yield />
+        </citations>
+    </xml>
+
+    <token name="@SET_JAVA_OPTS@">
+        export JAVA_OPTS="-Djava.awt.headless=true -Xmx\${GALAXY_MEMORY_MB:-1024}m"
+    </token>
+
+    <token name="@MASSAGE_OUTPUT@"><![CDATA[
+        sed 's|images_${report_name}/||g;s|css/||g' results/${report_name}.html > '$output_html' &&
+        mkdir '${output_html.extra_files_path}' &&
+        mv results/css/*.css '${output_html.extra_files_path}' &&
+        mv results/css/*.png '${output_html.extra_files_path}' &&
+        if [ -d results/images_${report_name} ]; then
+          mv results/images_${report_name}/* '${output_html.extra_files_path}' &&
+          for file in \$(ls -A results/raw_data_${report_name}); do mv "results/raw_data_${report_name}/\$file" `echo "results/\$file" | sed 's/(//;s/)//'`; done
+        fi
+        #if $summary_report:
+          && mv results/$summary_report results/summary_report.txt
+        #end if
+    ]]></token>
+</macros>
b
diff -r 000000000000 -r ac607906f10a test-data/coverage_across_reference.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/coverage_across_reference.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,326 @@
+#Position (bp) Coverage Std
+1.5 0.0 0.0
+3.5 0.0 0.0
+5.5 0.0 0.0
+7.5 0.0 0.0
+9.5 0.0 0.0
+11.5 0.0 0.0
+13.5 0.0 0.0
+15.5 0.0 0.0
+17.5 0.0 0.0
+19.5 0.0 0.0
+21.5 0.0 0.0
+23.5 0.0 0.0
+25.5 0.0 0.0
+27.5 0.0 0.0
+29.5 0.0 0.0
+31.5 0.0 0.0
+33.5 0.0 0.0
+35.5 0.0 0.0
+37.5 0.0 0.0
+39.5 0.0 0.0
+41.5 0.0 0.0
+43.5 0.0 0.0
+45.5 0.0 0.0
+47.5 0.0 0.0
+49.5 0.0 0.0
+51.5 0.0 0.0
+53.5 1.0 0.0
+55.5 2.5 1.5811388300841898
+57.5 3.5 1.8708286933869707
+59.5 4.0 0.0
+61.5 5.0 0.0
+63.5 5.0 0.0
+65.5 5.0 0.0
+67.5 5.0 0.0
+69.5 5.0 0.0
+71.5 5.5 2.345207879911715
+73.5 6.0 0.0
+75.5 6.0 0.0
+77.5 7.0 1.0
+79.5 8.0 0.0
+81.5 8.0 0.0
+83.5 9.0 0.0
+85.5 9.5 3.082207001484488
+87.5 10.5 3.24037034920393
+89.5 12.0 0.0
+91.5 13.5 3.6742346141747673
+93.5 15.0 1.0
+95.5 16.5 4.06201920231798
+97.5 17.5 4.183300132670378
+99.5 18.0 0.0
+101.5 21.5 4.636809247747852
+103.5 23.0 0.0
+105.5 23.5 4.847679857416329
+107.5 25.0 0.0
+109.5 25.0 0.0
+111.5 25.0 0.0
+113.5 25.0 0.0
+115.5 26.0 0.0
+117.5 26.0 0.0
+119.5 26.0 0.0
+121.5 26.5 5.1478150704935
+123.5 29.0 0.0
+125.5 32.0 0.0
+127.5 33.0 0.0
+129.5 32.5 5.70087712549569
+131.5 32.5 5.70087712549569
+133.5 32.0 0.0
+135.5 33.5 5.787918451395113
+137.5 33.0 0.0
+139.5 34.0 0.0
+141.5 36.0 1.0
+143.5 37.5 6.123724356957945
+145.5 38.0 0.0
+147.5 38.0 0.0
+149.5 38.0 0.0
+151.5 38.0 0.0
+153.5 36.0 0.0
+155.5 37.0 0.0
+157.5 36.5 6.041522986797286
+159.5 36.5 6.041522986797286
+161.5 36.5 6.041522986797286
+163.5 35.5 5.958187643906492
+165.5 36.5 6.041522986797286
+167.5 36.0 0.0
+169.5 35.0 1.0
+171.5 36.0 0.0
+173.5 35.0 0.0
+175.5 34.5 6.041522986797286
+177.5 32.0 0.0
+179.5 32.0 0.0
+181.5 31.5 5.612486080160912
+183.5 31.5 5.612486080160912
+185.5 32.0 0.0
+187.5 32.0 0.0
+189.5 31.5 5.612486080160912
+191.5 31.0 0.0
+193.5 31.5 5.612486080160912
+195.5 32.0 0.0
+197.5 31.0 1.0
+199.5 28.5 5.338539126015656
+201.5 30.5 5.522680508593631
+203.5 29.5 5.431390245600108
+205.5 31.0 0.0
+207.5 31.5 5.612486080160912
+209.5 31.0 0.0
+211.5 31.0 0.0
+213.5 31.0 0.0
+215.5 30.5 5.522680508593631
+217.5 28.0 0.0
+219.5 27.0 0.0
+221.5 28.0 0.0
+223.5 28.0 0.0
+225.5 29.0 0.0
+227.5 30.0 0.0
+229.5 30.5 5.522680508593631
+231.5 31.5 5.612486080160912
+233.5 33.5 5.787918451395113
+235.5 33.0 0.0
+237.5 32.0 0.0
+239.5 31.0 1.0
+241.5 29.5 5.431390245600108
+243.5 31.5 5.612486080160912
+245.5 29.5 5.431390245600108
+247.5 29.0 0.0
+249.5 28.5 5.338539126015656
+251.5 27.0 0.0
+253.5 27.0 0.0
+255.5 27.0 0.0
+257.5 27.0 0.0
+259.5 27.0 0.0
+261.5 27.0 0.0
+263.5 27.0 0.0
+265.5 27.0 0.0
+267.5 27.0 0.0
+269.5 27.0 0.0
+271.5 27.0 0.0
+273.5 27.0 0.0
+275.5 27.0 0.0
+277.5 27.0 0.0
+279.5 27.0 0.0
+281.5 27.0 0.0
+283.5 27.0 0.0
+285.5 27.0 0.0
+287.5 27.0 0.0
+289.5 27.0 0.0
+291.5 27.0 0.0
+293.5 27.0 0.0
+295.5 27.0 0.0
+297.5 27.0 0.0
+299.5 27.0 0.0
+301.5 0.0 0.0
+303.5 0.0 0.0
+305.5 0.0 0.0
+307.5 0.0 0.0
+309.5 0.0 0.0
+311.5 0.0 0.0
+313.5 0.0 0.0
+315.5 0.0 0.0
+317.5 0.0 0.0
+319.5 0.0 0.0
+321.5 0.0 0.0
+323.5 0.0 0.0
+325.5 0.0 0.0
+327.5 0.0 0.0
+329.5 0.0 0.0
+331.5 0.0 0.0
+333.5 0.0 0.0
+335.5 0.0 0.0
+337.5 0.0 0.0
+339.5 0.0 0.0
+341.5 0.0 0.0
+343.5 0.0 0.0
+345.5 0.0 0.0
+347.5 0.0 0.0
+349.5 0.0 0.0
+351.5 0.0 0.0
+353.5 0.0 0.0
+355.5 0.0 0.0
+357.5 0.0 0.0
+359.5 0.0 0.0
+361.5 0.0 0.0
+363.5 0.0 0.0
+365.5 0.0 0.0
+367.5 0.0 0.0
+369.5 0.0 0.0
+371.5 0.0 0.0
+373.5 0.0 0.0
+375.5 0.0 0.0
+377.5 0.0 0.0
+379.5 0.0 0.0
+381.5 0.0 0.0
+383.5 0.0 0.0
+385.5 0.0 0.0
+387.5 0.0 0.0
+389.5 0.0 0.0
+391.5 0.0 0.0
+393.5 0.0 0.0
+395.5 0.0 0.0
+397.5 0.0 0.0
+399.5 0.0 0.0
+401.5 0.0 0.0
+403.5 0.0 0.0
+405.5 0.0 0.0
+407.5 0.0 0.0
+409.5 0.0 0.0
+411.5 0.0 0.0
+413.5 0.0 0.0
+415.5 0.0 0.0
+417.5 0.0 0.0
+419.5 0.0 0.0
+421.5 0.0 0.0
+423.5 0.0 0.0
+425.5 0.0 0.0
+427.5 0.0 0.0
+429.5 0.0 0.0
+431.5 0.0 0.0
+433.5 0.0 0.0
+435.5 0.0 0.0
+437.5 0.0 0.0
+439.5 0.0 0.0
+441.5 0.0 0.0
+443.5 0.0 0.0
+445.5 0.0 0.0
+447.5 0.0 0.0
+449.5 0.0 0.0
+451.5 0.0 0.0
+453.5 0.0 0.0
+455.5 0.0 0.0
+457.5 0.0 0.0
+459.5 0.0 0.0
+461.5 0.0 0.0
+463.5 0.0 0.0
+465.5 0.0 0.0
+467.5 0.0 0.0
+469.5 0.0 0.0
+471.5 0.0 0.0
+473.5 0.0 0.0
+475.5 0.0 0.0
+477.5 0.0 0.0
+479.5 0.0 0.0
+481.5 0.0 0.0
+483.5 0.0 0.0
+485.5 0.0 0.0
+487.5 0.0 0.0
+489.5 0.0 0.0
+491.5 0.0 0.0
+493.5 0.0 0.0
+495.5 0.0 0.0
+497.5 0.0 0.0
+499.5 0.0 0.0
+501.5 0.0 0.0
+503.5 0.0 0.0
+505.5 0.0 0.0
+507.5 0.0 0.0
+509.5 0.0 0.0
+511.5 0.0 0.0
+513.5 0.0 0.0
+515.5 0.0 0.0
+517.5 0.0 0.0
+519.5 0.0 0.0
+521.5 0.0 0.0
+523.5 0.0 0.0
+525.5 0.0 0.0
+527.5 0.0 0.0
+529.5 0.0 0.0
+531.5 0.0 0.0
+533.5 0.0 0.0
+535.5 0.0 0.0
+537.5 0.0 0.0
+539.5 0.0 0.0
+541.5 0.0 0.0
+543.5 0.0 0.0
+545.5 0.0 0.0
+547.5 0.0 0.0
+549.5 0.0 0.0
+551.5 0.0 0.0
+553.5 0.0 0.0
+555.5 0.0 0.0
+557.5 0.0 0.0
+559.5 0.0 0.0
+561.5 0.0 0.0
+563.5 0.0 0.0
+565.5 0.0 0.0
+567.5 0.0 0.0
+569.5 0.0 0.0
+571.5 0.0 0.0
+573.5 0.0 0.0
+575.5 0.0 0.0
+577.5 0.0 0.0
+579.5 0.0 0.0
+581.5 0.0 0.0
+583.5 0.0 0.0
+585.5 0.0 0.0
+587.5 0.0 0.0
+589.5 0.0 0.0
+591.5 0.0 0.0
+593.5 0.0 0.0
+595.5 0.0 0.0
+597.5 0.0 0.0
+599.5 0.0 0.0
+601.5 0.0 0.0
+603.5 0.0 0.0
+605.5 0.0 0.0
+607.5 0.0 0.0
+609.5 0.0 0.0
+611.5 0.0 0.0
+613.5 0.0 0.0
+615.5 0.0 0.0
+617.5 0.0 0.0
+619.5 0.0 0.0
+621.5 0.0 0.0
+623.5 0.0 0.0
+625.5 0.0 0.0
+627.5 0.0 0.0
+629.5 0.0 0.0
+631.5 0.0 0.0
+633.5 0.0 0.0
+635.5 0.0 0.0
+637.5 0.0 0.0
+639.5 0.0 0.0
+641.5 0.0 0.0
+643.5 0.0 0.0
+645.5 0.0 0.0
+647.5 0.0 0.0
+649.5 0.0 0.0
b
diff -r 000000000000 -r ac607906f10a test-data/coverage_histogram.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/coverage_histogram.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,36 @@
+#Coverage Number of genomic locations
+0.0 102.0
+1.0 2.0
+2.0 1.0
+3.0 2.0
+4.0 3.0
+5.0 11.0
+6.0 6.0
+8.0 5.0
+9.0 3.0
+10.0 2.0
+11.0 1.0
+12.0 2.0
+13.0 1.0
+14.0 2.0
+16.0 2.0
+17.0 2.0
+18.0 3.0
+21.0 1.0
+22.0 1.0
+23.0 3.0
+24.0 1.0
+25.0 8.0
+26.0 7.0
+27.0 3.0
+28.0 8.0
+29.0 11.0
+30.0 10.0
+31.0 20.0
+32.0 27.0
+33.0 11.0
+34.0 5.0
+35.0 4.0
+36.0 13.0
+37.0 8.0
+38.0 9.0
b
diff -r 000000000000 -r ac607906f10a test-data/duplication_rate_histogram.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/duplication_rate_histogram.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,51 @@
+#Duplication rate Coverage
+1.0 57.0
+2.0 8.0
+3.0 3.0
+4.0 0.0
+5.0 0.0
+6.0 0.0
+7.0 0.0
+8.0 0.0
+9.0 0.0
+10.0 0.0
+11.0 0.0
+12.0 0.0
+13.0 0.0
+14.0 0.0
+15.0 0.0
+16.0 0.0
+17.0 0.0
+18.0 0.0
+19.0 0.0
+20.0 0.0
+21.0 0.0
+22.0 0.0
+23.0 0.0
+24.0 0.0
+25.0 0.0
+26.0 0.0
+27.0 0.0
+28.0 0.0
+29.0 0.0
+30.0 0.0
+31.0 0.0
+32.0 0.0
+33.0 0.0
+34.0 0.0
+35.0 0.0
+36.0 0.0
+37.0 0.0
+38.0 0.0
+39.0 0.0
+40.0 0.0
+41.0 0.0
+42.0 0.0
+43.0 0.0
+44.0 0.0
+45.0 0.0
+46.0 0.0
+47.0 0.0
+48.0 0.0
+49.0 0.0
+50.0 0.0
b
diff -r 000000000000 -r ac607906f10a test-data/features.gtf
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/features.gtf Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,4 @@
+test_chromosome test gene 1 300 . + . gene_id "GENE1"; gene_name "GENE1"; transcript_id "GENE1_t1";
+test_chromosome test transcript 1 300 . + . gene_id "GENE1"; gene_name "GENE1"; transcript_id "GENE1_t1";
+test_chromosome test exon 1 300 . + . gene_id "GENE1"; transcript_id "GENE1_t1"; exon_number "1"; gene_name "GENE1";
+test_chromosome test CDS 100 250 . + . gene_id "GENE1"; transcript_id "GENE1_t1"; exon_number "1"; gene_name "GENE1";
b
diff -r 000000000000 -r ac607906f10a test-data/genome_fraction_coverage.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/genome_fraction_coverage.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,52 @@
+#Coverage (X) Coverage
+1.0 66.0
+2.0 65.33333333333334
+3.0 65.0
+4.0 64.33333333333334
+5.0 63.333333333333336
+6.0 59.66666666666667
+7.0 57.66666666666667
+8.0 57.66666666666667
+9.0 56.00000000000001
+10.0 55.00000000000001
+11.0 54.33333333333334
+12.0 54.00000000000001
+13.0 53.33333333333334
+14.0 53.00000000000001
+15.0 52.33333333333334
+16.0 52.33333333333334
+17.0 51.66666666666668
+18.0 51.000000000000014
+19.0 50.000000000000014
+20.0 50.000000000000014
+21.0 50.000000000000014
+22.0 49.66666666666668
+23.0 49.33333333333334
+24.0 48.33333333333334
+25.0 48.00000000000001
+26.0 45.33333333333334
+27.0 43.00000000000001
+28.0 42.00000000000001
+29.0 39.33333333333334
+30.0 35.66666666666667
+31.0 32.33333333333334
+32.0 25.66666666666667
+33.0 16.66666666666667
+34.0 13.0
+35.0 11.333333333333329
+36.0 10.0
+37.0 5.666666666666671
+38.0 3.0
+39.0 0.0
+40.0 0.0
+41.0 0.0
+42.0 0.0
+43.0 0.0
+44.0 0.0
+45.0 0.0
+46.0 0.0
+47.0 0.0
+48.0 0.0
+49.0 0.0
+50.0 0.0
+51.0 0.0
b
diff -r 000000000000 -r ac607906f10a test-data/genome_results_default.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/genome_results_default.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,122 @@
+BamQC report
+-----------------------------------
+
+>>>>>>> Input
+
+     bam file = /tmp/tmprT4oZK/files/b/1/f/dataset_b1fbaad4-1b5a-4769-8a15-b4cec4de5731.dat
+     outfile = results/genome_results.txt
+
+
+>>>>>>> Reference
+
+     number of bases = 650 bp
+     number of contigs = 1
+
+
+>>>>>>> Globals
+
+     number of windows = 325
+
+     number of reads = 100
+     number of mapped reads = 100 (100%)
+     number of secondary alignments = 1
+
+     number of mapped bases = 12,665 bp
+     number of sequenced bases = 7,465 bp
+     number of aligned bases = 0 bp
+     number of duplicated reads (estimated) = 18
+     duplication rate = 15.85%
+
+
+>>>>>>> Insert size
+
+     mean insert size = 0
+     std insert size = 0
+     median insert size = 0
+
+
+>>>>>>> Mapping quality
+
+     mean mapping quality = 115.8402
+
+
+>>>>>>> ACTG content
+
+     number of A's = 1,910 bp (25.59%)
+     number of C's = 1,902 bp (25.48%)
+     number of T's = 1,497 bp (20.05%)
+     number of G's = 2,156 bp (28.88%)
+     number of N's = 0 bp (0%)
+
+     GC percentage = 54.36%
+
+
+>>>>>>> Mismatches and indels
+
+    general error rate = 0.0118
+    number of mismatches = 149
+
+
+>>>>>>> Coverage
+
+     mean coverageData = 19.4846X
+     std coverageData = 16.5813X
+
+     There is a 45.85% of reference with a coverageData >= 1X
+     There is a 45.38% of reference with a coverageData >= 2X
+     There is a 44.92% of reference with a coverageData >= 3X
+     There is a 44.62% of reference with a coverageData >= 4X
+     There is a 43.85% of reference with a coverageData >= 5X
+     There is a 41.85% of reference with a coverageData >= 6X
+     There is a 40.46% of reference with a coverageData >= 7X
+     There is a 39.85% of reference with a coverageData >= 8X
+     There is a 38.77% of reference with a coverageData >= 9X
+     There is a 38% of reference with a coverageData >= 10X
+     There is a 37.23% of reference with a coverageData >= 11X
+     There is a 36.92% of reference with a coverageData >= 12X
+     There is a 36.46% of reference with a coverageData >= 13X
+     There is a 36.15% of reference with a coverageData >= 14X
+     There is a 35.54% of reference with a coverageData >= 15X
+     There is a 35.23% of reference with a coverageData >= 16X
+     There is a 34.77% of reference with a coverageData >= 17X
+     There is a 34% of reference with a coverageData >= 18X
+     There is a 33.38% of reference with a coverageData >= 19X
+     There is a 33.38% of reference with a coverageData >= 20X
+     There is a 33.23% of reference with a coverageData >= 21X
+     There is a 31.85% of reference with a coverageData >= 22X
+     There is a 31.54% of reference with a coverageData >= 23X
+     There is a 30.92% of reference with a coverageData >= 24X
+     There is a 30.62% of reference with a coverageData >= 25X
+     There is a 28.62% of reference with a coverageData >= 26X
+     There is a 27.54% of reference with a coverageData >= 27X
+     There is a 26.92% of reference with a coverageData >= 28X
+     There is a 25.23% of reference with a coverageData >= 29X
+     There is a 23.38% of reference with a coverageData >= 30X
+     There is a 20.77% of reference with a coverageData >= 31X
+     There is a 16.77% of reference with a coverageData >= 32X
+     There is a 12.15% of reference with a coverageData >= 33X
+     There is a 10.15% of reference with a coverageData >= 34X
+     There is a 8.62% of reference with a coverageData >= 35X
+     There is a 7.38% of reference with a coverageData >= 36X
+     There is a 4% of reference with a coverageData >= 37X
+     There is a 1.85% of reference with a coverageData >= 38X
+     There is a 0.31% of reference with a coverageData >= 39X
+     There is a 0% of reference with a coverageData >= 40X
+     There is a 0% of reference with a coverageData >= 41X
+     There is a 0% of reference with a coverageData >= 42X
+     There is a 0% of reference with a coverageData >= 43X
+     There is a 0% of reference with a coverageData >= 44X
+     There is a 0% of reference with a coverageData >= 45X
+     There is a 0% of reference with a coverageData >= 46X
+     There is a 0% of reference with a coverageData >= 47X
+     There is a 0% of reference with a coverageData >= 48X
+     There is a 0% of reference with a coverageData >= 49X
+     There is a 0% of reference with a coverageData >= 50X
+     There is a 0% of reference with a coverageData >= 51X
+
+
+>>>>>>> Coverage per contig
+
+ test_chromosome 650 12665 19.484615384615385 16.581284452141826
+
+
b
diff -r 000000000000 -r ac607906f10a test-data/genome_results_inside_features.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/genome_results_inside_features.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,128 @@
+BamQC report
+-----------------------------------
+
+>>>>>>> Input
+
+     bam file = /tmp/tmpB4Sf0A/files/e/d/7/dataset_ed7510d4-c8f3-41fc-8cc1-3b854d9bdf14.dat
+     outfile = results/genome_results.txt
+
+
+>>>>>>> Reference
+
+     number of bases = 650 bp
+     number of contigs = 1
+
+
+>>>>>>> Globals
+
+     number of windows = 325
+
+     number of reads = 100
+     number of mapped reads = 100 (100%)
+     number of secondary alignments = 1
+
+     number of mapped bases = 6,451 bp
+     number of sequenced bases = 5,101 bp
+     number of aligned bases = 0 bp
+     number of duplicated reads (estimated) = 14
+     duplication rate = 16.18%
+
+
+>>>>>>> Globals inside
+
+     regions size = 300 (46.15%)
+     number of mapped reads = 82 (82%)
+
+
+>>>>>>> Insert size
+
+     mean insert size = 0
+     std insert size = 0
+     median insert size = 0
+
+
+>>>>>>> Mapping quality
+
+     mean mapping quality = 251.3687
+
+
+>>>>>>> ACTG content
+
+     number of A's = 1,243 bp (24.37%)
+     number of C's = 1,257 bp (24.64%)
+     number of T's = 985 bp (19.31%)
+     number of G's = 1,616 bp (31.68%)
+     number of N's = 0 bp (0%)
+
+     GC percentage = 56.32%
+
+
+>>>>>>> Mismatches and indels
+
+    general error rate = 0.0211
+    number of mismatches = 136
+
+
+>>>>>>> Coverage
+
+     mean coverageData = 21.5033X
+     std coverageData = 15.658X
+
+     There is a 66% of reference with a coverageData >= 1X
+     There is a 65.33% of reference with a coverageData >= 2X
+     There is a 65% of reference with a coverageData >= 3X
+     There is a 64.33% of reference with a coverageData >= 4X
+     There is a 63.33% of reference with a coverageData >= 5X
+     There is a 59.67% of reference with a coverageData >= 6X
+     There is a 57.67% of reference with a coverageData >= 7X
+     There is a 57.67% of reference with a coverageData >= 8X
+     There is a 56% of reference with a coverageData >= 9X
+     There is a 55% of reference with a coverageData >= 10X
+     There is a 54.33% of reference with a coverageData >= 11X
+     There is a 54% of reference with a coverageData >= 12X
+     There is a 53.33% of reference with a coverageData >= 13X
+     There is a 53% of reference with a coverageData >= 14X
+     There is a 52.33% of reference with a coverageData >= 15X
+     There is a 52.33% of reference with a coverageData >= 16X
+     There is a 51.67% of reference with a coverageData >= 17X
+     There is a 51% of reference with a coverageData >= 18X
+     There is a 50% of reference with a coverageData >= 19X
+     There is a 50% of reference with a coverageData >= 20X
+     There is a 50% of reference with a coverageData >= 21X
+     There is a 49.67% of reference with a coverageData >= 22X
+     There is a 49.33% of reference with a coverageData >= 23X
+     There is a 48.33% of reference with a coverageData >= 24X
+     There is a 48% of reference with a coverageData >= 25X
+     There is a 45.33% of reference with a coverageData >= 26X
+     There is a 43% of reference with a coverageData >= 27X
+     There is a 42% of reference with a coverageData >= 28X
+     There is a 39.33% of reference with a coverageData >= 29X
+     There is a 35.67% of reference with a coverageData >= 30X
+     There is a 32.33% of reference with a coverageData >= 31X
+     There is a 25.67% of reference with a coverageData >= 32X
+     There is a 16.67% of reference with a coverageData >= 33X
+     There is a 13% of reference with a coverageData >= 34X
+     There is a 11.33% of reference with a coverageData >= 35X
+     There is a 10% of reference with a coverageData >= 36X
+     There is a 5.67% of reference with a coverageData >= 37X
+     There is a 3% of reference with a coverageData >= 38X
+     There is a 0% of reference with a coverageData >= 39X
+     There is a 0% of reference with a coverageData >= 40X
+     There is a 0% of reference with a coverageData >= 41X
+     There is a 0% of reference with a coverageData >= 42X
+     There is a 0% of reference with a coverageData >= 43X
+     There is a 0% of reference with a coverageData >= 44X
+     There is a 0% of reference with a coverageData >= 45X
+     There is a 0% of reference with a coverageData >= 46X
+     There is a 0% of reference with a coverageData >= 47X
+     There is a 0% of reference with a coverageData >= 48X
+     There is a 0% of reference with a coverageData >= 49X
+     There is a 0% of reference with a coverageData >= 50X
+     There is a 0% of reference with a coverageData >= 51X
+
+
+>>>>>>> Coverage per contig
+
+ test_chromosome 300 6451 21.503333333333334 15.658011438947867
+
+
b
diff -r 000000000000 -r ac607906f10a test-data/genome_results_outside_features.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/genome_results_outside_features.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,128 @@
+BamQC report
+-----------------------------------
+
+>>>>>>> Input
+
+     bam file = /tmp/tmprT4oZK/files/6/4/1/dataset_6416c0ca-3ecd-4273-ae01-22f2ac60965d.dat
+     outfile = results/outside_results.txt
+
+
+>>>>>>> Reference
+
+     number of bases = 650 bp
+     number of contigs = 1
+
+
+>>>>>>> Globals
+
+     number of windows = 325
+
+     number of reads = 100
+     number of mapped reads = 100 (100%)
+     number of secondary alignments = 1
+
+     number of mapped bases = 6,214 bp
+     number of sequenced bases = 2,364 bp
+     number of aligned bases = 0 bp
+     number of duplicated reads (estimated) = 4
+     duplication rate = 14.29%
+
+
+>>>>>>> Globals outside
+
+     regions size = 350 (53.85%)
+     number of mapped reads = 18 (18%)
+
+
+>>>>>>> Insert size
+
+     mean insert size = 0
+     std insert size = 0
+     median insert size = 0
+
+
+>>>>>>> Mapping quality
+
+     mean mapping quality = 255
+
+
+>>>>>>> ACTG content
+
+     number of A's = 667 bp (28.21%)
+     number of C's = 645 bp (27.28%)
+     number of T's = 512 bp (21.66%)
+     number of G's = 540 bp (22.84%)
+     number of N's = 0 bp (0%)
+
+     GC percentage = 50.13%
+
+
+>>>>>>> Mismatches and indels
+
+    general error rate = 0.0021
+    number of mismatches = 13
+
+
+>>>>>>> Coverage
+
+     mean coverageData = 17.7543X
+     std coverageData = 16.4764X
+
+     There is a 10.58% of reference with a coverageData >= 1X
+     There is a 10.58% of reference with a coverageData >= 2X
+     There is a 10.58% of reference with a coverageData >= 3X
+     There is a 10.58% of reference with a coverageData >= 4X
+     There is a 10.58% of reference with a coverageData >= 5X
+     There is a 10.58% of reference with a coverageData >= 6X
+     There is a 10.58% of reference with a coverageData >= 7X
+     There is a 10.58% of reference with a coverageData >= 8X
+     There is a 10.58% of reference with a coverageData >= 9X
+     There is a 10.58% of reference with a coverageData >= 10X
+     There is a 10.58% of reference with a coverageData >= 11X
+     There is a 10.58% of reference with a coverageData >= 12X
+     There is a 10.58% of reference with a coverageData >= 13X
+     There is a 10.58% of reference with a coverageData >= 14X
+     There is a 10.58% of reference with a coverageData >= 15X
+     There is a 10.58% of reference with a coverageData >= 16X
+     There is a 10.58% of reference with a coverageData >= 17X
+     There is a 10.58% of reference with a coverageData >= 18X
+     There is a 10.58% of reference with a coverageData >= 19X
+     There is a 10.58% of reference with a coverageData >= 20X
+     There is a 10.58% of reference with a coverageData >= 21X
+     There is a 10.58% of reference with a coverageData >= 22X
+     There is a 10.58% of reference with a coverageData >= 23X
+     There is a 10.58% of reference with a coverageData >= 24X
+     There is a 10.58% of reference with a coverageData >= 25X
+     There is a 10.58% of reference with a coverageData >= 26X
+     There is a 10.58% of reference with a coverageData >= 27X
+     There is a 10.58% of reference with a coverageData >= 28X
+     There is a 10.58% of reference with a coverageData >= 29X
+     There is a 10.58% of reference with a coverageData >= 30X
+     There is a 10.58% of reference with a coverageData >= 31X
+     There is a 10.1% of reference with a coverageData >= 32X
+     There is a 9.62% of reference with a coverageData >= 33X
+     There is a 9.13% of reference with a coverageData >= 34X
+     There is a 7.21% of reference with a coverageData >= 35X
+     There is a 5.29% of reference with a coverageData >= 36X
+     There is a 3.37% of reference with a coverageData >= 37X
+     There is a 1.44% of reference with a coverageData >= 38X
+     There is a 0.96% of reference with a coverageData >= 39X
+     There is a 0% of reference with a coverageData >= 40X
+     There is a 0% of reference with a coverageData >= 41X
+     There is a 0% of reference with a coverageData >= 42X
+     There is a 0% of reference with a coverageData >= 43X
+     There is a 0% of reference with a coverageData >= 44X
+     There is a 0% of reference with a coverageData >= 45X
+     There is a 0% of reference with a coverageData >= 46X
+     There is a 0% of reference with a coverageData >= 47X
+     There is a 0% of reference with a coverageData >= 48X
+     There is a 0% of reference with a coverageData >= 49X
+     There is a 0% of reference with a coverageData >= 50X
+     There is a 0% of reference with a coverageData >= 51X
+
+
+>>>>>>> Coverage per contig
+
+ test_chromosome 4424 6214 1.404611211573237 3.9981701358234454
+
+
b
diff -r 000000000000 -r ac607906f10a test-data/mapped_reads_clipping_profile.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/mapped_reads_clipping_profile.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,76 @@
+#Read position (bp) Clipping profile
+0.0 16.666666666666664
+1.0 5.555555555555555
+2.0 0.0
+3.0 0.0
+4.0 0.0
+5.0 0.0
+6.0 0.0
+7.0 0.0
+8.0 0.0
+9.0 0.0
+10.0 0.0
+11.0 0.0
+12.0 0.0
+13.0 0.0
+14.0 0.0
+15.0 0.0
+16.0 0.0
+17.0 0.0
+18.0 0.0
+19.0 0.0
+20.0 0.0
+21.0 0.0
+22.0 0.0
+23.0 0.0
+24.0 0.0
+25.0 0.0
+26.0 0.0
+27.0 0.0
+28.0 0.0
+29.0 0.0
+30.0 0.0
+31.0 0.0
+32.0 0.0
+33.0 0.0
+34.0 0.0
+35.0 0.0
+36.0 0.0
+37.0 0.0
+38.0 0.0
+39.0 0.0
+40.0 0.0
+41.0 0.0
+42.0 0.0
+43.0 0.0
+44.0 0.0
+45.0 0.0
+46.0 0.0
+47.0 0.0
+48.0 0.0
+49.0 0.0
+50.0 0.0
+51.0 0.0
+52.0 0.0
+53.0 0.0
+54.0 0.0
+55.0 0.0
+56.0 0.0
+57.0 0.0
+58.0 0.0
+59.0 0.0
+60.0 0.0
+61.0 0.0
+62.0 0.0
+63.0 0.0
+64.0 0.0
+65.0 0.0
+66.0 0.0
+67.0 0.0
+68.0 0.0
+69.0 0.0
+70.0 0.0
+71.0 5.555555555555555
+72.0 11.11111111111111
+73.0 16.666666666666664
+74.0 44.44444444444444
b
diff -r 000000000000 -r ac607906f10a test-data/mapped_reads_gc-content_distribution.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/mapped_reads_gc-content_distribution.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,101 @@
+#GC Content (%) Sample HUMAN (hg19)
+1.0 0.0 0.0
+2.0 0.0 3.0E-6
+3.0 0.0 3.0E-6
+4.0 0.0 1.0E-5
+5.0 0.0 8.0E-6
+6.0 0.0 1.2000000000000002E-5
+7.0 0.0 6.999999999999999E-6
+8.0 0.0 9.999999999999999E-6
+9.0 0.0 1.6E-5
+10.0 0.0 1.5000000000000002E-5
+11.0 0.0 1.9999999999999998E-5
+12.0 0.0 1.6E-5
+13.0 0.0 1.6E-5
+14.0 0.0 2.4999999999999998E-5
+15.0 0.0 2.2000000000000003E-5
+16.0 0.0 2.6999999999999996E-5
+17.0 0.0 3.3E-5
+18.0 0.0 3.7000000000000005E-5
+19.0 0.0 4.7000000000000004E-5
+20.0 0.0 5.4999999999999995E-5
+21.0 0.0 6.400000000000001E-5
+22.0 0.0 9.3E-5
+23.0 0.0 1.66E-4
+24.0 0.0 3.03E-4
+25.0 0.0 7.059999999999999E-4
+26.0 0.0 0.0016020000000000001
+27.0 0.0 0.0033690000000000005
+28.0 0.0 0.006167999999999999
+29.0 0.0 0.010352
+30.0 0.0 0.015647
+31.0 0.0 0.022148
+32.0 0.0 0.028808
+33.0 0.0 0.035106000000000005
+34.0 0.0 0.041161
+35.0 0.0 0.047408
+36.0 0.0 0.052698999999999996
+37.0 0.0 0.057517000000000006
+38.0 0.0 0.061487999999999994
+39.0 0.0 0.063321
+40.0 0.0 0.060248
+41.0 0.0 0.055508
+42.0 0.0 0.051446000000000006
+43.0 0.0 0.048194
+44.0 0.017857142857142856 0.044423
+45.0 0.0 0.040655000000000004
+46.0 0.03571428571428571 0.03652
+47.0 0.0 0.032315
+48.0 0.017857142857142856 0.028137999999999996
+49.0 0.0 0.02458
+50.0 0.017857142857142856 0.020810000000000002
+51.0 0.0 0.017532
+52.0 0.05357142857142857 0.014506999999999997
+53.0 0.017857142857142856 0.012007
+54.0 0.14285714285714285 0.010018
+55.0 0.08928571428571429 0.008414999999999999
+56.0 0.17857142857142855 0.007258
+57.0 0.03571428571428571 0.006174999999999999
+58.0 0.08928571428571429 0.005363000000000001
+59.0 0.03571428571428571 0.0046159999999999994
+60.0 0.14285714285714285 0.003946
+61.0 0.0 0.0034579999999999993
+62.0 0.03571428571428571 0.002915
+63.0 0.017857142857142856 0.0024610000000000005
+64.0 0.05357142857142857 0.0019950000000000002
+65.0 0.0 0.0016779999999999998
+66.0 0.017857142857142856 0.0013500000000000003
+67.0 0.0 0.001069
+68.0 0.0 8.119999999999999E-4
+69.0 0.0 6.4E-4
+70.0 0.0 5.05E-4
+71.0 0.0 4.0299999999999993E-4
+72.0 0.0 3.2E-4
+73.0 0.0 2.79E-4
+74.0 0.0 2.0700000000000002E-4
+75.0 0.0 2.0099999999999998E-4
+76.0 0.0 1.54E-4
+77.0 0.0 1.12E-4
+78.0 0.0 7.699999999999999E-5
+79.0 0.0 6.3E-5
+80.0 0.0 3.7E-5
+81.0 0.0 2.8000000000000003E-5
+82.0 0.0 1.6E-5
+83.0 0.0 9.0E-6
+84.0 0.0 4.0E-6
+85.0 0.0 0.0
+86.0 0.0 1.0E-6
+87.0 0.0 0.0
+88.0 0.0 0.0
+89.0 0.0 0.0
+90.0 0.0 0.0
+91.0 0.0 0.0
+92.0 0.0 0.0
+93.0 0.0 0.0
+94.0 0.0 0.0
+95.0 0.0 0.0
+96.0 0.0 0.0
+97.0 0.0 0.0
+98.0 0.0 0.0
+99.0 0.0 0.0
+100.0 0.0 0.0
b
diff -r 000000000000 -r ac607906f10a test-data/mapped_reads_nucleotide_content.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/mapped_reads_nucleotide_content.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,76 @@
+# Position (bp) A C G T N
+0.0 34.177215189873415 21.518987341772153 26.582278481012654 17.72151898734177 0.0
+1.0 23.456790123456788 33.33333333333333 20.98765432098765 22.22222222222222 0.0
+2.0 25.609756097560975 20.73170731707317 35.36585365853659 18.29268292682927 0.0
+3.0 23.170731707317074 28.04878048780488 28.04878048780488 20.73170731707317 0.0
+4.0 31.70731707317073 26.82926829268293 24.390243902439025 17.073170731707318 0.0
+5.0 25.609756097560975 19.51219512195122 34.146341463414636 20.73170731707317 0.0
+6.0 25.609756097560975 23.170731707317074 34.146341463414636 17.073170731707318 0.0
+7.0 28.04878048780488 21.951219512195124 30.48780487804878 19.51219512195122 0.0
+8.0 23.170731707317074 18.29268292682927 39.02439024390244 19.51219512195122 0.0
+9.0 28.04878048780488 29.268292682926827 24.390243902439025 18.29268292682927 0.0
+10.0 24.390243902439025 23.170731707317074 34.146341463414636 18.29268292682927 0.0
+11.0 19.51219512195122 31.70731707317073 28.04878048780488 20.73170731707317 0.0
+12.0 34.146341463414636 23.170731707317074 28.04878048780488 14.634146341463413 0.0
+13.0 26.82926829268293 18.29268292682927 36.58536585365854 18.29268292682927 0.0
+14.0 24.390243902439025 26.82926829268293 29.268292682926827 19.51219512195122 0.0
+15.0 29.268292682926827 21.951219512195124 28.04878048780488 20.73170731707317 0.0
+16.0 31.70731707317073 25.609756097560975 29.268292682926827 13.414634146341465 0.0
+17.0 15.853658536585366 25.609756097560975 35.36585365853659 23.170731707317074 0.0
+18.0 31.70731707317073 17.073170731707318 31.70731707317073 19.51219512195122 0.0
+19.0 24.390243902439025 29.268292682926827 29.268292682926827 17.073170731707318 0.0
+20.0 21.951219512195124 25.609756097560975 29.268292682926827 23.170731707317074 0.0
+21.0 23.170731707317074 28.04878048780488 36.58536585365854 12.195121951219512 0.0
+22.0 28.04878048780488 30.48780487804878 18.29268292682927 23.170731707317074 0.0
+23.0 23.170731707317074 18.29268292682927 29.268292682926827 29.268292682926827 0.0
+24.0 25.609756097560975 19.51219512195122 35.36585365853659 19.51219512195122 0.0
+25.0 23.170731707317074 26.82926829268293 25.609756097560975 24.390243902439025 0.0
+26.0 21.951219512195124 19.51219512195122 41.46341463414634 17.073170731707318 0.0
+27.0 18.29268292682927 29.268292682926827 35.36585365853659 17.073170731707318 0.0
+28.0 29.268292682926827 30.48780487804878 25.609756097560975 14.634146341463413 0.0
+29.0 17.073170731707318 25.609756097560975 36.58536585365854 20.73170731707317 0.0
+30.0 25.609756097560975 21.951219512195124 32.926829268292686 19.51219512195122 0.0
+31.0 29.268292682926827 23.170731707317074 32.926829268292686 14.634146341463413 0.0
+32.0 29.268292682926827 21.951219512195124 24.390243902439025 24.390243902439025 0.0
+33.0 20.73170731707317 31.70731707317073 26.82926829268293 20.73170731707317 0.0
+34.0 18.29268292682927 20.73170731707317 40.243902439024396 20.73170731707317 0.0
+35.0 37.80487804878049 13.414634146341465 28.04878048780488 20.73170731707317 0.0
+36.0 28.04878048780488 28.04878048780488 28.04878048780488 15.853658536585366 0.0
+37.0 21.951219512195124 20.73170731707317 39.02439024390244 18.29268292682927 0.0
+38.0 34.146341463414636 26.82926829268293 14.634146341463413 24.390243902439025 0.0
+39.0 31.70731707317073 17.073170731707318 26.82926829268293 24.390243902439025 0.0
+40.0 26.82926829268293 24.390243902439025 31.70731707317073 17.073170731707318 0.0
+41.0 30.48780487804878 19.51219512195122 25.609756097560975 24.390243902439025 0.0
+42.0 15.853658536585366 25.609756097560975 34.146341463414636 24.390243902439025 0.0
+43.0 29.268292682926827 24.390243902439025 28.04878048780488 18.29268292682927 0.0
+44.0 26.82926829268293 24.390243902439025 31.70731707317073 17.073170731707318 0.0
+45.0 29.268292682926827 23.170731707317074 29.268292682926827 18.29268292682927 0.0
+46.0 24.390243902439025 23.170731707317074 30.48780487804878 21.951219512195124 0.0
+47.0 25.609756097560975 28.04878048780488 30.48780487804878 15.853658536585366 0.0
+48.0 29.268292682926827 25.609756097560975 29.268292682926827 15.853658536585366 0.0
+49.0 23.170731707317074 23.170731707317074 32.926829268292686 20.73170731707317 0.0
+50.0 24.390243902439025 28.04878048780488 31.70731707317073 15.853658536585366 0.0
+51.0 24.390243902439025 24.390243902439025 31.70731707317073 19.51219512195122 0.0
+52.0 18.29268292682927 26.82926829268293 34.146341463414636 20.73170731707317 0.0
+53.0 26.82926829268293 23.170731707317074 36.58536585365854 13.414634146341465 0.0
+54.0 26.82926829268293 23.170731707317074 31.70731707317073 18.29268292682927 0.0
+55.0 25.609756097560975 26.82926829268293 28.04878048780488 19.51219512195122 0.0
+56.0 10.975609756097562 32.926829268292686 40.243902439024396 15.853658536585366 0.0
+57.0 26.82926829268293 19.51219512195122 29.268292682926827 24.390243902439025 0.0
+58.0 26.82926829268293 30.48780487804878 25.609756097560975 17.073170731707318 0.0
+59.0 21.951219512195124 19.51219512195122 32.926829268292686 25.609756097560975 0.0
+60.0 20.73170731707317 28.04878048780488 32.926829268292686 18.29268292682927 0.0
+61.0 30.48780487804878 19.51219512195122 32.926829268292686 17.073170731707318 0.0
+62.0 24.390243902439025 23.170731707317074 30.48780487804878 21.951219512195124 0.0
+63.0 19.51219512195122 23.170731707317074 29.268292682926827 28.04878048780488 0.0
+64.0 25.609756097560975 31.70731707317073 24.390243902439025 18.29268292682927 0.0
+65.0 18.29268292682927 24.390243902439025 30.48780487804878 26.82926829268293 0.0
+66.0 20.73170731707317 19.51219512195122 36.58536585365854 23.170731707317074 0.0
+67.0 23.170731707317074 31.70731707317073 25.609756097560975 19.51219512195122 0.0
+68.0 15.853658536585366 24.390243902439025 42.68292682926829 17.073170731707318 0.0
+69.0 25.609756097560975 26.82926829268293 29.268292682926827 18.29268292682927 0.0
+70.0 24.390243902439025 34.146341463414636 31.70731707317073 9.75609756097561 0.0
+71.0 19.753086419753085 27.160493827160494 29.629629629629626 23.456790123456788 0.0
+72.0 25.0 37.5 17.5 20.0 0.0
+73.0 24.050632911392405 20.253164556962027 34.177215189873415 21.518987341772153 0.0
+74.0 17.56756756756757 31.08108108108108 28.37837837837838 22.972972972972975 0.0
b
diff -r 000000000000 -r ac607906f10a test-data/mapping_quality_across_reference.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/mapping_quality_across_reference.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,326 @@
+#Position (bp) mapping quality
+1.5 0.0
+3.5 0.0
+5.5 0.0
+7.5 0.0
+9.5 0.0
+11.5 0.0
+13.5 0.0
+15.5 0.0
+17.5 0.0
+19.5 0.0
+21.5 0.0
+23.5 0.0
+25.5 0.0
+27.5 0.0
+29.5 0.0
+31.5 0.0
+33.5 0.0
+35.5 0.0
+37.5 0.0
+39.5 0.0
+41.5 0.0
+43.5 0.0
+45.5 0.0
+47.5 0.0
+49.5 0.0
+51.5 0.0
+53.5 255.0
+55.5 255.0
+57.5 255.0
+59.5 255.0
+61.5 255.0
+63.5 255.0
+65.5 255.0
+67.5 255.0
+69.5 255.0
+71.5 255.0
+73.5 255.0
+75.5 255.0
+77.5 255.0
+79.5 255.0
+81.5 255.0
+83.5 255.0
+85.5 255.0
+87.5 255.0
+89.5 255.0
+91.5 236.33333333333334
+93.5 238.2
+95.5 239.72727272727272
+97.5 240.6
+99.5 241.0
+101.5 243.27906976744185
+103.5 244.04347826086956
+105.5 244.27659574468086
+107.5 244.92
+109.5 244.92
+111.5 244.92
+113.5 244.92
+115.5 245.30769230769232
+117.5 245.30769230769232
+119.5 245.30769230769232
+121.5 245.49056603773585
+123.5 246.31034482758622
+125.5 247.125
+127.5 247.36363636363637
+129.5 247.24615384615385
+131.5 247.24615384615385
+133.5 247.125
+135.5 247.47761194029852
+137.5 247.36363636363637
+139.5 247.58823529411765
+141.5 248.0
+143.5 248.28
+145.5 248.3684210526316
+147.5 248.3684210526316
+149.5 248.3684210526316
+151.5 248.3684210526316
+153.5 248.0
+155.5 248.1891891891892
+157.5 248.0958904109589
+159.5 248.0958904109589
+161.5 248.0958904109589
+163.5 247.90140845070422
+165.5 251.54794520547946
+167.5 255.0
+169.5 255.0
+171.5 255.0
+173.5 255.0
+175.5 255.0
+177.5 255.0
+179.5 255.0
+181.5 255.0
+183.5 255.0
+185.5 255.0
+187.5 255.0
+189.5 255.0
+191.5 255.0
+193.5 255.0
+195.5 255.0
+197.5 255.0
+199.5 255.0
+201.5 255.0
+203.5 255.0
+205.5 255.0
+207.5 255.0
+209.5 255.0
+211.5 255.0
+213.5 255.0
+215.5 255.0
+217.5 255.0
+219.5 255.0
+221.5 255.0
+223.5 255.0
+225.5 255.0
+227.5 255.0
+229.5 255.0
+231.5 255.0
+233.5 255.0
+235.5 255.0
+237.5 255.0
+239.5 255.0
+241.5 255.0
+243.5 255.0
+245.5 255.0
+247.5 255.0
+249.5 255.0
+251.5 0.0
+253.5 0.0
+255.5 0.0
+257.5 0.0
+259.5 0.0
+261.5 0.0
+263.5 0.0
+265.5 0.0
+267.5 0.0
+269.5 0.0
+271.5 0.0
+273.5 0.0
+275.5 0.0
+277.5 0.0
+279.5 0.0
+281.5 0.0
+283.5 0.0
+285.5 0.0
+287.5 0.0
+289.5 0.0
+291.5 0.0
+293.5 0.0
+295.5 0.0
+297.5 0.0
+299.5 0.0
+301.5 0.0
+303.5 0.0
+305.5 0.0
+307.5 0.0
+309.5 0.0
+311.5 0.0
+313.5 0.0
+315.5 0.0
+317.5 0.0
+319.5 0.0
+321.5 0.0
+323.5 0.0
+325.5 0.0
+327.5 0.0
+329.5 0.0
+331.5 0.0
+333.5 0.0
+335.5 0.0
+337.5 0.0
+339.5 0.0
+341.5 0.0
+343.5 0.0
+345.5 0.0
+347.5 0.0
+349.5 0.0
+351.5 0.0
+353.5 0.0
+355.5 0.0
+357.5 0.0
+359.5 0.0
+361.5 0.0
+363.5 0.0
+365.5 0.0
+367.5 0.0
+369.5 0.0
+371.5 0.0
+373.5 0.0
+375.5 0.0
+377.5 0.0
+379.5 0.0
+381.5 0.0
+383.5 0.0
+385.5 0.0
+387.5 0.0
+389.5 0.0
+391.5 0.0
+393.5 0.0
+395.5 0.0
+397.5 0.0
+399.5 0.0
+401.5 0.0
+403.5 0.0
+405.5 0.0
+407.5 0.0
+409.5 0.0
+411.5 0.0
+413.5 0.0
+415.5 0.0
+417.5 0.0
+419.5 0.0
+421.5 0.0
+423.5 0.0
+425.5 0.0
+427.5 0.0
+429.5 0.0
+431.5 0.0
+433.5 0.0
+435.5 0.0
+437.5 0.0
+439.5 0.0
+441.5 0.0
+443.5 0.0
+445.5 0.0
+447.5 0.0
+449.5 0.0
+451.5 0.0
+453.5 0.0
+455.5 0.0
+457.5 0.0
+459.5 0.0
+461.5 0.0
+463.5 0.0
+465.5 0.0
+467.5 0.0
+469.5 0.0
+471.5 0.0
+473.5 0.0
+475.5 0.0
+477.5 0.0
+479.5 0.0
+481.5 0.0
+483.5 0.0
+485.5 0.0
+487.5 0.0
+489.5 0.0
+491.5 0.0
+493.5 0.0
+495.5 0.0
+497.5 0.0
+499.5 0.0
+501.5 0.0
+503.5 0.0
+505.5 0.0
+507.5 0.0
+509.5 0.0
+511.5 0.0
+513.5 0.0
+515.5 0.0
+517.5 0.0
+519.5 0.0
+521.5 0.0
+523.5 0.0
+525.5 0.0
+527.5 0.0
+529.5 0.0
+531.5 0.0
+533.5 0.0
+535.5 0.0
+537.5 0.0
+539.5 0.0
+541.5 0.0
+543.5 0.0
+545.5 0.0
+547.5 0.0
+549.5 0.0
+551.5 0.0
+553.5 0.0
+555.5 0.0
+557.5 0.0
+559.5 0.0
+561.5 0.0
+563.5 0.0
+565.5 0.0
+567.5 0.0
+569.5 0.0
+571.5 0.0
+573.5 0.0
+575.5 0.0
+577.5 0.0
+579.5 0.0
+581.5 0.0
+583.5 0.0
+585.5 0.0
+587.5 0.0
+589.5 0.0
+591.5 0.0
+593.5 0.0
+595.5 0.0
+597.5 0.0
+599.5 0.0
+601.5 0.0
+603.5 0.0
+605.5 0.0
+607.5 0.0
+609.5 0.0
+611.5 0.0
+613.5 0.0
+615.5 0.0
+617.5 0.0
+619.5 0.0
+621.5 0.0
+623.5 0.0
+625.5 0.0
+627.5 0.0
+629.5 0.0
+631.5 0.0
+633.5 0.0
+635.5 0.0
+637.5 0.0
+639.5 0.0
+641.5 0.0
+643.5 0.0
+645.5 0.0
+647.5 0.0
+649.5 0.0
b
diff -r 000000000000 -r ac607906f10a test-data/mapping_quality_histogram.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/mapping_quality_histogram.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,13 @@
+#Mapping quality mapping quality
+235.0 1.0
+237.0 2.0
+239.0 2.0
+240.0 2.0
+241.0 3.0
+243.0 2.0
+244.0 12.0
+245.0 8.0
+246.0 2.0
+247.0 18.0
+248.0 23.0
+255.0 123.0
b
diff -r 000000000000 -r ac607906f10a test-data/mouse_counts_ensemble_1000_6.tsv
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/mouse_counts_ensemble_1000_6.tsv Thu Oct 10 17:42:04 2019 -0400
b
b'@@ -0,0 +1,1000 @@\n+#gene\t-GlcN01\t-GlcN02\t-GlcN03\t+GlcN01\t+GlcN02\t+GlcN03\n+ENSMUSG00000085214\t641\t542\t776\t491\t791\t368\n+ENSMUSG00000096956\t73\t17\t47\t32\t29\t63\n+ENSMUSG00000021252\t2359\t1563\t2497\t2126\t2149\t1296\n+ENSMUSG00000046840\t39\t165\t124\t83\t144\t31\n+ENSMUSG00000044062\t1\t0\t3\t6\t1\t0\n+ENSMUSG00000007777\t292\t184\t220\t228\t237\t270\n+ENSMUSG00000037361\t421\t300\t380\t316\t393\t310\n+ENSMUSG00000043644\t16\t5\t7\t6\t6\t8\n+ENSMUSG00000024442\t812\t566\t727\t569\t784\t692\n+ENSMUSG00000042208\t445\t273\t375\t305\t410\t399\n+ENSMUSG00000020831\t5\t4\t5\t6\t6\t3\n+ENSMUSG00000025731\t1259\t942\t1113\t888\t1585\t1076\n+ENSMUSG00000096232\t394\t309\t396\t261\t417\t306\n+ENSMUSG00000055312\t53\t53\t59\t43\t63\t35\n+ENSMUSG00000058706\t948\t581\t835\t716\t927\t622\n+ENSMUSG00000001418\t1440\t813\t1052\t1092\t1044\t1168\n+ENSMUSG00000028608\t305\t208\t262\t229\t273\t225\n+ENSMUSG00000089900\t3\t3\t2\t7\t4\t5\n+ENSMUSG00000058812\t51\t30\t45\t41\t32\t35\n+ENSMUSG00000089889\t15\t5\t21\t12\t22\t9\n+ENSMUSG00000087341\t82\t84\t69\t145\t96\t93\n+ENSMUSG00000060512\t621\t495\t565\t399\t691\t591\n+ENSMUSG00000087361\t26\t13\t20\t18\t26\t10\n+ENSMUSG00000051748\t5290\t4503\t5667\t8817\t5797\t1277\n+ENSMUSG00000062691\t135\t92\t121\t113\t141\t105\n+ENSMUSG00000019689\t341\t270\t361\t194\t400\t200\n+ENSMUSG00000071456\t38\t19\t25\t22\t29\t43\n+ENSMUSG00000022972\t221\t148\t167\t172\t215\t203\n+ENSMUSG00000030663\t1070\t597\t798\t671\t862\t776\n+ENSMUSG00000037960\t194\t88\t114\t106\t114\t119\n+ENSMUSG00000027637\t779\t632\t672\t570\t820\t868\n+ENSMUSG00000029600\t6\t3\t1\t6\t1\t5\n+ENSMUSG00000021023\t249\t234\t221\t184\t237\t194\n+ENSMUSG00000039195\t163\t90\t166\t216\t145\t199\n+ENSMUSG00000045237\t267\t222\t288\t198\t217\t183\n+ENSMUSG00000028441\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000097616\t41\t27\t32\t33\t33\t48\n+ENSMUSG00000097047\t12\t15\t23\t9\t12\t13\n+ENSMUSG00000085972\t24\t7\t25\t19\t25\t104\n+ENSMUSG00000037971\t480\t419\t491\t343\t577\t361\n+ENSMUSG00000046999\t0\t0\t1\t1\t0\t0\n+ENSMUSG00000044991\t176\t58\t126\t103\t123\t256\n+ENSMUSG00000086191\t5\t3\t4\t2\t6\t2\n+ENSMUSG00000097036\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000040720\t807\t491\t596\t568\t621\t542\n+ENSMUSG00000092203\t75\t46\t59\t28\t40\t65\n+ENSMUSG00000063236\t196\t91\t125\t110\t129\t127\n+ENSMUSG00000085457\t10\t5\t8\t9\t4\t4\n+ENSMUSG00000040591\t52\t34\t59\t39\t39\t24\n+ENSMUSG00000097867\t125\t152\t200\t139\t205\t82\n+ENSMUSG00000037669\t63\t43\t39\t35\t45\t45\n+ENSMUSG00000026154\t298\t233\t285\t231\t296\t211\n+ENSMUSG00000035171\t488\t227\t316\t257\t378\t300\n+ENSMUSG00000032551\t328\t126\t195\t153\t278\t192\n+ENSMUSG00000037086\t0\t0\t0\t0\t0\t1\n+ENSMUSG00000078570\t129\t50\t72\t44\t86\t109\n+ENSMUSG00000045414\t1179\t854\t813\t640\t1108\t894\n+ENSMUSG00000096537\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000043687\t72\t59\t69\t42\t80\t64\n+ENSMUSG00000063320\t47\t29\t43\t25\t43\t29\n+ENSMUSG00000020783\t491\t222\t312\t300\t407\t381\n+ENSMUSG00000097290\t253\t122\t210\t126\t150\t160\n+ENSMUSG00000028327\t2650\t1496\t2279\t2932\t2340\t1618\n+ENSMUSG00000087382\t433\t112\t307\t114\t140\t121\n+ENSMUSG00000033688\t11175\t8162\t11448\t9133\t11522\t6616\n+ENSMUSG00000087264\t1\t0\t0\t0\t0\t0\n+ENSMUSG00000068099\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000087651\t1\t1\t0\t0\t0\t0\n+ENSMUSG00000072694\t58\t13\t20\t26\t15\t17\n+ENSMUSG00000051319\t186\t177\t170\t164\t149\t229\n+ENSMUSG00000074578\t356\t143\t234\t122\t179\t142\n+ENSMUSG00000097247\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000094732\t0\t2\t1\t0\t2\t0\n+ENSMUSG00000026051\t1\t0\t0\t1\t1\t1\n+ENSMUSG00000087211\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000087579\t62\t81\t134\t28\t343\t91\n+ENSMUSG00000086899\t1\t0\t0\t0\t0\t0\n+ENSMUSG00000024118\t1225\t363\t383\t459\t2160\t664\n+ENSMUSG00000035595\t61\t46\t59\t59\t54\t39\n+ENSMUSG00000050088\t51\t20\t29\t24\t37\t42\n+ENSMUSG00000054676\t321\t178\t294\t187\t249\t260\n+ENSMUSG00000094690\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000035491\t0\t1\t1\t0\t0\t0\n+ENSMUSG00000029813\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000025500\t54\t15\t39\t44\t28\t36\n+ENSMUSG00000097048\t22\t17\t11\t12\t23\t23\n+ENSMUSG00000085114\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000032068\t95\t33\t60\t58\t83\t32\n+ENSMUSG00000046242\t1\t0\t1\t1\t1\t0\n+ENSMUSG00000057818\t0\t0\t1\t0\t1\t1\n+ENSMUSG00000029182\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000047150\t71\t57\t82\t58\t79\t34\n+ENSMUSG00000095384\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000085984\t1\t0\t0\t0\t0\t0\n+ENSMUSG00000044744\t1\t1\t1\t1\t1\t4\n+ENSMUSG00000095040\t1\t0\t0\t0\t0\t0\n+ENSMUSG00000056508\t2\t0\t1\t1\t2\t0\n+ENSMUSG00000075511\t369\t208\t462\t271\t474\t273\n+ENSMUSG00000021534\t2\t0\t1\t0\t1\t0\n+ENSMUSG00000044320\t2\t0\t0\t0\t1\t0\n+ENSMUSG00000018543\t0\t0\t0\t0\t0\t1\n+ENSMUSG0'..b'00000013668\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000067848\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000044084\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000079460\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000085122\t1\t0\t2\t0\t1\t0\n+ENSMUSG00000097908\t74\t144\t79\t100\t101\t47\n+ENSMUSG00000097852\t2\t1\t4\t3\t3\t2\n+ENSMUSG00000086209\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000013158\t0\t1\t1\t1\t2\t0\n+ENSMUSG00000084234\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000097494\t4\t3\t10\t1\t0\t6\n+ENSMUSG00000085796\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000090585\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000086025\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000087475\t9\t7\t6\t3\t9\t5\n+ENSMUSG00000027229\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000086690\t1\t0\t0\t0\t0\t0\n+ENSMUSG00000085820\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000050526\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000073759\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000087396\t46\t21\t25\t26\t24\t29\n+ENSMUSG00000026224\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000049357\t2\t3\t7\t4\t3\t8\n+ENSMUSG00000097407\t6\t2\t7\t5\t2\t2\n+ENSMUSG00000084848\t5\t2\t3\t3\t3\t2\n+ENSMUSG00000053896\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000062224\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000029586\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000090369\t108\t83\t112\t75\t91\t37\n+ENSMUSG00000031631\t435\t262\t313\t297\t409\t306\n+ENSMUSG00000071749\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000079304\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000021874\t7\t0\t0\t1\t0\t1\n+ENSMUSG00000079048\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000073730\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000074734\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000085198\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000071763\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000097262\t1\t0\t0\t0\t0\t0\n+ENSMUSG00000021415\t0\t31\t22\t22\t0\t0\n+ENSMUSG00000097637\t40\t28\t38\t30\t50\t21\n+ENSMUSG00000030463\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000085037\t119\t51\t74\t59\t96\t96\n+ENSMUSG00000070332\t32\t49\t68\t42\t42\t33\n+ENSMUSG00000085785\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000021718\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000086446\t0\t0\t1\t0\t0\t0\n+ENSMUSG00000021133\t27\t11\t17\t15\t22\t18\n+ENSMUSG00000055403\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000020807\t193\t75\t109\t79\t79\t117\n+ENSMUSG00000044276\t1\t1\t2\t1\t1\t1\n+ENSMUSG00000033015\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000030877\t2\t1\t0\t1\t0\t1\n+ENSMUSG00000085928\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000086355\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000058046\t2\t4\t3\t2\t3\t5\n+ENSMUSG00000086892\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000086968\t102\t43\t78\t96\t70\t142\n+ENSMUSG00000036770\t1\t0\t0\t0\t0\t0\n+ENSMUSG00000097044\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000097021\t4\t1\t0\t2\t0\t2\n+ENSMUSG00000027942\t311\t193\t281\t250\t336\t312\n+ENSMUSG00000035057\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000031736\t2\t4\t4\t3\t1\t0\n+ENSMUSG00000025288\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000093617\t6\t1\t0\t0\t5\t7\n+ENSMUSG00000072893\t73\t44\t47\t40\t57\t46\n+ENSMUSG00000085184\t4\t3\t1\t0\t2\t1\n+ENSMUSG00000097082\t1\t0\t4\t0\t3\t4\n+ENSMUSG00000045928\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000051940\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000034959\t64\t16\t36\t30\t13\t18\n+ENSMUSG00000073147\t217\t93\t166\t118\t178\t192\n+ENSMUSG00000085129\t9\t4\t4\t11\t6\t6\n+ENSMUSG00000044574\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000097885\t6\t1\t2\t2\t2\t0\n+ENSMUSG00000036046\t1374\t858\t1055\t838\t917\t1196\n+ENSMUSG00000055137\t558\t652\t687\t315\t846\t515\n+ENSMUSG00000087455\t8\t0\t2\t3\t0\t1\n+ENSMUSG00000086822\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000040412\t51\t14\t52\t23\t18\t16\n+ENSMUSG00000022639\t41\t9\t13\t19\t24\t18\n+ENSMUSG00000087620\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000094113\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000044121\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000085776\t1\t0\t0\t0\t1\t2\n+ENSMUSG00000085741\t2\t3\t2\t3\t1\t1\n+ENSMUSG00000097772\t73\t47\t49\t27\t79\t45\n+ENSMUSG00000028475\t20\t6\t14\t8\t7\t13\n+ENSMUSG00000074918\t173\t111\t115\t170\t185\t175\n+ENSMUSG00000006204\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000067613\t0\t0\t1\t1\t0\t0\n+ENSMUSG00000021566\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000085004\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000025058\t87\t10\t32\t22\t18\t21\n+ENSMUSG00000052688\t155\t30\t54\t73\t33\t63\n+ENSMUSG00000055594\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000020703\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000087174\t82\t44\t68\t45\t60\t46\n+ENSMUSG00000087424\t10\t10\t11\t9\t9\t14\n+ENSMUSG00000073755\t2\t0\t2\t1\t1\t2\n+ENSMUSG00000087553\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000081793\t0\t0\t1\t0\t0\t0\n+ENSMUSG00000097766\t0\t0\t0\t0\t0\t1\n+ENSMUSG00000097191\t0\t0\t0\t0\t0\t1\n+ENSMUSG00000091370\t2\t4\t1\t3\t7\t2\n+ENSMUSG00000057181\t540\t306\t400\t355\t504\t467\n+ENSMUSG00000086126\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000062046\t0\t0\t0\t2\t0\t0\n+ENSMUSG00000093459\t40\t16\t25\t12\t28\t13\n+ENSMUSG00000073197\t0\t0\t0\t0\t0\t0\n+ENSMUSG00000050549\t110\t55\t68\t54\t72\t95\n+ENSMUSG00000073101\t6\t2\t3\t3\t6\t7\n+ENSMUSG00000041605\t14\t3\t2\t1\t2\t13\n+ENSMUSG00000062470\t17\t6\t16\t9\t13\t13\n+ENSMUSG00000045968\t4\t3\t5\t3\t5\t63\n+ENSMUSG00000086259\t0\t0\t0\t0\t2\t0\n+ENSMUSG00000054672\t2\t0\t0\t0\t0\t0\n'
b
diff -r 000000000000 -r ac607906f10a test-data/per_base_coverage_default.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/per_base_coverage_default.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,299 @@
+#chr pos coverage
+test_chromosome 53 1
+test_chromosome 54 1
+test_chromosome 55 2
+test_chromosome 56 3
+test_chromosome 57 3
+test_chromosome 58 4
+test_chromosome 59 4
+test_chromosome 60 4
+test_chromosome 61 5
+test_chromosome 62 5
+test_chromosome 63 5
+test_chromosome 64 5
+test_chromosome 65 5
+test_chromosome 66 5
+test_chromosome 67 5
+test_chromosome 68 5
+test_chromosome 69 5
+test_chromosome 70 5
+test_chromosome 71 5
+test_chromosome 72 6
+test_chromosome 73 6
+test_chromosome 74 6
+test_chromosome 75 6
+test_chromosome 76 6
+test_chromosome 77 6
+test_chromosome 78 8
+test_chromosome 79 8
+test_chromosome 80 8
+test_chromosome 81 8
+test_chromosome 82 8
+test_chromosome 83 9
+test_chromosome 84 9
+test_chromosome 85 9
+test_chromosome 86 10
+test_chromosome 87 10
+test_chromosome 88 11
+test_chromosome 89 12
+test_chromosome 90 12
+test_chromosome 91 13
+test_chromosome 92 14
+test_chromosome 93 14
+test_chromosome 94 16
+test_chromosome 95 16
+test_chromosome 96 17
+test_chromosome 97 17
+test_chromosome 98 18
+test_chromosome 99 18
+test_chromosome 100 18
+test_chromosome 101 21
+test_chromosome 102 22
+test_chromosome 103 23
+test_chromosome 104 23
+test_chromosome 105 23
+test_chromosome 106 24
+test_chromosome 107 25
+test_chromosome 108 25
+test_chromosome 109 25
+test_chromosome 110 25
+test_chromosome 111 25
+test_chromosome 112 25
+test_chromosome 113 25
+test_chromosome 114 25
+test_chromosome 115 26
+test_chromosome 116 26
+test_chromosome 117 26
+test_chromosome 118 26
+test_chromosome 119 26
+test_chromosome 120 26
+test_chromosome 121 26
+test_chromosome 122 27
+test_chromosome 123 29
+test_chromosome 124 29
+test_chromosome 125 32
+test_chromosome 126 32
+test_chromosome 127 33
+test_chromosome 128 33
+test_chromosome 129 33
+test_chromosome 130 32
+test_chromosome 131 32
+test_chromosome 132 33
+test_chromosome 133 32
+test_chromosome 134 32
+test_chromosome 135 34
+test_chromosome 136 33
+test_chromosome 137 33
+test_chromosome 138 33
+test_chromosome 139 34
+test_chromosome 140 34
+test_chromosome 141 35
+test_chromosome 142 37
+test_chromosome 143 37
+test_chromosome 144 38
+test_chromosome 145 38
+test_chromosome 146 38
+test_chromosome 147 38
+test_chromosome 148 38
+test_chromosome 149 38
+test_chromosome 150 38
+test_chromosome 151 38
+test_chromosome 152 38
+test_chromosome 153 36
+test_chromosome 154 36
+test_chromosome 155 37
+test_chromosome 156 37
+test_chromosome 157 37
+test_chromosome 158 36
+test_chromosome 159 36
+test_chromosome 160 37
+test_chromosome 161 36
+test_chromosome 162 37
+test_chromosome 163 36
+test_chromosome 164 35
+test_chromosome 165 37
+test_chromosome 166 36
+test_chromosome 167 36
+test_chromosome 168 36
+test_chromosome 169 34
+test_chromosome 170 36
+test_chromosome 171 36
+test_chromosome 172 36
+test_chromosome 173 35
+test_chromosome 174 35
+test_chromosome 175 36
+test_chromosome 176 33
+test_chromosome 177 32
+test_chromosome 178 32
+test_chromosome 179 32
+test_chromosome 180 32
+test_chromosome 181 32
+test_chromosome 182 31
+test_chromosome 183 31
+test_chromosome 184 32
+test_chromosome 185 32
+test_chromosome 186 32
+test_chromosome 187 32
+test_chromosome 188 32
+test_chromosome 189 32
+test_chromosome 190 31
+test_chromosome 191 31
+test_chromosome 192 31
+test_chromosome 193 31
+test_chromosome 194 32
+test_chromosome 195 32
+test_chromosome 196 32
+test_chromosome 197 32
+test_chromosome 198 30
+test_chromosome 199 29
+test_chromosome 200 28
+test_chromosome 201 31
+test_chromosome 202 30
+test_chromosome 203 29
+test_chromosome 204 30
+test_chromosome 205 31
+test_chromosome 206 31
+test_chromosome 207 31
+test_chromosome 208 32
+test_chromosome 209 31
+test_chromosome 210 31
+test_chromosome 211 31
+test_chromosome 212 31
+test_chromosome 213 31
+test_chromosome 214 31
+test_chromosome 215 31
+test_chromosome 216 30
+test_chromosome 217 28
+test_chromosome 218 28
+test_chromosome 219 27
+test_chromosome 220 27
+test_chromosome 221 28
+test_chromosome 222 28
+test_chromosome 223 28
+test_chromosome 224 28
+test_chromosome 225 29
+test_chromosome 226 29
+test_chromosome 227 30
+test_chromosome 228 30
+test_chromosome 229 31
+test_chromosome 230 30
+test_chromosome 231 31
+test_chromosome 232 32
+test_chromosome 233 33
+test_chromosome 234 34
+test_chromosome 235 33
+test_chromosome 236 33
+test_chromosome 237 32
+test_chromosome 238 32
+test_chromosome 239 32
+test_chromosome 240 30
+test_chromosome 241 30
+test_chromosome 242 29
+test_chromosome 243 31
+test_chromosome 244 32
+test_chromosome 245 30
+test_chromosome 246 29
+test_chromosome 247 29
+test_chromosome 248 29
+test_chromosome 249 29
+test_chromosome 250 28
+test_chromosome 351 31
+test_chromosome 352 32
+test_chromosome 353 33
+test_chromosome 354 34
+test_chromosome 355 35
+test_chromosome 356 34
+test_chromosome 357 34
+test_chromosome 358 35
+test_chromosome 359 34
+test_chromosome 360 35
+test_chromosome 361 35
+test_chromosome 362 36
+test_chromosome 363 36
+test_chromosome 364 36
+test_chromosome 365 36
+test_chromosome 366 37
+test_chromosome 367 37
+test_chromosome 368 37
+test_chromosome 369 37
+test_chromosome 370 36
+test_chromosome 371 37
+test_chromosome 372 37
+test_chromosome 373 39
+test_chromosome 374 39
+test_chromosome 375 38
+test_chromosome 376 36
+test_chromosome 377 36
+test_chromosome 378 36
+test_chromosome 379 36
+test_chromosome 380 34
+test_chromosome 381 33
+test_chromosome 382 32
+test_chromosome 383 32
+test_chromosome 384 31
+test_chromosome 385 31
+test_chromosome 386 31
+test_chromosome 387 31
+test_chromosome 388 31
+test_chromosome 389 30
+test_chromosome 390 30
+test_chromosome 391 30
+test_chromosome 392 30
+test_chromosome 393 30
+test_chromosome 394 30
+test_chromosome 395 30
+test_chromosome 396 29
+test_chromosome 397 28
+test_chromosome 398 28
+test_chromosome 399 28
+test_chromosome 400 27
+test_chromosome 501 25
+test_chromosome 502 25
+test_chromosome 503 25
+test_chromosome 504 25
+test_chromosome 505 25
+test_chromosome 506 24
+test_chromosome 507 23
+test_chromosome 508 22
+test_chromosome 509 21
+test_chromosome 510 21
+test_chromosome 511 21
+test_chromosome 512 21
+test_chromosome 513 21
+test_chromosome 514 21
+test_chromosome 515 21
+test_chromosome 516 21
+test_chromosome 517 20
+test_chromosome 518 18
+test_chromosome 519 17
+test_chromosome 520 17
+test_chromosome 521 17
+test_chromosome 522 16
+test_chromosome 523 15
+test_chromosome 524 15
+test_chromosome 525 14
+test_chromosome 526 14
+test_chromosome 527 13
+test_chromosome 528 12
+test_chromosome 529 11
+test_chromosome 530 10
+test_chromosome 531 10
+test_chromosome 532 10
+test_chromosome 533 9
+test_chromosome 534 9
+test_chromosome 535 8
+test_chromosome 536 8
+test_chromosome 537 7
+test_chromosome 538 7
+test_chromosome 539 7
+test_chromosome 540 7
+test_chromosome 541 6
+test_chromosome 542 6
+test_chromosome 543 6
+test_chromosome 544 5
+test_chromosome 545 5
+test_chromosome 546 4
+test_chromosome 547 4
+test_chromosome 548 2
+test_chromosome 549 2
+test_chromosome 550 1
b
diff -r 000000000000 -r ac607906f10a test-data/per_base_coverage_inside_features.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/per_base_coverage_inside_features.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,199 @@
+#chr pos coverage
+test_chromosome 53 1
+test_chromosome 54 1
+test_chromosome 55 2
+test_chromosome 56 3
+test_chromosome 57 3
+test_chromosome 58 4
+test_chromosome 59 4
+test_chromosome 60 4
+test_chromosome 61 5
+test_chromosome 62 5
+test_chromosome 63 5
+test_chromosome 64 5
+test_chromosome 65 5
+test_chromosome 66 5
+test_chromosome 67 5
+test_chromosome 68 5
+test_chromosome 69 5
+test_chromosome 70 5
+test_chromosome 71 5
+test_chromosome 72 6
+test_chromosome 73 6
+test_chromosome 74 6
+test_chromosome 75 6
+test_chromosome 76 6
+test_chromosome 77 6
+test_chromosome 78 8
+test_chromosome 79 8
+test_chromosome 80 8
+test_chromosome 81 8
+test_chromosome 82 8
+test_chromosome 83 9
+test_chromosome 84 9
+test_chromosome 85 9
+test_chromosome 86 10
+test_chromosome 87 10
+test_chromosome 88 11
+test_chromosome 89 12
+test_chromosome 90 12
+test_chromosome 91 13
+test_chromosome 92 14
+test_chromosome 93 14
+test_chromosome 94 16
+test_chromosome 95 16
+test_chromosome 96 17
+test_chromosome 97 17
+test_chromosome 98 18
+test_chromosome 99 18
+test_chromosome 100 18
+test_chromosome 101 21
+test_chromosome 102 22
+test_chromosome 103 23
+test_chromosome 104 23
+test_chromosome 105 23
+test_chromosome 106 24
+test_chromosome 107 25
+test_chromosome 108 25
+test_chromosome 109 25
+test_chromosome 110 25
+test_chromosome 111 25
+test_chromosome 112 25
+test_chromosome 113 25
+test_chromosome 114 25
+test_chromosome 115 26
+test_chromosome 116 26
+test_chromosome 117 26
+test_chromosome 118 26
+test_chromosome 119 26
+test_chromosome 120 26
+test_chromosome 121 26
+test_chromosome 122 27
+test_chromosome 123 29
+test_chromosome 124 29
+test_chromosome 125 32
+test_chromosome 126 32
+test_chromosome 127 33
+test_chromosome 128 33
+test_chromosome 129 33
+test_chromosome 130 32
+test_chromosome 131 32
+test_chromosome 132 33
+test_chromosome 133 32
+test_chromosome 134 32
+test_chromosome 135 34
+test_chromosome 136 33
+test_chromosome 137 33
+test_chromosome 138 33
+test_chromosome 139 34
+test_chromosome 140 34
+test_chromosome 141 35
+test_chromosome 142 37
+test_chromosome 143 37
+test_chromosome 144 38
+test_chromosome 145 38
+test_chromosome 146 38
+test_chromosome 147 38
+test_chromosome 148 38
+test_chromosome 149 38
+test_chromosome 150 38
+test_chromosome 151 38
+test_chromosome 152 38
+test_chromosome 153 36
+test_chromosome 154 36
+test_chromosome 155 37
+test_chromosome 156 37
+test_chromosome 157 37
+test_chromosome 158 36
+test_chromosome 159 36
+test_chromosome 160 37
+test_chromosome 161 36
+test_chromosome 162 37
+test_chromosome 163 36
+test_chromosome 164 35
+test_chromosome 165 37
+test_chromosome 166 36
+test_chromosome 167 36
+test_chromosome 168 36
+test_chromosome 169 34
+test_chromosome 170 36
+test_chromosome 171 36
+test_chromosome 172 36
+test_chromosome 173 35
+test_chromosome 174 35
+test_chromosome 175 36
+test_chromosome 176 33
+test_chromosome 177 32
+test_chromosome 178 32
+test_chromosome 179 32
+test_chromosome 180 32
+test_chromosome 181 32
+test_chromosome 182 31
+test_chromosome 183 31
+test_chromosome 184 32
+test_chromosome 185 32
+test_chromosome 186 32
+test_chromosome 187 32
+test_chromosome 188 32
+test_chromosome 189 32
+test_chromosome 190 31
+test_chromosome 191 31
+test_chromosome 192 31
+test_chromosome 193 31
+test_chromosome 194 32
+test_chromosome 195 32
+test_chromosome 196 32
+test_chromosome 197 32
+test_chromosome 198 30
+test_chromosome 199 29
+test_chromosome 200 28
+test_chromosome 201 31
+test_chromosome 202 30
+test_chromosome 203 29
+test_chromosome 204 30
+test_chromosome 205 31
+test_chromosome 206 31
+test_chromosome 207 31
+test_chromosome 208 32
+test_chromosome 209 31
+test_chromosome 210 31
+test_chromosome 211 31
+test_chromosome 212 31
+test_chromosome 213 31
+test_chromosome 214 31
+test_chromosome 215 31
+test_chromosome 216 30
+test_chromosome 217 28
+test_chromosome 218 28
+test_chromosome 219 27
+test_chromosome 220 27
+test_chromosome 221 28
+test_chromosome 222 28
+test_chromosome 223 28
+test_chromosome 224 28
+test_chromosome 225 29
+test_chromosome 226 29
+test_chromosome 227 30
+test_chromosome 228 30
+test_chromosome 229 31
+test_chromosome 230 30
+test_chromosome 231 31
+test_chromosome 232 32
+test_chromosome 233 33
+test_chromosome 234 34
+test_chromosome 235 33
+test_chromosome 236 33
+test_chromosome 237 32
+test_chromosome 238 32
+test_chromosome 239 32
+test_chromosome 240 30
+test_chromosome 241 30
+test_chromosome 242 29
+test_chromosome 243 31
+test_chromosome 244 32
+test_chromosome 245 30
+test_chromosome 246 29
+test_chromosome 247 29
+test_chromosome 248 29
+test_chromosome 249 29
+test_chromosome 250 28
b
diff -r 000000000000 -r ac607906f10a test-data/per_base_coverage_outside_features.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/per_base_coverage_outside_features.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,101 @@
+#chr pos coverage
+test_chromosome 351 31
+test_chromosome 352 32
+test_chromosome 353 33
+test_chromosome 354 34
+test_chromosome 355 35
+test_chromosome 356 34
+test_chromosome 357 34
+test_chromosome 358 35
+test_chromosome 359 34
+test_chromosome 360 35
+test_chromosome 361 35
+test_chromosome 362 36
+test_chromosome 363 36
+test_chromosome 364 36
+test_chromosome 365 36
+test_chromosome 366 37
+test_chromosome 367 37
+test_chromosome 368 37
+test_chromosome 369 37
+test_chromosome 370 36
+test_chromosome 371 37
+test_chromosome 372 37
+test_chromosome 373 39
+test_chromosome 374 39
+test_chromosome 375 38
+test_chromosome 376 36
+test_chromosome 377 36
+test_chromosome 378 36
+test_chromosome 379 36
+test_chromosome 380 34
+test_chromosome 381 33
+test_chromosome 382 32
+test_chromosome 383 32
+test_chromosome 384 31
+test_chromosome 385 31
+test_chromosome 386 31
+test_chromosome 387 31
+test_chromosome 388 31
+test_chromosome 389 30
+test_chromosome 390 30
+test_chromosome 391 30
+test_chromosome 392 30
+test_chromosome 393 30
+test_chromosome 394 30
+test_chromosome 395 30
+test_chromosome 396 29
+test_chromosome 397 28
+test_chromosome 398 28
+test_chromosome 399 28
+test_chromosome 400 27
+test_chromosome 501 25
+test_chromosome 502 25
+test_chromosome 503 25
+test_chromosome 504 25
+test_chromosome 505 25
+test_chromosome 506 24
+test_chromosome 507 23
+test_chromosome 508 22
+test_chromosome 509 21
+test_chromosome 510 21
+test_chromosome 511 21
+test_chromosome 512 21
+test_chromosome 513 21
+test_chromosome 514 21
+test_chromosome 515 21
+test_chromosome 516 21
+test_chromosome 517 20
+test_chromosome 518 18
+test_chromosome 519 17
+test_chromosome 520 17
+test_chromosome 521 17
+test_chromosome 522 16
+test_chromosome 523 15
+test_chromosome 524 15
+test_chromosome 525 14
+test_chromosome 526 14
+test_chromosome 527 13
+test_chromosome 528 12
+test_chromosome 529 11
+test_chromosome 530 10
+test_chromosome 531 10
+test_chromosome 532 10
+test_chromosome 533 9
+test_chromosome 534 9
+test_chromosome 535 8
+test_chromosome 536 8
+test_chromosome 537 7
+test_chromosome 538 7
+test_chromosome 539 7
+test_chromosome 540 7
+test_chromosome 541 6
+test_chromosome 542 6
+test_chromosome 543 6
+test_chromosome 544 5
+test_chromosome 545 5
+test_chromosome 546 4
+test_chromosome 547 4
+test_chromosome 548 2
+test_chromosome 549 2
+test_chromosome 550 1
b
diff -r 000000000000 -r ac607906f10a test-data/rnaseq_qc_counts_custom.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaseq_qc_counts_custom.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,2 @@
+#GeneID try_this
+GENE1 45
b
diff -r 000000000000 -r ac607906f10a test-data/rnaseq_qc_results_custom.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaseq_qc_results_custom.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,48 @@
+RNA-Seq QC report
+-----------------------------------
+
+>>>>>>> Input
+
+    bam file = /tmp/tmp_CFvEQ/files/000/dataset_49.dat
+    gff file = /tmp/tmp_CFvEQ/files/000/dataset_2.dat
+    counting algorithm = proportional
+    protocol = strand-specific-forward
+
+
+>>>>>>> Reads alignment
+
+    reads aligned (left/right) = 0 / 0
+    read pairs aligned  = 0
+    total alignments = 101
+    secondary alignments = 1
+    non-unique alignments = 0
+    aligned to genes  = 46
+    ambiguous alignments = 0
+    no feature assigned = 55
+    not aligned = 0
+
+
+>>>>>>> Reads genomic origin
+
+    exonic =  46 (45.54%)
+    intronic = 0 (0%)
+    intergenic = 55 (54.46%)
+    overlapping exon = 0 (0%)
+
+
+>>>>>>> Transcript coverage profile
+
+    5' bias = NaN
+    3' bias = NaN
+    5'-3' bias = NaN
+
+
+>>>>>>> Junction analysis
+
+    reads at junctions = 53
+
+    CTTT : 47.17%
+    CCAC : 47.17%
+    GCAC : 1.89%
+    CTAC : 1.89%
+    GAAT : 1.89%
b
diff -r 000000000000 -r ac607906f10a test-data/rnaseq_qc_results_default.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rnaseq_qc_results_default.txt Thu Oct 10 17:42:04 2019 -0400
b
@@ -0,0 +1,47 @@
+RNA-Seq QC report
+-----------------------------------
+
+>>>>>>> Input
+
+    bam file = /tmp/tmp_CFvEQ/files/000/dataset_1.dat
+    gff file = /tmp/tmp_CFvEQ/files/000/dataset_2.dat
+    counting algorithm = uniquely-mapped-reads
+    protocol = non-strand-specific
+
+
+>>>>>>> Reads alignment
+
+    reads aligned  = 100
+    total alignments = 101
+    secondary alignments = 1
+    non-unique alignments = 2
+    aligned to genes  = 54
+    ambiguous alignments = 0
+    no feature assigned = 45
+    not aligned = 0
+    SSP estimation (fwd/rev) = 0.59 / 0.41
+
+
+>>>>>>> Reads genomic origin
+
+    exonic =  54 (54.55%)
+    intronic = 0 (0%)
+    intergenic = 45 (45.45%)
+    overlapping exon = 0 (0%)
+
+
+>>>>>>> Transcript coverage profile
+
+    5' bias = NaN
+    3' bias = NaN
+    5'-3' bias = NaN
+
+
+>>>>>>> Junction analysis
+
+    reads at junctions = 52
+
+    CTTT : 48.08%
+    CCAC : 48.08%
+    GCAC : 1.92%
+    CTAC : 1.92%
b
diff -r 000000000000 -r ac607906f10a test-data/test_mapped_reads.bam
b
Binary file test-data/test_mapped_reads.bam has changed