.command.sh - Workdir - NFCORE_RNAFUSION:RNAFUSION:QC_WORKFLOW:SEQTOOL_QC (aih-tih-sc-5ceb9a-R1_A23YTGFLT4_1)

File Info

Filename: .command.sh
Full Path: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/5d/036cbfdd1d70687d9e258765fe5ca9/.command.sh
Size: 963 bytes
Task: NFCORE_RNAFUSION:RNAFUSION:QC_WORKFLOW:SEQTOOL_QC (aih-tih-sc-5ceb9a-R1_A23YTGFLT4_1)
Attempt

Content

#!/bin/bash -euo pipefail
set -euxo pipefail

mkdir aih-tih-sc-5ceb9a-R1_A23YTGFLT4_1

all_qc \
-bam "aih-tih-sc-5ceb9a-R1_A23YTGFLT4_1.md.bam" \
-beds Twist_rna_exome_AND_fusion_hg38.bed \
-dbsnp_vcf "Homo_sapiens_assembly38.dbsnp138.vcf.gz" \
-contamination_bed "dbsnp_hg38_contam_sites.bed" \
-output_folder aih-tih-sc-5ceb9a-R1_A23YTGFLT4_1 \
-reference "ref_genome.fa" \
-region 'chr20'

merge_qc -folders aih-tih-sc-5ceb9a-R1_A23YTGFLT4_1/* -output aih-tih-sc-5ceb9a-R1_A23YTGFLT4_1/aih-tih-sc-5ceb9a-R1_A23YTGFLT4_1_full_basic_stats.tsv

head -n1 aih-tih-sc-5ceb9a-R1_A23YTGFLT4_1/aih-tih-sc-5ceb9a-R1_A23YTGFLT4_1_full_basic_stats.tsv > aih-tih-sc-5ceb9a-R1_A23YTGFLT4_1/aih-tih-sc-5ceb9a-R1_A23YTGFLT4_1_basic_stats.tsv
grep -e "$(basename Twist_rna_exome_AND_fusion_hg38.bed)" aih-tih-sc-5ceb9a-R1_A23YTGFLT4_1/aih-tih-sc-5ceb9a-R1_A23YTGFLT4_1_full_basic_stats.tsv >> aih-tih-sc-5ceb9a-R1_A23YTGFLT4_1/aih-tih-sc-5ceb9a-R1_A23YTGFLT4_1_basic_stats.tsv