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

Changeset 1:4a89c6f84425 (2019-10-29)
Previous changeset 0:ac607906f10a (2019-10-10) Next changeset 2:e7fd6754d093 (2020-01-13)
Commit message:
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/qualimap commit 08beeabf7a8e91d0e95f5cdb71249dd6efe2bd46"
modified:
qualimap_bamqc.xml
test-data/features.gtf
b
diff -r ac607906f10a -r 4a89c6f84425 qualimap_bamqc.xml
--- a/qualimap_bamqc.xml Thu Oct 10 17:42:04 2019 -0400
+++ b/qualimap_bamqc.xml Tue Oct 29 09:14:20 2019 -0400
[
@@ -1,4 +1,4 @@
-<tool id="qualimap_bamqc" name="QualiMap BamQC" version="@VERSION@">
+<tool id="qualimap_bamqc" name="QualiMap BamQC" version="@VERSION@+galaxy1">
     <macros>
         <import>qualimap_macros.xml</import>
     </macros>
@@ -21,11 +21,12 @@
         ## from $out_dir.
         #set $per_base_coverage_target = $coverage_file
 
-        qualimap bamqc
-        -bam '$input1' -outdir results -outformat html
-        --collect-overlap-pairs
         #if str($stats_regions.region_select) == 'custom_regions':
-            -gff ${stats_regions.regions}
+            ## Have qualimap detect file format by suffix
+            ## since its autodetection is unreliable.
+            #set $regions_file = '.'.join(['regions', str($stats_regions.regions.ext)])
+            ln -s '${stats_regions.regions}' ${regions_file} &&
+
             #if $stats_regions.outside_stats:
                 #set $report_name = 'qualimapReportOutsideRegions'
                 #set $summary_report = 'outside_results.txt'
@@ -35,9 +36,15 @@
                 #if $per_base_coverage:
                     #set $per_base_coverage_target = '/dev/null'
                 #end if
-                ${stats_regions.outside_stats}
             #end if
         #end if
+
+        qualimap bamqc
+        -bam '$input1' -outdir results -outformat html
+        --collect-overlap-pairs
+        #if str($stats_regions.region_select) == 'custom_regions':
+            -gff ${regions_file} ${stats_regions.outside_stats}
+        #end if
         #if $per_base_coverage:
             $per_base_coverage $per_base_coverage_target
         #end if
b
diff -r ac607906f10a -r 4a89c6f84425 test-data/features.gtf
--- a/test-data/features.gtf Thu Oct 10 17:42:04 2019 -0400
+++ b/test-data/features.gtf Tue Oct 29 09:14:20 2019 -0400
b
@@ -1,3 +1,4 @@
+#comment lines (like this one) break qualimap's regions file autodetection
 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";