Question: Extracting heterozygous genotype from vcf file and compare 2 files
0
gravatar for priya
16 months ago by
priya0
priya0 wrote:

Hi all,

I have a vcf file and would like to extract certain chromosomes, say chr16 and only GT = 1/1. I tried using grep and it worked fine. However when I try to use these vcf files created by using grep, to compare the SNPs between 2 such files using vcftools, it shows me an error in the file format, which makes me think grep is not being able to keep the same format of vcf. Could anybody suggest the best possible way to do that?

My end goal is to compare 2 vcf files (from different samples but same chromosome 16) and find an intersection of SNPs between the files containing 1/1 genotype from say chr16.

Best, Priya

genotype snps vcf • 531 views
ADD COMMENTlink modified 16 months ago by Dan Gaston6.8k • written 16 months ago by priya0
0
gravatar for Dan Gaston
16 months ago by
Dan Gaston6.8k
Canada
Dan Gaston6.8k wrote:

You should use one of the many different tool sets out there that have been developed for working with VCF files and doing filtering on them. A variety to check out would be:

VCFTools BCFTools SNPSift

ADD COMMENTlink written 16 months ago by Dan Gaston6.8k
Please log in to add an answer.

Help
Access

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