gatk --java-options '-Xmx40g' HaplotypeCaller -R $REF    --max-alternate-alleles 4  --sample-ploidy 2 -O  all.raw.vcf.gz    -I 1.sorted.dedup.bam  -I 2.sorted.dedup.bam  --tmp-dir tmp