how to extract the RD, AD, VAF information for the normal and tumor samples from the MuTect2 output vcf file?

Dear all:

I have the tumor-normal pair Mutect2 output VCF file. I want to get the DP, RD, AD and VAF information for both tumor and normal samples. Can anyone please give me some suggestions? Thanks very much.

best

exom-seq mutect2

• 3,294 views • link

updated 9.0 years ago by Biostar

• written 9.0 years ago by lhaiyan3

Please post example VCF records. In the mean time, refer to bcftools or snpsift manuals to extract info and format information from VCF.

• 0 views • link

written 9.0 years ago by cpad0112

0 answers

No answers yet.

More posts like this

filtering variants in a Strelka2 VCF file based on AD and AF
written 2.9 years ago by Bogdan

Dear all, I would appreciate having your suggestions on the following. I am working with a VCF file that was produced by Strelka on Tumor-Normal …
Which parameter can be used to set up cutoff not to call peaks <150bp for the macs2 broad peak call…
written 7.4 years ago by lhaiyan3

Dear all, I don't want to call peaks <150bp in the `macs2` broadpeak calling. Can you please tell me some suggestions which parameter I can …
Filtering VCF with python
written 8.4 years ago by miss •

Hi, I extracted informations that I need from my vcf file and one line of it looks like this: X 126165205 . T G DCAF12L2 …
ngsplot-how to get multiplot around the IP peak region and the genes nearest to peak region?
written 8.4 years ago by lhaiyan3

Dear all, I have the 3 IP (TF1, TF2, TF3) data, and I use mac2 callpeak to get the peak files (TF1.peak). I want to …
How to use GSVA for the RNA-seq RPKM file
written 8.6 years ago by lhaiyan3

Dear all, I am trying to use GSVA for the RNA-seq analysis. I have the refseq annotation based RPKM files for the RNA-seq. RPKM <- …
How to analyze the targeted RNA-panel data?
written 9.1 years ago by lhaiyan3

Hi, We have the targeted RNA-seq fastq file for 9 genes. The first 12 bp of any given read in the FASTQ file is a …
How to convert the bam file to the aligned fasta file and perform the phylogenetic analysis?
written 9.6 years ago by lhaiyan3

Dear all, We have a bunch of bam file for the targeted tumor genes sequencing. We want to perform the phylogenetic analysis for these data. …
How to build the STAR Genome index from the .gbk file?
written 9.8 years ago by lhaiyan3

Dear all, I need to build the STAR alignment Genome index file for my analysis. I have the .GBK file for my target sequence. I …
ChrM in the Varscan2 copynumber output
written 10.4 years ago by lhaiyan3

Hi, I am using the varscan2 to call the copynumber for the WES. But my output.copynumber file just has the chrM call. I used the …
The plot generated by sciClone
written 10.4 years ago by lhaiyan3

Hi, I am using the sciClone for the tumor clonality analysis. I got the plot, but my plot is very clean. It did not show …