diff --git a/variants/germline_vc.sh b/variants/germline_vc.sh index ac2c1d70c0be82583ba7251d36a18c4c78c8c571..c81b05ff6bb04616e863b3262f96f81f305f7ef0 100755 --- a/variants/germline_vc.sh +++ b/variants/germline_vc.sh @@ -86,7 +86,9 @@ then for i in *.bam; do bamlist="$bamlist --bam ${PWD}/${i}" done - cut -f 1 $fbsplit | parallel --memfree 2G --delay 1 --jobs 0 "freebayes -f ${index_path}/genome.fa --min-mapping-quality 0 --min-base-quality 20 --min-coverage 10 --min-alternate-fraction 0.01 -C 3 --use-best-n-alleles 3 -r {} ${bamlist} > fb.{}.vcf" + #--memfree 2G for DNANexus + + cut -f 1 $fbsplit | parallel --delay 1 --jobs 0 "freebayes -f ${index_path}/genome.fa --min-mapping-quality 0 --min-base-quality 20 --min-coverage 10 --min-alternate-fraction 0.01 -C 3 --use-best-n-alleles 3 -r {} ${bamlist} > fb.{}.vcf" vcf-concat fb.*.vcf | vcf-sort | vcf-annotate -n --fill-type | bcftools norm -c s -f ${reffa} -w 10 -O z -o ${pair_id}.fb.vcf.gz - elif [[ $algo == 'platypus' ]] then diff --git a/variants/somatic_vc.sh b/variants/somatic_vc.sh index d3701aab91542be8219e7cce34bb0fe6ca6d91f2..f3461f0fc65441eaaf0c98a0d668c3514397caa0 100755 --- a/variants/somatic_vc.sh +++ b/variants/somatic_vc.sh @@ -142,6 +142,7 @@ then elif [ $algo == 'shimmer' ] then module load R/3.6.1-gccmkl + module rm perl/5.18.2 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