Question: can VarScan2 do joint somatic calling?
0
gravatar for Zoey Zou
4 months ago by
Zoey Zou40
Australia
Zoey Zou40 wrote:

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 somatic calling varscan • 193 views
ADD COMMENTlink modified 4 months ago • written 4 months ago by Zoey Zou40
0
gravatar for ATpoint
4 months ago by
ATpoint22k
Germany
ATpoint22k wrote:

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 COMMENTlink written 4 months ago by ATpoint22k

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 REPLYlink written 4 months ago by Zoey Zou40
Please log in to add an answer.

Help
Access

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