Tagged with #samtools
0 documentation articles | 0 announcements | 2 forum discussions


No posts found with the requested search criteria.
No posts found with the requested search criteria.
Comments (2)

Hello All,

I'm a new user of haplotypecaller. I am currently comparing it to unifiedgenotyper and samtools. When I do a count of SNPs for all DPs, it looks like the is a multiple of 3 issues. There seems to be a lot more SNPs with a DP of 3, 6, 9, ... compared to 1, 2, 4, 5, 7, 8 ....

DP #ofSNP
1 43
2 41
3 105
4 24
5 18
6 232
7 38
8 66
9 3749

I was wondering if this was a normal behavior and what caused it.

Yannick

Comments (1)

Hi all,

I am comparing the variant calls from samtools and GATK. For samtools, I have been using quality score cut-offs 100 for SNP and 1000 for INDEL (quite stringent) and as a result, many variants are excluded after filtering. In case of GATK, I have been using our default setting (99% sensitivity for SNPs and 95% sensitivity for INDEL) and included only the variants with FILTER field "PASS". I was wondering if there is any more stringent filters that I can apply and that could be equivalent to samtools QS thresholds since it does not look like this is a fair comparison. Any of your suggestions will be appreciated.