File Info

Filename
.command.sh
Full Path
s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/8e/4b00ec874028e170f2610139bcc7e8/.command.sh
Size
1.6 KB
Task
DAQ:BAM_POSTDEDUP_QC:SEQTOOL (659_Lx-T1-TDNA-1_A23MHWWLT4_1)
Attempt

Content

View in new tab 
#!/bin/bash -Ceuo pipefail
# Matplotlib cache dir must be writable; default location may not be
export MPLCONFIGDIR=/tmp/matplotlib

all_qc \
    -bam 659_Lx-T1-TDNA-1_A23MHWWLT4_1.md.bam \
    -reference Homo_sapiens_assembly38.fasta \
    -beds xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed xgen-exome-hyb-panel-v2_AND_altera_v3_probes_short_hg38_AND_cdx_spike_custom30off4x_hg38_AND_foresight_clarity_footprint.bed altera_v3_targets_postQC_hg38.bed altera_v3_targets_coding_postQC_hg38.bed altera_v3_targets_noncoding_postQC_hg38.bed xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed FINAL_hotspot_list_07302025_hg38_sorted.bed foresight_clarity_targets_hg38.bed \
    -dbsnp_bed dbsnp138_hg38_common.bed \
    -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz \
    -contamination_bed dbsnp_hg38_contam_sites.bed \
    -output_folder 659_Lx-T1-TDNA-1_A23MHWWLT4_1 \
    -threads 31 \


merge_qc -folders 659_Lx-T1-TDNA-1_A23MHWWLT4_1/* -output 659_Lx-T1-TDNA-1_A23MHWWLT4_1/659_Lx-T1-TDNA-1_A23MHWWLT4_1_full_basic_stats.tsv

head -n1 659_Lx-T1-TDNA-1_A23MHWWLT4_1/659_Lx-T1-TDNA-1_A23MHWWLT4_1_full_basic_stats.tsv > 659_Lx-T1-TDNA-1_A23MHWWLT4_1/659_Lx-T1-TDNA-1_A23MHWWLT4_1_basic_stats.tsv
grep -e "$(basename xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed)" 659_Lx-T1-TDNA-1_A23MHWWLT4_1/659_Lx-T1-TDNA-1_A23MHWWLT4_1_full_basic_stats.tsv >> 659_Lx-T1-TDNA-1_A23MHWWLT4_1/659_Lx-T1-TDNA-1_A23MHWWLT4_1_basic_stats.tsv

cat <<-END_VERSIONS > versions.yml
"DAQ:BAM_POSTDEDUP_QC:SEQTOOL":
    seqtool: 2.0
END_VERSIONS