-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathalter_rg_804-805-806-807-808.sh
38 lines (29 loc) · 3.99 KB
/
alter_rg_804-805-806-807-808.sh
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
#!/bin/bash
#PBS -N alter_rg804_805_806_807_808
#PBS -q largeq
#PBS -e alter_rg804_805_806_807_808_script_err
#PBS -o alter_rg804_805_806_807_808_script_log
#PBS -l mem=100gb,walltime=100:00:00,nodes=1:ppn=8
#HSCHED -s hschedd
#############################################################################################################################
cd /leofs/mishl_group/liangl/rare_dis
ref=/leofs/mishl_group/liangl/rare_dis/hg19_gatkorder
bwa=/software/biosoft/software/bwa-0.7.10/bwa
java=/software/biosoft/software/java1.7/jdk1.7.0_45/bin/java
gatk=/software/biosoft/software/gatk/GenomeAnalysisTK.jar
picard_dir=/software/biosoft/software/picard-tools-1.86
#samtools view -h WGC034804U_dedup_realign_recal.bam | sed 's/:WGC/:WGC804/g' | samtools view -bS - > WGC804.bam
#samtools index WGC804.bam
#samtools view -h WGC034805U_dedup_realign_recal.bam | sed 's/:WGC/:WGC805/g' | samtools view -bS - > WGC805.bam
#samtools index WGC805.bam
#samtools view -h WGC034806U_dedup_realign_recal.bam | sed 's/:WGC/:WGC806/g' | samtools view -bS - > WGC806.bam
#samtools index WGC806.bam
#samtools view -h WGC034807U_dedup_realign_recal.bam | sed 's/:WGC/:WGC807/g' | samtools view -bS - > WGC807.bam
#samtools index WGC807.bam
#samtools view -h WGC034808U_dedup_realign_recal.bam | sed 's/:WGC/:WGC808/g' | samtools view -bS - > WGC808.bam
#samtools index WGC808.bam
${java} -Xmx75g -jar ${gatk} -T UnifiedGenotyper -R ${ref}.fa -glm BOTH -I WGC804.bam -I WGC805.bam -I WGC806.bam -I WGC807.bam -I WGC808.bam --dbsnp: ./hg19/dbsnp_138.hg19.vcf -stand_call_conf 30 -stand_emit_conf 10 -o RG804_805_806_807_808Fam_dedup_realign_recal_rawU_allvar.vcf -minIndelCnt 4 -minIndelFrac 0.05 -dcov 1000 2>>RG804_805_806_807_808_err
${java} -Xmx75g -jar ${gatk} -T VariantRecalibrator -R ${ref}.fa -input RG804_805_806_807_808Fam_dedup_realign_recal_rawU_allvar.vcf -mode SNP -resource:hapmap,known=false,training=true,truth=true,prior=15.0 ./hg19/hapmap_3.3.hg19.vcf -resource:omni,known=false,training=true,truth=false,prior=12.0 ./hg19/1000G_omni2.5.hg19.vcf -resource:1000G,known=false,training=true,truth=false,prior=10.0 ./hg19/1000G_phase1.indels.hg19.vcf -resource:dbsnp,known=true,training=false,truth=false,prior=2.0 ./hg19/dbsnp_138.hg19.vcf -an DP -an FS -an MQRankSum -an ReadPosRankSum -tranche 100.0 -tranche 99.9 -tranche 99.0 -tranche 90.0 --maxGaussians 4 -recalFile RG804_805_806_807_808Fam_dedup_realign_recal_rawU_snp.recal -tranchesFile RG804_805_806_807_808Fam_dedup_realign_recal_rawU_snp.tranches -rscriptFile RG804_805_806_807_808Fam_dedup_realign_recal_rawU_snp_plots.R 2>>RG804_805_806_807_808_err
${java} -Xmx75g -jar ${gatk} -T VariantRecalibrator -R ${ref}.fa -input RG804_805_806_807_808Fam_dedup_realign_recal_rawU_allvar.vcf -mode INDEL -resource:mills,known=true,training=true,truth=true,prior=12.0 ./hg19/Mills_and_1000G_gold_standard.indels.hg19.vcf -an DP -an FS -an MQRankSum -an ReadPosRankSum -tranche 100.0 -tranche 99.9 -tranche 99.0 -tranche 90.0 --maxGaussians 4 -recalFile RG804_805_806_807_808Fam_dedup_realign_recal_rawU_indel.recal -tranchesFile RG804_805_806_807_808Fam_dedup_realign_recal_rawU_indel.tranches -rscriptFile RG804_805_806_807_808Fam_dedup_realign_recal_rawU_indel_plots.R 2>>RG804_805_806_807_808_err
${java} -Xmx75g -jar ${gatk} -T ApplyRecalibration -R ${ref}.fa -input RG804_805_806_807_808Fam_dedup_realign_recal_rawU_allvar.vcf -mode SNP --ts_filter_level 99.0 -recalFile RG804_805_806_807_808Fam_dedup_realign_recal_rawU_snp.recal -tranchesFile RG804_805_806_807_808Fam_dedup_realign_recal_rawU_snp.tranches -o RG804_805_806_807_808Fam_dedup_realign_recal_rawU_snpF.vcf 2>>RG804_805_806_807_808_err
${java} -Xmx75g -jar ${gatk} -T ApplyRecalibration -R ${ref}.fa -input RG804_805_806_807_808Fam_dedup_realign_recal_rawU_snpF.vcf -mode INDEL --ts_filter_level 99.0 -recalFile RG804_805_806_807_808Fam_dedup_realign_recal_rawU_indel.recal -tranchesFile RG804_805_806_807_808Fam_dedup_realign_recal_rawU_indel.tranches -o RG804_805_806_807_808Fam_dedup_realign_recal_rawU_snpF_indelF.vcf 2>>RG804_805_806_807_808_err