Biostar Test Site

This is site is used for testing only. Visit: https://www.biostars.org to ask a question.

Vg Call, not calling SNP in all samples
0
3
Entering edit mode
9 months ago

Hi,

I have been trying to use Vg call to call variants from pooled samples. I am interested in using it for indels but I first compared the vg call frequency for two biaalelic SNPs whose frequency was also estimated using a different method (PoolSNP). The estimates from PoolSNP seem quite reliable but only for biallelic SNPs. Vg call produces similar AF for the variants it has called, see plots below. However, it does not call the variant in most samples, again in plot. The minimum sample depth in the PoolSNP calls for both sites is >40. My feeling it has something to do with the fact that Vg call work on samples individually without joint calling. I also attached my approach below, I am using vg-v1.26.1.

vg map -x x.xg -g x.gcsa -f 1.fq.gz  -f  2.fq.gz -t 32 -Z 100 > mapped.gam

vg gamsort mapped.gam -t 32  -i mapped_sorted.gam.gai >mapped_sorted.gam

vg chunk -x x.xg -c 10 -p 3:3716672-3718990 -g -a mapped_sorted.gam -O pg >chunk.pg 2>chunk_err

vg augment chunk.pg chunk_0_3_3716352_3719336.gam -s -A chunk_1_aug.gam > chunk_1_aug.pg

vg snarls chunk_1_aug.pg > chunk_1_aug.snarls

vg pack -x chunk_1_aug.pg -g chunk_1_aug.gam -o chunk_1_aug.pack

vg call chunk_1_aug.pg -a -r chunk_1_aug.snarls -k chunk_1_aug.pack -s sample > sample_calls.vcf


I only show for two SNPs but this happens for many others as well and it never the same samples that are called by Vg call. Please let me know if you have any insights into what is going on. Thanks so much for your time, Arun.

Vg call vg tools vg • 406 views