To celebrate the release of GATK 4.0, we are giving away free credits for running the GATK4 Best Practices pipelines in FireCloud, our secure online analysis portal. It’s first come first serve, so sign up now to claim your free credits worth $250. Sponsored by Google Cloud. Learn more at

variant analysis -calling

I am working with plants, and so cannot follow most of your human/cancer specific workflows etc.
Also my question is probably not GATK specifc, but I would be very grateful for any kind of answer. I have used BWA to align a set of reads and the contigs assembled from these reads (de novo). I have used samtools mpileup and GATK unifiedGenotyper on the aligned reads. This results for both in a list of ONLY SNPs, the same ones (aside of sensitivity).
I have also used the samtools steps on the aligned contigs, resulting in a list of ONLY indels.
I can clearly see both SNPs and indels in IGV, both looking at the reads or the contigs.
generally I see that reads are being used for variant calling, so is it "wrong" to use contigs? and why do I not get ANY indels in either tool's output when using reads? any ideas?

Best Answer


Sign In or Register to comment.