This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Best pratices for variant calling on deep sequencing TruSeq Custom Amplicon data

Hi，

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

amplicon sequencing variant calling
• 3,528 views • link
updated by igor igor • written by adp7 • adp7

1 answer

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

• 0 views • link
written by igor igor

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

• 0 views • link
written by fanli.gcb fanli.gcb

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.

• 0 views • link
written by igor igor

Log in to answer this question.

More posts like this