diff --git a/variants/cnvkit.sh b/variants/cnvkit.sh index 84903ed49f14bd172c9897298b38a082af9a143d..b98711c5f9598ff18a12c10bbce3ccd199ad4e2a 100755 --- a/variants/cnvkit.sh +++ b/variants/cnvkit.sh @@ -56,6 +56,7 @@ echo "${targets}antitargets.bed" source /etc/profile.d/modules.sh module load cnvkit/0.9.0 bedtools/2.26.0 +unset DISPLAY cnvkit.py coverage ${sbam} ${targets}targets.bed -o ${pair_id}.targetcoverage.cnn cnvkit.py coverage ${sbam} ${targets}antitargets.bed -o ${pair_id}.antitargetcoverage.cnn cnvkit.py fix ${pair_id}.targetcoverage.cnn ${pair_id}.antitargetcoverage.cnn ${normals} -o ${pair_id}.cnr diff --git a/variants/somatic_vc.sh b/variants/somatic_vc.sh index a45ad2de35a444556e75c6fa158b92264bd02be7..3e43d778333131c712a7922ca334845d97e03928 100755 --- a/variants/somatic_vc.sh +++ b/variants/somatic_vc.sh @@ -99,6 +99,7 @@ if [ $algo == 'virmid' ] virmid -R ${reffa} -D ${tumor} -N ${normal} -s ${cosmic} -t $SLURM_CPUS_ON_NODE -M 2000 -c1 10 -c2 10 perl $baseDir/addgt_virmid.pl ${tumor}.virmid.som.passed.vcf perl $baseDir/addgt_virmid.pl ${tumor}.virmid.loh.passed.vcf + module rm java/oracle/jdk1.7.0_51 module load snpeff/4.3q vcf-concat *gt.vcf | vcf-sort | vcf-annotate -n --fill-type -n | java -jar $SNPEFF_HOME/SnpSift.jar filter '((NDP >= 10) & (DDP >= 10))' | perl -pe "s/TUMOR/${tid}/g" | perl -pe "s/NORMAL/${nid}/g" | bgzip > ${pair_id}.virmid.vcf.gz fi @@ -129,6 +130,7 @@ then sambamba mpileup --tmpdir=./ -t $SLURM_CPUS_ON_NODE ${normal} --samtools "-C 50 -f ${reffa}" > n.mpileup VarScan somatic n.mpileup t.mpileup vscan --output-vcf 1 VarScan copynumber n.mpileup t.mpileup vscancnv + module rm java/oracle/jdk1.7.0_51 module load snpeff/4.3q vcf-concat vscan*.vcf | vcf-sort | vcf-annotate -n --fill-type -n | java -jar $SNPEFF_HOME/SnpSift.jar filter '((exists SOMATIC) & (GEN[*].DP >= 10))' | perl -pe "s/TUMOR/${tid}/" | perl -pe "s/NORMAL/${nid}/g" | bgzip > ${pair_id}.varscan.vcf.gz fi @@ -138,6 +140,7 @@ then module load shimmer/0.1.1 samtools/1.6 vcftools/0.1.14 shimmer.pl --minqual 25 --ref ${reffa} ${normal} ${tumor} --outdir shimmer 2> shimmer.err perl $baseDir/add_readct_shimmer.pl + module rm java/oracle/jdk1.7.0_51 module load snpeff/4.3q vcf-annotate -n --fill-type shimmer/somatic_diffs.readct.vcf | java -jar $SNPEFF_HOME/SnpSift.jar filter '(GEN[*].DP >= 10)' | perl -pe "s/TUMOR/${tid}/" | perl -pe "s/NORMAL/${nid}/g" | bgzip > ${pair_id}.shimmer.vcf.gz fi