This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to create a data of genotypes for each SNP for each individual

Hello everyone,

I would like to create a data followed by the genotype on each SNP for each individual.

Does someone knows how to get the next followed structure? or Is there any existing tool that provides the below data?

Individuals ID_SPNs Genotype
  ind_1     rs01    AA
  ind_1     rs02    TT
  ind_1     rs03    CC
  :
  .     
  ind_2     rs01    AG
  ind_2     rs02    TC
  ind_2     rs03    AC
  :
  .
  ind_n     rs01    GG
  ind_n     rs02    CC
  ind_n     rs03    AA

I have vcf and plink (bed, fam and bim) format files.

Hopefully someone can give me a hint!

Thanks,

Ale.

geno plink genotype
• 1,559 views • link
updated by Ram Ram • written by Ale Lope • Ale Lope

It looks like the format is 23andme format?

• 0 views • link
written by Shicheng Guo Shicheng Guo

1 answer

If you are fine with the records being grouped by SNP first and individuals second, rather than the other way around,

plink --bfile ... --recode lgen

generates almost exactly what you want (you'd then just need to get rid of the first chromosome-# column, and the space between the two alleles; this can be done with a bash one-liner).

• 0 views • link
written by chrchang523 chrchang523

Thank you so much, it works perfectly!

• 0 views • link
written by Ale Lope • Ale Lope

Log in to answer this question.

More posts like this