How to extract mutation from...?
0
0
Entering edit mode
6.0 years ago
ruhollah ▴ 10

I have a VCF file in which the genome of the normal tissue is compared with that from the tumor in the same sample.

Using this VCF file, how can I get all of the mutations in this sample? (I'm using R 3.4.4)

Any idea & comment, even if they seem very basic to you highly appreciated.

Cancer genome • 1.7k views
ADD COMMENT
1
Entering edit mode

Is this a technical question like how to load a VCF into R? Please elaborate.

ADD REPLY
1
Entering edit mode

I have a VCF file in which the genome of the normal tissue is compared with that from the tumor

Probably not, most commonly both are compared to the reference genome for variant calling, not to each other.

If you take a look at the vcf specifications you'll see that alleles are encoded as 0 (ref) or 1 (alt). If you want to use R then that's something that you can filter on in a daraframe.

ADD REPLY
1
Entering edit mode

The OP is likely performing somatic variant calling, since they are talking about tumors and normals.

You previously clarified the concept: Somatic vs Germline Variant Calling

ADD REPLY
1
Entering edit mode

Isn't this VCF actually a text file of the differences (mutations if you want to call them that) you are looking for?

ADD REPLY

Login before adding your answer.

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