Question: What Agilent Interval Files (.Bed) Should I Use For Exome Variant Calling With Gatk?
gravatar for newDNASeqer
7.2 years ago by
United States
newDNASeqer710 wrote:

I've just read this post on Biostars, but I'm still unsure what intervals from Agilent I should use for my human exome variant calling with GATK. The exome enrichment kit our lab uses is Agilent Exome Capture Kit Sure Select XT Target Enrichment System for Illumina Paired –End Sequencing Library. Version XTV4.

I was able to find a list of different designs on the Agilent website eArray I am assuming this design named "SureSelect Human All Exon V4" is for the kit we are using, there're a few .bed files in this design such as "regions, alltracks, covered, etc". Should I use the one named "covered" for my variant calling?

ps: my reference genome is human hg19

Sorry for the naive questions, please be gentle. thanks for reply.

agilent bed gatk • 7.2k views
ADD COMMENTlink modified 11 months ago by Karma300 • written 7.2 years ago by newDNASeqer710

Why cant we use the default interval list as mentioned here here

ADD REPLYlink written 3.3 years ago by Karma300

You can use regions mentioned in the website

ADD REPLYlink written 11 months ago by Karma300

I'd use any of them + 100bp to the left and to the right

ADD REPLYlink written 11 months ago by German.M.Demidov1.9k

Why is this, is this because you can enrich pretty effectively for sequences that are 100bp on wither side of the probe because of some kind of overhang effect?

ADD REPLYlink written 11 months ago by curious470

sure, usually regions in some area around the probe is also covered with reads and may be efficiently used for SNVs detection, some people use 200 bp, some people calculate effective enrichment looking at the results instead of provided regions file

ADD REPLYlink written 11 months ago by German.M.Demidov1.9k
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: 1323 users visited in the last hour