This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How To Get The Genotype From Sequence Data?

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!

genotyping nucleotide sequence haplotype
• 5,665 views • link
updated by Biostar Biostar • written by Stingery • Stingery

2 answers

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:

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

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.

• 0 views • link
written by Khader Shameer Khader Shameer

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

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

Log in to answer this question.

More posts like this