Genotype calls for a list of variable sites
1
0
Entering edit mode
7.2 years ago
Jautis ▴ 580

Hi,

I have a list of variable sites identified in a previous study of the same population. I would like to call genotypes in a set of new samples specifically at those sites. How can I do so?

I think this is a mix between -L commands to select a particular chromosome region and using a dbSNP file to give prior information about known variants. The difference is that I am only interested in testing the known variants, not in identifying new ones, and I don't care whether the site would pass filtering in the new dataset, only that the genotype call for an individual is reliable assuming the site is known to be variable.

Thank you very much!

genotype-calling GATK variant-calling • 1.3k views
ADD COMMENT
0
Entering edit mode
7.2 years ago
ADD COMMENT
0
Entering edit mode

That's useful, but is there a way to always call the sites, even if there isn't enough information to identify a variable site in the new dataset?

ADD REPLY

Login before adding your answer.

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