Question: Joint variant calling
gravatar for arash.iranzadeh1980
19 months ago by
arash.iranzadeh198010 wrote:


I have 1500 bacterial clean and sorted alignment bam file and I want to call variant and produce a multisample vcf file including 1500 samples. I am going to use GATK and bcftools and then compare results. Bacteria are haploid. In the following commands, -ploidy must be set to 1 or 1500?

bcftools call -m -v --ploidy 1 -O z -o variant.vcf.gz mpileup.bcf

java -jar $GATK -T UnifiedGenotyper -R reference.fa -I bam.list -glm BOTH -stand_call_conf 30 -ploidy 1 -o call.vcf

joint variant calling • 1.0k views
ADD COMMENTlink modified 19 months ago • written 19 months ago by arash.iranzadeh198010

why you are not using HaplotypeCaller?

The HaplotypeCaller is a more recent and sophisticated tool than the UnifiedGenotyper. Its ability to call SNPs is equivalent to that of the UnifiedGenotyper, its ability to call indels is far superior, and it is now capable of calling non-diploid samples. As of GATK version

3.3, we recommend using HaplotypeCaller in all cases, with no exceptions.

--sample_ploidy / -ploidy

Ploidy per sample. For pooled data, set to (Number of samples in each pool * Sample Ploidy). Sample ploidy - equivalent to number of chromosome copies per pool. For pooled experiments this should be set to the number of samples in pool multiplied by individual sample ploidy.

ADD REPLYlink modified 18 months ago • written 18 months ago by Medhat8.2k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2183 users visited in the last hour