File Info

Filename
.command.sh
Full Path
s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/f1/3549e6f0a3cd6883239fe71c82922c/.command.sh
Size
1.1 KB
Task
DAQ:SENTIEON_BQSR:QUALCAL_APPLY (1136_6TL-T1-TDNA-1_A23MHWWLT4_1)
Attempt

Content

View in new tab 
#!/bin/bash -Ceuo pipefail
sentieon driver \
    -t 31 \
    -r Homo_sapiens_assembly38.fasta \
    -i 1136_6TL-T1-TDNA-1_A23MHWWLT4_1.md.bam \
    --read_filter QualCalFilter,table=1136_6TL-T1-TDNA-1_A23MHWWLT4_1.table,indel=false \
    --algo QualCal \
     \
    -k dbsnp_146.hg38.vcf.gz -k Mills_and_1000G_gold_standard.indels.hg38.vcf.gz -k Homo_sapiens_assembly38.known_indels.vcf.gz \
    1136_6TL-T1-TDNA-1_A23MHWWLT4_1.table.post \
    --algo ReadWriter 1136_6TL-T1-TDNA-1_A23MHWWLT4_1.recalibrated.bam

# Sentieon ReadWriter doesn't create index - use samtools
samtools index 1136_6TL-T1-TDNA-1_A23MHWWLT4_1.recalibrated.bam

sentieon driver \
    -t 31 \
    --algo QualCal \
    --plot \
    --before 1136_6TL-T1-TDNA-1_A23MHWWLT4_1.table \
    --after 1136_6TL-T1-TDNA-1_A23MHWWLT4_1.table.post \
     \
    1136_6TL-T1-TDNA-1_A23MHWWLT4_1.csv

sentieon plot QualCal -o 1136_6TL-T1-TDNA-1_A23MHWWLT4_1.pdf 1136_6TL-T1-TDNA-1_A23MHWWLT4_1.csv

cat <<-END_VERSIONS > versions.yml
"DAQ:SENTIEON_BQSR:QUALCAL_APPLY":
    sentieon: $(echo $(sentieon driver --version 2>&1) | sed -e "s/sentieon-genomics-//g")
END_VERSIONS