/share/work/biosoft/bcftools/latest/bin/bcftools mpileup -a 'AD,DP' -Ou -f /share/work/database/ref/Triticum_aestivum/IWGSC_v2.1/iwgsc_refseqv2.1_gene_annotation_200916/iwgsc_refseqv2.1_assembly.fa /share/nas1/renzx/project/zx-20210524-20_BSR_2xiaomai/ref/ann/9.SNP_Analysis/2.bam_process/CH.bam.split.bam /share/nas1/renzx/project/zx-20210524-20_BSR_2xiaomai/ref/ann/9.SNP_Analysis/2.bam_process/CL.bam.split.bam /share/nas1/renzx/project/zx-20210524-20_BSR_2xiaomai/ref/ann/9.SNP_Analysis/2.bam_process/PH.bam.split.bam /share/nas1/renzx/project/zx-20210524-20_BSR_2xiaomai/ref/ann/9.SNP_Analysis/2.bam_process/PL.bam.split.bam|/share/work/biosoft/bcftools/latest/bin/bcftools call -mv -Oz -o ~/calls.vcf.gz
