SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
Exome enrichment, How does Nextera Exome Enrichment compares to TrueseqNimbleAgilent? dottomarco Sample Prep / Library Generation 5 12-19-2013 06:20 AM
Too many Variants called by HaplotypeCaller GATK drmaly Bioinformatics 2 12-05-2013 08:56 PM
Nextera Exome Enrichment kit dummyseq Illumina/Solexa 8 10-03-2013 05:46 AM
truseq exome enrichment niceday General 1 11-04-2011 02:52 AM
TruSeq Custom Enrichment Kits aleferna Illumina/Solexa 3 08-04-2011 01:16 PM

Reply
 
Thread Tools
Old 11-01-2016, 01:09 AM   #1
bioinf newbie
Member
 
Location: India

Join Date: Feb 2012
Posts: 13
Default Variants called using GATK from samples of Nextera and TruSeq exome enrichment kits

I am doing variant calling for exomes generated from Illumina HiSeq 2000, with Nextera Exome Enrichment kit.
Earlier I used the same variant calling pipeline to call variants from sequences based on TruSeq Exome Enrichment kit.
Now with the Nextera Exome Enrichment kit, I find that the number of variants called using GATK is reduced to almost half as compared to TruSeq.
I am not able to reach a conclusion as to why this drastic reduction in number.

TruSeq sample Nextera sample
Mean coverage 40X 47.5X
GATK variants 53917 27596

Here is the GATK command I used for calling variants for both TruSeq and Nextera samples:
java -Xmx2g -Djava.io.tmpdir=temp -jar $gatk_jar -T UnifiedGenotyper -I $in -R $ucsc_hg19_fasta -D:name,VCF $dbsnp146_hg19_vcf -o $o -dcov 1000 -A AlleleBalance -A BaseCounts -A VariantType -baq CALCULATE_AS_NECESSARY -stand_call_conf 30.0 -stand_emit_conf 10.0 -glm BOTH -L $nextera_target_bed(or $truseq_target_bed )
bioinf newbie is offline   Reply With Quote
Old 11-01-2016, 04:00 AM   #2
HESmith
Senior Member
 
Location: Bethesda MD

Join Date: Oct 2009
Posts: 498
Default

A few possible explanations:

1) higher error rate in the sequence data from the TruSeq sample, producing false-positive variant calls. You can use FastQC to assess the base quality metrics of the two data sets.

2) capture of highly polymorphic loci in the TruSeq but not Nextera enrichment. Compare the covered intervals with Bedtools, and see if there's a large difference in variant frequency between the shared vs unique regions.
HESmith is offline   Reply With Quote
Reply

Tags
gatk, nextera exomes, variant calling

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 04:28 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO