How to get all variants in original BAM file with HaplotypeCaller? (No correction and BAMout)
0
0
Entering edit mode
6.4 years ago

Dear All,

As you know when we use HaplotypeCaller in the variant calling pipeline, the script makes corrections and focuses on the true positive variants. and eliminates many vairant so called wrong one or caused by the artifact sequencing. But no matter what I would like to get all the variants comes directly from the original BAM file. How should I use HaplotypeCaller so that my original BAM and the final BAMout have exactly the same variants?

Thanks in advance.

variant calling haplotypecaller sequencing bamout • 1.8k views
ADD COMMENT
1
Entering edit mode

If you want that then why not use samtools mpileup or something like that instead? If you don't want to output likely variants but all differences then don't use a variant caller.

ADD REPLY

Login before adding your answer.

Traffic: 1881 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6