Creating a VCF file with A and B allele
0
0
Entering edit mode
2.9 years ago
PeterKW ▴ 90

Hi,

I have an lgen file that looks like this

FID    ID     SNP    allele1_AB     allele2_AB
Y     131   ARS-5       A                  A
Y     131   ARS-2       A                  B
Y     131   ARS-7       B                  B

From this I want to create a VCF file with the genotypes as 0/0, 0/1 and 1/1.

I am using plink --cow --file wagyu --recode vcf but I get a warning

Warning: At least one VCF allele code violates the official specification;
other tools may not accept the file.  (Valid codes must either start with a
'<', only contain characters in {A,C,G,T,N,a,c,g,t,n}, be an isolated '*', or
represent a breakend.)

How can I go about this without losing data. I can't use --snps-only just-acgt because I have alleles A and B.

Thanks.

VCF plink allele • 791 views
ADD COMMENT

Login before adding your answer.

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