Question: How To Get The Genotype From Sequence Data?
2
gravatar for Stingery
6.8 years ago by
Stingery20
Stingery20 wrote:

I have a data set of unphased sequence data, i.e. one nucleotide sequence per locus per individual. From there, I would like to compute the genotype, because a wide range of statistical software works with the genotype, but not the sequence. How can I do that? Is there a programme, which I could employ?

Thanks!

ADD COMMENTlink modified 4.0 years ago by Biostar ♦♦ 20 • written 6.8 years ago by Stingery20
3
gravatar for Giovanni M Dall'Olio
6.8 years ago by
London, UK
Giovanni M Dall'Olio26k wrote:

You have to do a type of analysis called SNP calling. There are many tools to do it, depending on the type of data you have. A good starting point may be bowtie used in combination with samtools; have a look at its tutorial, it should be a good primer. Another tool is gatk.

Have also a look at the other discussions in this website. For example:

ADD COMMENTlink modified 6.8 years ago • written 6.8 years ago by Giovanni M Dall'Olio26k
1
gravatar for Khader Shameer
6.8 years ago by
Manhattan, NY
Khader Shameer18k wrote:

If you have a VCF file after QC, Alignment/Mapping and Variant calling steps you can get the genotype data from it (using utility or use vcftools to filter based on various conditions (Read about VCF specifications here and VCFtools manuscript here) . If you have any missing genotypes, you can perform genotype imputation. If you are new to analysis of NGS data refer to Giovanni's answer and other sequencing related discussions at BioStar.

ADD COMMENTlink modified 6.8 years ago • written 6.8 years ago by Khader Shameer18k

+1 for the nice answer and for the new photo plus hair cut! :-)

ADD REPLYlink modified 6.8 years ago • written 6.8 years ago by Giovanni M Dall'Olio26k

Thanks Giovanni :) !

ADD REPLYlink written 6.8 years ago by Khader Shameer18k
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: 676 users visited in the last hour