Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Multi-sample vs Single sample SNP calling for Linkage analysis meher Bioinformatics 0 10-23-2013 05:13 AM
Multi sample SNP calling meher Bioinformatics 3 11-15-2012 12:28 PM
Multi-sample VCF generated in GATK is not working with novoutil IUPAC sdarko Bioinformatics 0 07-19-2012 04:55 AM
Large-Scale, Multi-Sample SNP Analysis Video DNASTAR Vendor Forum 0 10-26-2011 07:11 AM

Thread Tools
Old 06-04-2018, 12:44 AM   #1
Dr. Deepak
Junior Member
Location: Delhi

Join Date: Aug 2015
Posts: 3
Lightbulb GATK SNP multi-sample without known SNP database ?

I am working on identifying the SNPs in the genome of an insect which is sequenced in three samples. Following the GATK best practice pipeline i have completed the following task.

Quality filtered reads.
Read Alignment for each sample eg. S1, S2, S3.
Sort and mark duplicates in each sample S1, S2, S3.
Realign target and Indel realinger on each sample S1, S2, S3. Now i am stuck with next steps. As it require the base recalibration (BSQR) step. which require the known database of snps. however i dont have any known snp database so i am following the alternative steps mentioned the GATK for non-model organism. That is
I call the raw SNPs in gvcf model from each sample and then make joint calling for a final combined vcf file.
Next, I applied the hard filters on this file and extract good quality SNPs and indels as a reference database of BSQR step.

Now the question is i have three realigned files S1_realigned.bam, S2_realigned.bam, S3_realigned.bam from step 4 and a reference database of SNP and Indels (if i am right) from the step 6. So how i would proceed further. Did i use against each sample separately or make a combine re-calibration table by providing all three in the same command.
Thanks Dr. Deepak
Dr. Deepak is offline   Reply With Quote

base recaibration, gatk, snp call

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 03:59 PM.

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