Question: Is this a good workflow to find denovo mutations in trios?
gravatar for lait
3.4 years ago by
lait160 wrote:

I have around 100 trios for which WES was done. My goal is to find denovo mutations in the child associated with each trios. So first I will do the following steps:

  1. Alignment to reference genome
  2. marking duplicates
  3. base recalibration
  4. realigning indels
  5. Haplotype caller per sample with the -ERC GVCF option (this will call the ReadBackedPhasing, correct?)
  6. Joint genotyping
  7. Varinat recalibration
  8. Genotype refinement workflow, where pedegree information is used and de novos are annotated using VariantAnnotator.

    1- Do you think this workflow is efficient to find denovos?

    2- Are the variants in the output vcf file produced after step 8 already phased ? (because ReadBackedPhasing was already used in step 5)

    3- Do I need to use PhaseByTransmission afterwards after step 8 ?

Thanks for all your feedback!

Many thanks

ADD COMMENTlink written 3.4 years ago by lait160

Hello lait!

It appears that your post has been cross-posted to another site:

This is typically not recommended as it runs the risk of annoying people in both communities.

ADD REPLYlink written 3.4 years ago by Pierre Lindenbaum134k

Hi Lalit, are there any updates for this question?

ADD REPLYlink written 7 weeks ago by dare_devil1.4k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1784 users visited in the last hour