Qualimap rnaseq (#74)

* first version * complete script for qualimap * add escaping character before leading hashtag (#50) * add escaping character before leading hashtag * update changelog * Update CHANGELOG.md Co-authored-by: Robrecht Cannoodt <[email protected]> * replace escaping \ by \\ --------- Co-authored-by: Robrecht Cannoodt <[email protected]> * Samtools collate (#49) * initial commit dedup * Revert "initial commit dedup" This reverts commit 38f586b. * Initial commit, whole component is functional * Update viash (#51) * update viash * update readme * update changelog * update changelog * fix incorrect heading detection * update again * clean up readme * Samtools view (#48) * initial commit dedup * Revert "initial commit dedup" This reverts commit 38f586b. * initial version with a few tests, script, and config file * update changelog, add one test * add a 4th test, fix option names in the script * Fix name of component in config * remove option named with a number * add must_exist to input file argument * removed "default: null" from one of the arguments in config * remove utf8 characters from config * Update CHANGELOG.md --------- Co-authored-by: Robrecht Cannoodt <[email protected]> * Samtools fastq (#52) * initial commit dedup * Revert "initial commit dedup" This reverts commit 38f586b. * Initial commit, config, script, help and test_data * Update changelog, add tests, fix argument naming errors, add test data * update changelog, remove gffread namespace field --------- Co-authored-by: Robrecht Cannoodt <[email protected]> * format URL in the description (#55) * format URL in the description * update changelog * Change name in _viash.yaml (#60) * Update operational code (#63) * update readme * switch ci to toolbox * update to viash 0.9.0-RC6 * edit keywords * fix version * update biobox * cutadapt (#7) * First commit, clone of cutadapt in htrnaseq + help.txt * Add config * Don't allow multiple: true when providing a FASTA file with adapters * First version of script * Updates and fixes - se/pe * Add tests and fix --json argument * Add software version * Better consistency in using snake_case * Update src/cutadapt/config.vsh.yaml Co-authored-by: Robrecht Cannoodt <[email protected]> * Update src/cutadapt/config.vsh.yaml Co-authored-by: Robrecht Cannoodt <[email protected]> * Update src/cutadapt/config.vsh.yaml Co-authored-by: Robrecht Cannoodt <[email protected]> * Specify --input and --input_r2 as separate arguments * Avoid specifying default arg values * Add more information to `--minimum_length` and `maximum_length` * Add --cpus by means of $meta_cpus and set proper default * Allow multiple for adapters/fasta and add test * change multiple_sep to ';' * add example * simplify code with a helper function * create directories in test * use a different output extension if --fasta is provided * decrease code duplication by separating optional outputs from paired/unpaired output arguments * write custom tests for cutadapt * fix _r2 arguments * add debug flag as not to always print the cli command * remove comment * Update to Viash 0.9.0-RC4 * Ability to specify output globbing patterns * Avoid the need for both output_dir and output * Move fields from `info` to `links` Co-authored-by: Robrecht Cannoodt <[email protected]> * Move references back to the info field * apologies, I proposed a wrong syntax --------- Co-authored-by: Robrecht Cannoodt <[email protected]> * update changelog * update readme * Update salmon quant arguments (#57) * Make index an optional argument * FIx argument type and add optional argument * FEAT: add bedtools getfasta. (#59) * FEAT: add bedtools getfasta. * Add PR number to CHANGELOG * Add star genomegenerate component (#58) * Add star genomegenerate component * Update changelog * Rename component * Update test * Update CHANGELOG.md --------- Co-authored-by: Robrecht Cannoodt <[email protected]> * fix package config (#65) * Delete src/bgzip directory (#64) It was moved to toolbox * Output alignments to the transcriptome (#56) * Output alignments to the transcriptome * Change argument name * BUG: pear component failure is ignored (#70) * FEAT + BUG: cutadapt; allowing disabling demultiplexing and fix par_quality_cutoff_r2 (#69) * FEAT: Disable cutadapt demultiplexing by default * Cutadapt: fix --par_quality_cutoff_r2 * FEAT: update busco to 5.7.1 (#72) * FEAT: update busco to 5.7.1 * Typo * Samtools fasta (#53) * initial commit dedup * Revert "initial commit dedup" This reverts commit 38f586b. * Fasta component * change script resource to samtools_fastq script, with dummy argument to specify the command * add dummy argument to samtools_fastq to share the script with samtools_fasta * fix path to script in config * Update src/samtools/samtools_fastq/script.sh Co-authored-by: Robrecht Cannoodt <[email protected]> * Change default fields to examples * Two more default fields changed to examples * Minor formatting changes * Markdown formatting changes in configs --------- Co-authored-by: Robrecht Cannoodt <[email protected]> * Umi tools dedup (#54) * initial commit dedup * Revert "initial commit dedup" This reverts commit 38f586b. * inital commit dedup * Working component with one test * Update test 1 and test data, fix some arg types in config and script * test data files and changes to script * Add third test and test data * Fix typo in script * remove utf8 characters in config * Add choices fields and change default fields to exampels * Minor formatting changes * md formatting changes in config * Fix typo (#79) * add vscode to gitignore * update multiple separator (#81) * update multiple separator * update changelog * Update src/multiqc/config.vsh.yaml Co-authored-by: Robrecht Cannoodt <[email protected]> * Update src/multiqc/config.vsh.yaml Co-authored-by: Robrecht Cannoodt <[email protected]> * Update src/multiqc/config.vsh.yaml Co-authored-by: Robrecht Cannoodt <[email protected]> * Update src/multiqc/config.vsh.yaml Co-authored-by: Robrecht Cannoodt <[email protected]> * update ifs --------- Co-authored-by: Robrecht Cannoodt <[email protected]> * add test data * add tests * update changelog * remove unrequired test data * update descriptions * update changelog * update help text * Update src/qualimap/qualimap_rnaseq/script.sh Co-authored-by: Robrecht Cannoodt <[email protected]> * update unit tests * update unit tests * addres pr changes request * add version * remove whitespace multiqc * Apply suggestions from code review Co-authored-by: Robrecht Cannoodt <[email protected]> * address pr comments * Update CHANGELOG.md * fix doi * Fix name * update version and container image * write software version to file --------- Co-authored-by: dorien-er <[email protected]> Co-authored-by: Leila011 <[email protected]> Co-authored-by: Robrecht Cannoodt <[email protected]> Co-authored-by: emmarousseau <[email protected]> Co-authored-by: Sai Nirmayi Yasa <[email protected]> Co-authored-by: Dries Schaumont <[email protected]> Co-authored-by: Dorien <[email protected]>
viash-hub · Aug 21, 2024 · 766ab6c · 766ab6c
1 parent 923a6da
commit 766ab6c
Show file tree

Hide file tree

Showing 8 changed files with 339 additions and 0 deletions.
diff --git a/CHANGELOG.md b/CHANGELOG.md
@@ -31,6 +31,8 @@
   - `bedtools/bedtools_sort`: Sorts a feature file (bed/gff/vcf) by chromosome and other criteria (PR #98).
   - `bedtools/bedtools_bamtofastq`: Convert BAM alignments to FASTQ files (PR #101).
   - `bedtools/bedtools_bedtobam`: Converts genomic feature records (bed/gff/vcf) to BAM format (PR #111).
+
+* `qualimap/qualimap_rnaseq`: RNA-seq QC analysis using qualimap (PR #74). 
 
 ## MINOR CHANGES
 

diff --git a/src/qualimap/qualimap_rnaseq/config.vsh.yaml b/src/qualimap/qualimap_rnaseq/config.vsh.yaml
@@ -0,0 +1,103 @@
+name: qualimap_rnaseq
+namespace: qualimap
+keywords: [RNA-seq, quality control, QC Report]
+description: |
+  Qualimap RNA-seq QC reports quality control metrics and bias estimations 
+  which are specific for whole transcriptome sequencing, including reads genomic 
+  origin, junction analysis, transcript coverage and 5’-3’ bias computation.
+links:
+  homepage: http://qualimap.conesalab.org/
+  documentation: http://qualimap.conesalab.org/doc_html/analysis.html#rna-seq-qc
+  issue_tracker: https://bitbucket.org/kokonech/qualimap/issues?status=new&status=open
+  repository: https://bitbucket.org/kokonech/qualimap/commits/branch/master
+references:
+  doi: 10.1093/bioinformatics/btv566
+license: GPL-2.0
+authors:
+  - __merge__: /src/_authors/dorien_roosen.yaml
+    roles: [ author, maintainer ]
+argument_groups:
+  - name: "Input"
+    arguments: 
+    - name: "--bam"
+      type: file
+      required: true
+      example: alignment.bam
+      description: Path to the sequence alignment file in BAM format, produced by a splicing-aware aligner.
+    - name: "--gtf"
+      type: file
+      required: true
+      example: annotations.gtf
+      description: Path to genomic annotations in Ensembl GTF format.
+
+  - name: "Output"
+    arguments: 
+    - name: "--qc_results"
+      direction: output
+      type: file
+      required: true
+      example: rnaseq_qc_results.txt
+      description: Text file containing the RNAseq QC results.
+    - name: "--counts"
+      type: file
+      required: false
+      direction: output
+      description: Output file for computed counts.
+    - name: "--report"
+      type: file
+      direction: output
+      required: false
+      example: report.html
+      description: Report output file. Supported formats are PDF or HTML.
+
+  - name: "Optional"
+    arguments: 
+    - name: "--num_pr_bases"
+      type: integer
+      required: false
+      min: 1
+      description: Number of upstream/downstream nucleotide bases to compute 5'-3' bias (default = 100).
+    - name: "--num_tr_bias"
+      type: integer
+      required: false
+      min: 1
+      description: Number of top highly expressed transcripts to compute 5'-3' bias (default = 1000).
+    - name: "--algorithm"
+      type: string
+      required: false
+      choices: ["uniquely-mapped-reads", "proportional"]
+      description: Counting algorithm (uniquely-mapped-reads (default) or proportional).
+    - name: "--sequencing_protocol"
+      type: string
+      required: false
+      choices: ["non-strand-specific", "strand-specific-reverse", "strand-specific-forward"]
+      description: Sequencing library protocol (strand-specific-forward, strand-specific-reverse or non-strand-specific (default)).
+    - name: "--paired"
+      type: boolean_true
+      description: Setting this flag for paired-end experiments will result in counting fragments instead of reads.
+    - name: "--sorted"
+      type: boolean_true
+      description: Setting this flag indicates that the input file is already sorted by name. If flag is not set, additional sorting by name will be performed. Only requiredfor paired-end analysis.
+    - name: "--java_memory_size"
+      type: string
+      required: false
+      description: maximum Java heap memory size, default = 4G.
+
+resources:
+  - type: bash_script
+    path: script.sh
+test_resources:
+  - type: bash_script
+    path: test.sh
+  - path: test_data/
+
+engines:
+  - type: docker
+    image: quay.io/biocontainers/qualimap:2.3--hdfd78af_0
+    setup:   
+      - type: docker
+        run: |
+          echo QualiMap: $(qualimap 2>&1 | grep QualiMap | sed 's/^.*QualiMap//') > /var/software_versions.txt
+runners: 
+  - type: executable
+  - type: nextflow
diff --git a/src/qualimap/qualimap_rnaseq/help.txt b/src/qualimap/qualimap_rnaseq/help.txt
@@ -0,0 +1,52 @@
+QualiMap v.2.3
+Built on 2023-05-19 16:57
+
+usage: qualimap <tool> [options]
+
+To launch GUI leave <tool> empty.
+
+Available tools:
+
+    bamqc            Evaluate NGS mapping to a reference genome
+    rnaseq           Evaluate RNA-seq alignment data
+    counts           Counts data analysis (further RNA-seq data evaluation)
+    multi-bamqc      Compare QC reports from multiple NGS mappings
+    clustering       Cluster epigenomic signals
+    comp-counts      Compute feature counts
+
+Special arguments: 
+
+    --java-mem-size  Use this argument to set Java memory heap size. Example:
+                     qualimap bamqc -bam very_large_alignment.bam --java-mem-size=4G
+
+usage: qualimap rnaseq [-a <arg>] -bam <arg> -gtf <arg> [-npb <arg>] [-ntb
+       <arg>] [-oc <arg>] [-outdir <arg>] [-outfile <arg>] [-outformat <arg>]
+       [-p <arg>] [-pe] [-s]
+ -a,--algorithm <arg>             Counting algorithm:
+                                  uniquely-mapped-reads(default) or
+                                  proportional.
+ -bam <arg>                       Input mapping file in BAM format.
+ -gtf <arg>                       Annotations file in Ensembl GTF format.
+ -npb,--num-pr-bases <arg>        Number of upstream/downstream nucleotide bases
+                                  to compute 5'-3' bias (default is 100).
+ -ntb,--num-tr-bias <arg>         Number of top highly expressed transcripts to
+                                  compute 5'-3' bias (default is 1000).
+ -oc <arg>                        Output file for computed counts. If only name
+                                  of the file is provided, then the file will be
+                                  saved in the output folder.
+ -outdir <arg>                    Output folder for HTML report and raw data.
+ -outfile <arg>                   Output file for PDF report (default value is
+                                  report.pdf).
+ -outformat <arg>                 Format of the output report (PDF, HTML or both
+                                  PDF:HTML, default is HTML).
+ -p,--sequencing-protocol <arg>   Sequencing library protocol:
+                                  strand-specific-forward,
+                                  strand-specific-reverse or non-strand-specific
+                                  (default)
+ -pe,--paired                     Setting this flag for paired-end experiments
+                                  will result in counting fragments instead of
+                                  reads
+ -s,--sorted                      This flag indicates that the input file is
+                                  already sorted by name. If not set, additional
+                                  sorting by name will be performed. Only
+                                  required for paired-end analysis.
diff --git a/src/qualimap/qualimap_rnaseq/script.sh b/src/qualimap/qualimap_rnaseq/script.sh
@@ -0,0 +1,50 @@
+#!/bin/bash
+
+set -eo pipefail
+
+tmp_dir=$(mktemp -d -p "$meta_temp_dir" qualimap_XXXXXXXXX)
+
+# Handle output parameters
+if [ -n "$par_report" ]; then
+    outfile=$(basename "$par_report")
+    report_extension="${outfile##*.}"
+fi
+
+if [ -n "$par_counts" ]; then
+    counts=$(basename "$par_counts")
+fi
+
+# disable flags
+[[ "$par_paired" == "false" ]] && unset par_paired
+[[ "$par_sorted" == "false" ]] && unset par_sorted
+
+# Run qualimap
+qualimap rnaseq \
+    ${meta_memory_mb:+--java-mem-size=${meta_memory_mb}M} \
+    ${par_algorithm:+--algorithm $par_algorithm} \
+    ${par_sequencing_protocol:+--sequencing-protocol $par_sequencing_protocol} \
+    -bam $par_bam \
+    -gtf $par_gtf \
+    -outdir "$tmp_dir" \
+    ${par_num_pr_bases:+--num-pr-bases $par_num_pr_bases} \
+    ${par_num_tr_bias:+--num-tr-bias $par_num_tr_bias} \
+    ${par_report:+-outformat $report_extension} \
+    ${par_paired:+--paired} \
+    ${par_sorted:+--sorted} \
+    ${par_report:+-outfile "$outfile"} \
+    ${par_counts:+-oc "$counts"}
+
+# Move output files
+mv "$tmp_dir/rnaseq_qc_results.txt" "$par_qc_results"
+
+if [ -n "$par_report" ] && [ $report_extension = "html" ]; then
+    mv "$tmp_dir/qualimapReport.html" "$par_report"
+fi
+
+if [ -n "$par_report" ] && [ $report_extension = "pdf" ]; then
+    mv "$tmp_dir/$outfile" "$par_report"
+fi
+
+if [ -n "$par_counts" ]; then
+    mv "$tmp_dir/$counts" "$par_counts"
+fi
diff --git a/src/qualimap/qualimap_rnaseq/test.sh b/src/qualimap/qualimap_rnaseq/test.sh
@@ -0,0 +1,112 @@
+set -e
+
+#############################################
+# helper functions
+assert_file_exists() {
+  [ -f "$1" ] || { echo "File '$1' does not exist" && exit 1; }
+}
+assert_file_doesnt_exist() {
+  [ ! -f "$1" ] || { echo "File '$1' exists but shouldn't" && exit 1; }
+}
+assert_file_not_empty() {
+  [ -s "$1" ] || { echo "File '$1' is empty but shouldn't be" && exit 1; }
+}
+assert_file_contains() {
+  grep -q "$2" "$1" || { echo "File '$1' does not contain '$2'" && exit 1; }
+}
+#############################################
+
+
+test_dir="$meta_resources_dir/test_data"
+
+mkdir "run_qualimap_rnaseq_html"
+cd "run_qualimap_rnaseq_html"
+
+echo "> Running qualimap with html output report"
+
+"$meta_executable" \
+    --bam $test_dir/a.bam \
+    --gtf $test_dir/annotation.gtf \
+    --report report.html \
+    --counts counts.txt \
+    --qc_results output.txt
+
+echo ">> Checking output"
+assert_file_exists "report.html"
+assert_file_exists "counts.txt"
+assert_file_exists "output.txt"
+assert_file_doesnt_exist "report.pdf"
+
+echo ">> Checking if output is empty"
+assert_file_not_empty "report.html"
+assert_file_not_empty "counts.txt"
+assert_file_not_empty "output.txt"
+
+echo ">> Checking output contents"
+assert_file_contains "output.txt" ">>>>>>> Input"
+assert_file_contains "output.txt" ">>>>>>> Reads alignment"
+assert_file_contains "output.txt" ">>>>>>> Reads genomic origin"
+assert_file_contains "output.txt" ">>>>>>> Transcript coverage profile"
+assert_file_contains "output.txt" ">>>>>>> Junction analysis"
+assert_file_contains "output.txt" ">>>>>>> Transcript coverage profile"
+
+assert_file_contains "counts.txt" "ENSG00000125841.12"
+
+assert_file_contains "report.html" "<title>Qualimap report: RNA Seq QC</title>"
+assert_file_contains "report.html" "<h3>Input</h3>"
+assert_file_contains "report.html" "<h3>Reads alignment</h3>"
+assert_file_contains "report.html" "<h3>Reads genomic origin</h3>"
+assert_file_contains "report.html" "<h3>Transcript coverage profile</h3>"
+assert_file_contains "report.html" "<h3>Junction analysis</h3>"
+
+
+cd ..
+rm -r run_qualimap_rnaseq_html
+
+mkdir "run_qualimap_rnaseq_pdf"
+cd "run_qualimap_rnaseq_pdf"
+
+echo "> Running qualimap with pdf output report"
+
+"$meta_executable" \
+    --bam $test_dir/a.bam \
+    --gtf $test_dir/annotation.gtf \
+    --report report.pdf \
+    --counts counts.txt \
+    --qc_results output.txt
+
+echo ">> Checking output"
+assert_file_exists "report.pdf"
+assert_file_exists "counts.txt"
+assert_file_exists "output.txt"
+assert_file_doesnt_exist "report.html"
+
+echo ">> Checking if output is empty"
+assert_file_not_empty "report.pdf"
+assert_file_not_empty "counts.txt"
+assert_file_not_empty "output.txt"
+
+cd ..
+rm -r run_qualimap_rnaseq_pdf
+
+mkdir "run_qualimap_rnaseq"
+cd "run_qualimap_rnaseq"
+
+echo "> Running qualimap without report and counts output"
+
+"$meta_executable" \
+    --bam $test_dir/a.bam \
+    --gtf $test_dir/annotation.gtf \
+    --qc_results output.txt
+
+echo ">> Checking output"
+assert_file_doesnt_exist "report.pdf"
+assert_file_doesnt_exist "report.html"
+assert_file_doesnt_exist "counts.txt"
+assert_file_exists "output.txt"
+
+echo ">> Checking if output is empty"
+assert_file_not_empty "output.txt"
+
+cd ..
+rm -r run_qualimap_rnaseq
diff --git a/src/qualimap/qualimap_rnaseq/test_data/a.bam b/src/qualimap/qualimap_rnaseq/test_data/a.bam
diff --git a/src/qualimap/qualimap_rnaseq/test_data/annotation.gtf b/src/qualimap/qualimap_rnaseq/test_data/annotation.gtf
@@ -0,0 +1,10 @@
+chr20	HAVANA	transcript	347024	354868	.	+	.	gene_id "ENSG00000125841.12"; transcript_id "ENST00000382291.7"; gene_type "protein_coding"; gene_name "NRSN2"; transcript_type "protein_coding"; transcript_name "NRSN2-202"; level 2; protein_id "ENSP00000371728.3"; transcript_support_level "2"; tag "basic"; tag "appris_principal_1"; tag "CCDS"; ccdsid "CCDS12996.1"; havana_gene "OTTHUMG00000031628.5"; havana_transcript "OTTHUMT00000077446.1";
+chr20	HAVANA	exon	347024	347142	.	+	.	gene_id "ENSG00000125841.12"; transcript_id "ENST00000382291.7"; gene_type "protein_coding"; gene_name "NRSN2"; transcript_type "protein_coding"; transcript_name "NRSN2-202"; exon_number 1; exon_id "ENSE00001831391.1"; level 2; protein_id "ENSP00000371728.3"; transcript_support_level "2"; tag "basic"; tag "appris_principal_1"; tag "CCDS"; ccdsid "CCDS12996.1"; havana_gene "OTTHUMG00000031628.5"; havana_transcript "OTTHUMT00000077446.1";
+chr20	HAVANA	exon	349249	349363	.	+	.	gene_id "ENSG00000125841.12"; transcript_id "ENST00000382291.7"; gene_type "protein_coding"; gene_name "NRSN2"; transcript_type "protein_coding"; transcript_name "NRSN2-202"; exon_number 2; exon_id "ENSE00001491647.1"; level 2; protein_id "ENSP00000371728.3"; transcript_support_level "2"; tag "basic"; tag "appris_principal_1"; tag "CCDS"; ccdsid "CCDS12996.1"; havana_gene "OTTHUMG00000031628.5"; havana_transcript "OTTHUMT00000077446.1";
+chr20	HAVANA	exon	349638	349832	.	+	.	gene_id "ENSG00000125841.12"; transcript_id "ENST00000382291.7"; gene_type "protein_coding"; gene_name "NRSN2"; transcript_type "protein_coding"; transcript_name "NRSN2-202"; exon_number 3; exon_id "ENSE00003710328.1"; level 2; protein_id "ENSP00000371728.3"; transcript_support_level "2"; tag "basic"; tag "appris_principal_1"; tag "CCDS"; ccdsid "CCDS12996.1"; havana_gene "OTTHUMG00000031628.5"; havana_transcript "OTTHUMT00000077446.1";
+chr20	HAVANA	CDS	349644	349832	.	+	0	gene_id "ENSG00000125841.12"; transcript_id "ENST00000382291.7"; gene_type "protein_coding"; gene_name "NRSN2"; transcript_type "protein_coding"; transcript_name "NRSN2-202"; exon_number 3; exon_id "ENSE00003710328.1"; level 2; protein_id "ENSP00000371728.3"; transcript_support_level "2"; tag "basic"; tag "appris_principal_1"; tag "CCDS"; ccdsid "CCDS12996.1"; havana_gene "OTTHUMG00000031628.5"; havana_transcript "OTTHUMT00000077446.1";
+chr20	HAVANA	start_codon	349644	349646	.	+	0	gene_id "ENSG00000125841.12"; transcript_id "ENST00000382291.7"; gene_type "protein_coding"; gene_name "NRSN2"; transcript_type "protein_coding"; transcript_name "NRSN2-202"; exon_number 3; exon_id "ENSE00003710328.1"; level 2; protein_id "ENSP00000371728.3"; transcript_support_level "2"; tag "basic"; tag "appris_principal_1"; tag "CCDS"; ccdsid "CCDS12996.1"; havana_gene "OTTHUMG00000031628.5"; havana_transcript "OTTHUMT00000077446.1";
+chr20	HAVANA	exon	353210	354868	.	+	.	gene_id "ENSG00000125841.12"; transcript_id "ENST00000382291.7"; gene_type "protein_coding"; gene_name "NRSN2"; transcript_type "protein_coding"; transcript_name "NRSN2-202"; exon_number 4; exon_id "ENSE00001822456.1"; level 2; protein_id "ENSP00000371728.3"; transcript_support_level "2"; tag "basic"; tag "appris_principal_1"; tag "CCDS"; ccdsid "CCDS12996.1"; havana_gene "OTTHUMG00000031628.5"; havana_transcript "OTTHUMT00000077446.1";
+chr20	HAVANA	CDS	353210	353632	.	+	0	gene_id "ENSG00000125841.12"; transcript_id "ENST00000382291.7"; gene_type "protein_coding"; gene_name "NRSN2"; transcript_type "protein_coding"; transcript_name "NRSN2-202"; exon_number 4; exon_id "ENSE00001822456.1"; level 2; protein_id "ENSP00000371728.3"; transcript_support_level "2"; tag "basic"; tag "appris_principal_1"; tag "CCDS"; ccdsid "CCDS12996.1"; havana_gene "OTTHUMG00000031628.5"; havana_transcript "OTTHUMT00000077446.1";
+chr20	HAVANA	stop_codon	353633	353635	.	+	0	gene_id "ENSG00000125841.12"; transcript_id "ENST00000382291.7"; gene_type "protein_coding"; gene_name "NRSN2"; transcript_type "protein_coding"; transcript_name "NRSN2-202"; exon_number 4; exon_id "ENSE00001822456.1"; level 2; protein_id "ENSP00000371728.3"; transcript_support_level "2"; tag "basic"; tag "appris_principal_1"; tag "CCDS"; ccdsid "CCDS12996.1"; havana_gene "OTTHUMG00000031628.5"; havana_transcript "OTTHUMT00000077446.1";
+chr20	HAVANA	UTR	347024	347142	.	+	.	gene_id "ENSG00000125841.12"; transcript_id "ENST00000382291.7"; gene_type "protein_coding"; gene_name "NRSN2"; transcript_type "protein_coding"; transcript_name "NRSN2-202"; exon_number 1; exon_id "ENSE00001831391.1"; level 2; protein_id "ENSP00000371728.3"; transcript_support_level "2"; tag "basic"; tag "appris_principal_1"; tag "CCDS"; ccdsid "CCDS12996.1"; havana_gene "OTTHUMG00000031628.5"; havana_transcript "OTTHUMT00000077446.1";
diff --git a/src/qualimap/qualimap_rnaseq/test_data/script.sh b/src/qualimap/qualimap_rnaseq/test_data/script.sh
@@ -0,0 +1,10 @@
+# qualimap test data
+
+# Test data was obtained from https://github.com/snakemake/snakemake-wrappers/raw/master/bio/qualimap/rnaseq/test
+
+if [ ! -d /tmp/snakemake-wrappers ]; then
+  git clone --depth 1 --single-branch --branch master https://github.com/snakemake/snakemake-wrappers /tmp/snakemake-wrappers
+fi
+
+cp -r /tmp/snakemake-wrappers/bio/qualimap/rnaseq/test/mapped/a.bam src/qualimap/qualimap_rnaseq/test_data
+cp -r /tmp/snakemake-wrappers/bio/qualimap/rnaseq/test/annotation.gtf src/qualimap/qualimap_rnaseq/test_data