Question: How to get all variants in original BAM file with HaplotypeCaller? (No correction and BAMout)
gravatar for mhmtgenc85
20 months ago by
mhmtgenc8540 wrote:

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.

ADD COMMENTlink written 20 months ago by mhmtgenc8540

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 REPLYlink written 20 months ago by Devon Ryan94k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1100 users visited in the last hour