Hi all,
I'm going through some exome-seq data (10 samples), and following a "standard" pipeline - BWA, HaplotypeCaller, GenotypeGVCF, etc.. nothing too fancy.
My question is on the VQSR step. From my understanding, this works best when you have lots of samples to look at, but since I only have 10 samples, would it be better to simply do some hard-filtering on the dataset? Alternatively, I was thinking of bringing in some additional samples (such as a bunch of 1K genome exomes) and add these into the mix. This would increase the number of samples. Does anyone have any thoughts on this? Also, how many samples would be truly need to do the variant filtering with the VQSR rather than hard-filtering? (I'm having a little bit of a hard time finding what the minimum number of samples for it to work reliably).
Thanks for any advice
I'm going through some exome-seq data (10 samples), and following a "standard" pipeline - BWA, HaplotypeCaller, GenotypeGVCF, etc.. nothing too fancy.
My question is on the VQSR step. From my understanding, this works best when you have lots of samples to look at, but since I only have 10 samples, would it be better to simply do some hard-filtering on the dataset? Alternatively, I was thinking of bringing in some additional samples (such as a bunch of 1K genome exomes) and add these into the mix. This would increase the number of samples. Does anyone have any thoughts on this? Also, how many samples would be truly need to do the variant filtering with the VQSR rather than hard-filtering? (I'm having a little bit of a hard time finding what the minimum number of samples for it to work reliably).
Thanks for any advice
Comment