This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Genotype Info From 1000 Genomes - Confused By Forward Strand Naming

Hi, I'm trying to calculate the genotype frequencies of a couple of CYP2C19 alleles (rs4244285 + rs12248560) using the 1000 genomes data.

However, I'm confused by the "forward strand" naming convention.

I want to determine the % of people with *2 / *17 but every list I generate has the same person repeated twice it seems.

For eg. I get this on a list for re12248560:

  • HG00099 (F) A|A
  • HG00099 (F) T|C

The same individual it seems but which is it? A|T or A|C or A|A? How can I tell what the "typical" genotype would be?

I want to simply note the diplotype.

Many thanks for your help!

1000genomes strand genotyping
• 2,688 views • link
updated by Biostar Biostar • written by winners • winners

1 answer

Can you describe which steps you do to get the genotypes?

It seems that close to this SNPs, there is a deletion:

# note: you need the latest tabix version from svn for this to work, 
# otherwise you will get an error complaining that the index is too old.
$: tabix <ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20110521/ALL.chr10.phase1_release_v3.20101123.snps_indels_svs.genotypes.vcf.gz> 10:96521657-96521657| cut -f 1,2,3,4,5,12
10  96497183    MERGED_DEL_2_60889  A   <DEL>   0|0:0.000:0,0,0
10  96521657    rs12248560  C   T   1|0:1.000:-3.85,-0.00,-5.00

Maybe your script is confused by the deletion. Try to filter out anything that is not a SNP.

• 0 views • link
written by Giovanni M Dall'Olio Giovanni M Dall'Olio

A description of how you are generating your list (even a code example) would help a lot in trying to figure out what might be going on

• 0 views • link
written by Laura Laura

Log in to answer this question.

More posts like this