#!/bin/bash -euo pipefail
expression_qc_and_signatures.py \
    --signatures_master signatures_master.json \
    --tx2gene Homo_sapiens.GRCh38.102.tx2gene_full.tsv \
    --target_panel twist_rna_exome_target_regions_hg38_annotated.bed \
    --sample_id 659_bCu-T1-TRNA-1_B23WHTKLT4_1 \
    --abundance 659_bCu-T1-TRNA-1_B23WHTKLT4_1/abundance.tsv \
    --output_qc 659_bCu-T1-TRNA-1_B23WHTKLT4_1.sample_qc.tsv \
    --output_scores 659_bCu-T1-TRNA-1_B23WHTKLT4_1.sample_scores.tsv \
    --output_gene_table 659_bCu-T1-TRNA-1_B23WHTKLT4_1.gene_table.tsv \


cat <<-END_VERSIONS > versions.yml
"NFCORE_RNAFUSION:RNAFUSION:KALLISTO_TO_TX2GENE_QC":
    python: $(python3 --version | sed 's/Python //g')
    expression_qc_and_signatures: 1.0.0
END_VERSIONS