How can I call both somatic and germline variants from bam file?
1
2
Entering edit mode
3.8 years ago
svp ▴ 680

I have 'analysis ready bam file created using 'Best Practice Workflows'. My next plan is to get both somatic and germline mutation from bam file. Is there any tool that calls both somatic and germline variants from the bam file.

I have seen Somatic and Germline from broad institute that calls germline and somatic mutation separately.

Does Haplotyecaller outputs both germline and somatic variants ?

exome variant calling somatic germline ngs • 1.7k views
ADD COMMENT
4
Entering edit mode
3.8 years ago
godth13teen ▴ 70

they are called separately because they serve different purposes. The germline variant only need 1 sample to compare with the reference genome, while the somatic variants only found in tumor tissue, they need to be compared with other normal tissue of the same person.

So the answer to your question "Does Haplotyecaller outputs both germline and somatic variants ?" is: No, HaplotypeCaller only work for germline variants, because it compared one genome with the reference genome.

If you think that: "so if the tissue used for sequencing is the tumor one, can HaplotypeCaller detect it?" The answer is still no, because the tissue can also be heterogeneous, the somatic variants occur at quite low frequency, so it's likely to be filtered out by the QC of germline variant. You will need mutect with control and tumor tissue for the somatic variants.

ADD COMMENT

Login before adding your answer.

Traffic: 2711 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