number of SNPs per sample in multi-sample VCF file
0
0
Entering edit mode
12 days ago
reza ▴ 300

I run following command on my multi-sample vcf file including 10 samples

bcftools stats --samples "-" in.vcf.gz

in a part of output, there are a table with following column titles

PSC [2]id   [3]sample   [4]nRefHom  [5]nNonRefHom   [6]nHets    [7]nTransitions [8]nTransversions   [9]nIndels  [10]average depth   [11]nSingletons [12]nHapRef [13]nHapAlt [14]nMissing

sum of which columns is SNP count for each sample.

Thanks

VCF bcftools • 243 views
ADD COMMENT
0
Entering edit mode

it should be nNonRefHom isn't it ?

ADD REPLY
0
Entering edit mode

What's the ploidy of your samples? I would argue SNP count should be nNonRefHom and nHets since heterozygote SNPs are still informative in most cases so long as your samples are at least diploid.

ADD REPLY

Login before adding your answer.

Traffic: 2840 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