Dear all,
I'm wondering how to do a good SNP and indel calling. There are a lot of aligners each contains different parameters and gives different results.
There is always the dilemma between results that is significant but partial and results that are wide but contain false-positives.
In 1000-genome project site they say about SNP calling methods:
"We found empirically that, given the state of current methods, the consensus of multiple primary call sets from different methods proved of higher quality than any of the primary call sets themselves. ... Overall, genotypes obtained through a consensus procedure are estimated to have 30% fewer errors than those generated by any single caller."
- I wonder if there's a tool that does alignments with different parameters/ tools and then merges the results using statistical methods?
- What would you say are the best SNP calling tools? (bowtie? GATK? samtools? others?...)
Thanks,
Rachelly.
I'm wondering how to do a good SNP and indel calling. There are a lot of aligners each contains different parameters and gives different results.
There is always the dilemma between results that is significant but partial and results that are wide but contain false-positives.
In 1000-genome project site they say about SNP calling methods:
"We found empirically that, given the state of current methods, the consensus of multiple primary call sets from different methods proved of higher quality than any of the primary call sets themselves. ... Overall, genotypes obtained through a consensus procedure are estimated to have 30% fewer errors than those generated by any single caller."
- I wonder if there's a tool that does alignments with different parameters/ tools and then merges the results using statistical methods?
- What would you say are the best SNP calling tools? (bowtie? GATK? samtools? others?...)
Thanks,
Rachelly.
Comment