Skip to content
Snippets Groups Projects
Commit e5a54f52 authored by Brandi Cantarel's avatar Brandi Cantarel
Browse files

update small docker/dnanexus incompatobility

parent ba8a20f1
No related merge requests found
......@@ -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
......
......@@ -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
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment