SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Hapmap wrongly called APOE SNP rs429358??? ymc Bioinformatics 0 09-27-2011 02:14 AM
Pileup shows variants - Not called as SNPs BertieWooster Bioinformatics 10 05-12-2011 11:57 AM
Softwear to visualise snps called from NGS mrxcm3 Bioinformatics 3 01-14-2011 12:25 PM
How to count indels called by SAMtools yuanzhi Bioinformatics 3 10-31-2010 04:05 PM
Determining the number of bases and percent coverage in an aligned sequence kz26 Bioinformatics 0 06-28-2010 09:22 PM

Reply
 
Thread Tools
Old 10-06-2011, 01:04 AM   #1
claratsm
Junior Member
 
Location: Australia

Join Date: Sep 2011
Posts: 1
Default GATK : Drop in number of confidently called bases

We called the SNPs using GATK v1.0.5506 before and decided to recall using the recent version v1.2-4; however, we found a dramatical decrease in number of confidently called bases from ~90 to ~20%.

Can anyone suggest what's the problem?

Here is the old command and summary (taking a small region as an example),

--------------------------------------------------------------------------------------
java -Xmx4g -jar GenomeAnalysisTK-1.2-4-gd9ea764/GenomeAnalysisTK.jar \
-I recalibrated.bam \
-T UnifiedGenotyper \
-l INFO \
-R human_g1k_v37.fa \
--dbsnp dbsnp_132.b37.vcf \
-o test.raw.SNP.vcf \
-metrics test.raw.SNP.variants \
-glm SNP \
-stand_call_conf 50 \
-stand_emit_conf 10 \
-dcov 500 \
-nt 1 \
-L 11:76260000-76270000 \
-S LENIENT

================================
Visited bases 10001
Callable bases 9992
Confidently called bases 152
% callable bases of all loci 99.910
% confidently called bases of all loci 1.520
% confidently called bases of callable loci 1.521
Actual calls made 5
--------------------------------------------------------------------------------------
And the new command
java -Xmx4g -jar GenomeAnalysisTK-1.0.5506/GenomeAnalysisTK.jar \
-I recalibrated.bam \
-T UnifiedGenotyper \
-l INFO \
-R human_g1k_v37.fa
-B:dbsnp,vcf dbsnp_132.b37.vcf \
-o test.oldUG.raw.SNP.vcf \
-metrics test.oldUG.variants \
-glm SNP \
-stand_call_conf 50 \
-stand_emit_conf 10 \
-dcov 500 \
-nt 1 \
-L 11:76260000-76270000 \
-S LENIENT

==================================
Visited bases 10001
Callable bases 9990
Confidently called bases 9988
% callable bases of all loci 99.890
% confidently called bases of all loci 99.870
% confidently called bases of callable loci 99.980
Actual calls made 5
claratsm is offline   Reply With Quote
Reply

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 01:29 AM.


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