This is a test version of Biostars. For the public version, visit https://www.biostars.org.
BQSR Recalibration for Mycobacterium Genome Data

Hi, We are trying to run an SNP Analysis on Mtb whole genome. Currently, we are stuck at the Base Recalibration Step. It says that we require a DbSNP file to generate the BQSR.table for Applying BQSR. The DbSNP file for Mycobacterium tuberculosis does not exist.

Any advice on how to move forward with the Base Recalibration without the DBSNP would be greatly appreciated. Thanks!

snp bqsr
• 1,634 views • link
updated by 543541656 • 543541656 • written by arjun.murthy779 • arjun.murthy779

1 answer

According to what the gatk forum said,you can bootstrap a database of known SNPs. Here's how it works:

First do an initial round of SNP calling on your original, unrecalibrated data. Then take the SNPs that you have the highest confidence in and use that set as the database of known SNPs by feeding it as a VCF file to the base quality score recalibrator. Finally, do a real round of SNP calling with the recalibrated data. These steps could be repeated several times until convergence.

the site will help you https://gatkforums.broadinstitute.org/gatk/discussion/1706/best-recommendation-for-base-recalibration-on-non-human-data

• 0 views • link
written by 543541656 • 543541656

Log in to answer this question.

More posts like this