How to calculate indels per individual?
1
0
Entering edit mode
2.2 years ago

I have a multi sample VCF file from whole genome sequencing and would like to generate the following table of average variant count per individual: enter image description here

Is there a tool with which you can specify variant type, like deletion, SNP, insertion, Alu, SVA etc and count number of occurance of such variants in each individual?

vcf • 425 views
ADD COMMENT
1
Entering edit mode
2.2 years ago
 $ bcftools stats --samples - input.vcf.gz | grep -w PSC |  cut -f3,9
[3]sample  [9]nIndels
S1         2
S2         0
S3         0
S4         3
S5         0
ADD COMMENT

Login before adding your answer.

Traffic: 1966 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6