How to interpret snp frequency
Hi, I am running VCFtools --freq on a vcf files that contains 17 individuals the output is :
CHROM POS N_ALLELES N_CHR {FREQ}
1 13504 2 2 0.5 0.5
1 19776 2 2 0.5 0.5
1 19788 2 2 0.5 0.5
1 21580 2 6 0 1
1 22158 2 3 0 1
1 23975 2 13 0 1
1 30860 2 1 0 1
1 35510 2 2 0 1
1 38232 2 30 0 1
It doesn't make any sense to me, I manually checked the first variants in the vcf file I get:
13504 present in 1 individual out of 17
19776 present in 1 individual out of 17
19978 present in 1 individual out of 17
21580 present in 3 individual out of 17
22158 present in 3 individual out of 17
1 answer
Log in to answer this question.
More posts like this
-
allele frequency from vcf
written by niyomiw88 •
-
VCF tools --freq output file
written by Inquisitive8995
-
Calculate average parental genotype
written by coleman_jonathan
-
VCF allele frequency is not sensible
written by wrab425 •
-
Calculate the genotype frequency
written by SOHAIL
-
The meaning of '0' allele frequency in vcftools output
written by Joel Wallenius
-
How to solve this error IndexError: index 2 is out of bounds for axis 1 with size 0
written by python •
-
VCF format query
written by AW
-
Piloting population frequencies of different populations
written by kv.shamsudheen •
-
Meaning Of N_Alleles And N_Chr In .Frq Files
written by abhinkraj7 •
You can try with the varscan, for the readcount. Please follow the link here
/tool/Varscan/VarScan.v2.2.11.jar readcounts <mpileup> --min-coverage 20 --min-base-qual 30 --output-file <output file name>OP is looking at variant frequencies in a vcf, are you sure your question about getting read counts is a relevant answer?
The whole idea is to check for the variants at a particular locus. When getting the read counts for the position, what can be done is you can map the read counts with VCF file and check the exact reads supporting your variant calls.
I am not familiar with varScan but I have only access to the variant vcf files not the mpileup