Question: BQSR Recalibration for Mycobacterium Genome Data
0
gravatar for arjun.murthy779
6 months ago by
arjun.murthy77910 wrote:

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 • 248 views
ADD COMMENTlink modified 6 months ago by 54354165620 • written 6 months ago by arjun.murthy77910
2
gravatar for 543541656
6 months ago by
54354165620
54354165620 wrote:

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

ADD COMMENTlink written 6 months ago by 54354165620

Thank you! This seems like it should work

ADD REPLYlink written 6 months ago by arjun.murthy77910
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: 1013 users visited in the last hour