Question: Best pratices for variant calling on deep sequencing TruSeq Custom Amplicon data
gravatar for adp7
2.0 years ago by
adp710 wrote:


I'm trying to analyze a panel of genes in tumor samples using truseq custom amplicon (TSCA) with average depth ~2000X. I was wondering if classic tools like GATK HC, UG, MUTECT are adapted to this kind of data? Is there any extra filtering step to perform because of the PCR amplification? And if is there variant calling settings more suited for deep sequencing approach?

Thanks in advance

ADD COMMENTlink modified 2.0 years ago by igor7.6k • written 2.0 years ago by adp710
gravatar for igor
2.0 years ago by
United States
igor7.6k wrote:

Classic tools should be fine. However, you have to skip duplicate removal, since these are amplicons, so most reads are expected to be duplicates.

ADD COMMENTlink written 2.0 years ago by igor7.6k

Would you consider this to also be the case for exome sequencing?

ADD REPLYlink written 2.0 years ago by fanli.gcb660

Depends on the library type. Most exomes are capture-based, so the fragments should be randomly selected. There are some exceptions. For example, Agilent SureSelect is amplicon-based, so it's more like TSCA.

GATK Best Practices has the same section for WGS and WES, since the two are usually so similar.

ADD REPLYlink modified 2.0 years ago • written 2.0 years ago by igor7.6k
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: 1156 users visited in the last hour