Question: Would I want to call variants on hs37d5 for any reason?
0
gravatar for Manuel
3.7 years ago by
Manuel370
Germany
Manuel370 wrote:

Of course, one should use hs37d5.fa for read mapping.

However, should I bother to call variants on the decoy contig hs37d5? GATK Haplotype Caller is spending ages on the contig and I will never use the variant calls.

Would this contig be of any use downstream? Would it make sense to use the GATK variant recalibration on the variant calls of this contig or could they be safely ignored?

Thanks!

variant calling genome • 1.3k views
ADD COMMENTlink modified 3.7 years ago by mjaeger0 • written 3.7 years ago by Manuel370
0
gravatar for mjaeger
3.7 years ago by
mjaeger0
mjaeger0 wrote:

Hi,

we included hs37d5 for the mapping step to get rid of mismapped reads and to increase the SNV-calling sensitivity. See also GATK forum decoy discussion

For downstream analysis we ignore scaffolds and contigs, since most databases and tools are unaware of these.

ADD COMMENTlink written 3.7 years ago by mjaeger0
Please log in to add an answer.

Help
Access

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