From ff59f96d8ead6a9828b259544c39180ce063ec0b Mon Sep 17 00:00:00 2001
From: Brandi Cantarel <brandi.cantarel@utsouthwestern.edu>
Date: Tue, 5 Dec 2017 10:55:59 -0600
Subject: [PATCH] adding updates to sv and somatic coding
---
variants/somatic_vc.sh | 7 ++++---
variants/svcalling.sh | 2 --
2 files changed, 4 insertions(+), 5 deletions(-)
diff --git a/variants/somatic_vc.sh b/variants/somatic_vc.sh
index 929dea2..2c1dccb 100644
--- a/variants/somatic_vc.sh
+++ b/variants/somatic_vc.sh
@@ -105,7 +105,7 @@ fi
if [ $algo == 'mutect2' ]
then
module load parallel gatk/3.7 snpeff/4.3q vcftools/0.1.14
- cut -f 1 ${index_path}/genomefile.5M.txt | parallel --delay 2 -j 10 "java -Xmx20g -jar \$GATK_JAR -R ${reffa} -D ${dbsnp} -T MuTect2 -stand_call_conf 30 -stand_emit_conf 10.0 -A FisherStrand -A QualByDepth -A VariantType -A DepthPerAlleleBySample -A HaplotypeScore -A AlleleBalance -I:tumor ${tumor} -I:normal ${normal} --cosmic ${cosmic} -o ${tid}.{}.mutect.vcf -L {}"
+ cut -f 1 ${index_path}/genomefile.5M.txt | parallel --delay 2 -j 10 "java -Xmx20g -jar \$GATK_JAR -R ${reffa} -D ${dbsnp} -T MuTect2 -stand_call_conf 10 -A FisherStrand -A QualByDepth -A VariantType -A DepthPerAlleleBySample -A HaplotypeScore -A AlleleBalance -I:tumor ${tumor} -I:normal ${normal} --cosmic ${cosmic} -o ${tid}.{}.mutect.vcf -L {}"
vcf-concat ${tid}*.vcf | vcf-sort | vcf-annotate -n --fill-type | java -jar \$SNPEFF_HOME/SnpSift.jar filter -p '((FS <= 60) & GEN[*].DP >= 10)' | perl -pe 's/TUMOR/${tid}/' | perl -pe 's/NORMAL/${nid}/g' |bgzip > ${pair_id}.pmutect.vcf.gz
fi
@@ -130,5 +130,6 @@ fi
if [ $algo == 'lancet' ]
then
module load snpeff/4.3q lancet vcftools/0.1.14
- lancet --tumor ${tumor} --normal ${normal} --ref $reffa --bed $target_panel --num-threads 16 > out.vcf
- vcf-concat out.vcf | vcf-sort | vcf-annotate -n --fill-type -n | perl -pe 's/TUMOR/${tid}/' | perl -pe 's/NORMAL/${nid}/g' |bedtools intersect -header -a stdin -b $target_panel |bgzip > ${tid}_${nid}.lancet.vcf.gz
+ lancet --tumor ${tumor} --normal ${normal} --ref $reffa --bed $target_panel --num-threads 16 > lancet.vcf
+ vcf-sort lancet.vcf | vcf-annotate -n --fill-type -n | perl -pe 's/TUMOR/${tid}/' | perl -pe 's/NORMAL/${nid}/g' |bgzip > ${pair_id}.lancet.vcf.gz
+fi
diff --git a/variants/svcalling.sh b/variants/svcalling.sh
index 78029b1..504ff88 100644
--- a/variants/svcalling.sh
+++ b/variants/svcalling.sh
@@ -135,5 +135,3 @@ java -Xmx10g -jar $SNPEFF_HOME/snpEff.jar -no-intergenic -lof -c $SNPEFF_HOME/sn
perl $baseDir/svannot.pl -i ${pair_id}.sv.annot.vcf
bgzip ${pair_id}.sv.annot.vcf
tabix ${pair_id}.sv.annot.vcf.gz
-
-fi
--
GitLab