can VarScan2 do joint somatic calling?
1
0
Entering edit mode
3.4 years ago
Zoey Zou ▴ 40

Hi,

In my project, I am going to use the SNVs of different tumor regions to build the phylogenetic tree for that tumor. I am doing joint calling in different software separately. Now, I want to try VarScan2. I have one normal bam file and 10 tumor bam files from that patient. But I find that the import format for VarScan2 is:

samtools mpileup -f reference -q 1 -B normal.bam tumor.bam >normal-tumor.mpileup 
java -jar VarScan.jar somatic normal-tumor.mpileup normal-tumor.varScan.output --mpileup 1

Does it mean that I can only do somatic calling for one pair of normal-tumor sequence each time? How can I get the joint calling result for that?

I would appreciate any help.

Regards,

Zoey

SNP VarScan snp somatic calling • 916 views
ADD COMMENT
0
Entering edit mode
3.4 years ago
ATpoint 65k

Does it mean that I can only do somatic calling for one pair of normal-tumor sequence each time

Yes, from what I understand VarScan2 is not designed for joint calling.

ADD COMMENT
0
Entering edit mode

Hi, Appreciate for your help. I just realized that I did not have a good understanding of somatic calling. And I found that in the GATK website it explained that "there is currently no way to perform joint calling for somatic variant discovery". Maybe I can find another method to put the snvs together to build the tree.

ADD REPLY

Login before adding your answer.

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