Question: Mutect Vcf Output
0
gravatar for KS
5.4 years ago by
KS340
KS340 wrote:

Hello Everyone,

I just started using muTect to detect somatic mutations in matching tumor normal bam files. I would like to output in VCF format for downstream annovar annotation process. Can anyone suggest how to output MuTect in VCF format

Below is the command I am using currently:

java -Xmx2g -jar muTect-1.0.27783.jar --analysis_type MuTect --reference_sequence hg19.fasta --dbsnp dbsnp_132_b37.leftAligned.vcf --input_file:normal Normal.bam --input_file:tumor Tumor.bam --out example.call_stats.txt --coverage_file example.coverage.wig.txt

Thanks Suz

vcf • 6.3k views
ADD COMMENTlink modified 5.4 years ago by henryvuong720 • written 5.4 years ago by KS340

Is "--coverage_file example.coverage.wig.txt" input or output parameter ?

ADD REPLYlink written 3.8 years ago by Chirag Nepal2.1k
3
gravatar for henryvuong
5.4 years ago by
henryvuong720
USA
henryvuong720 wrote:

the latest version of MuTect (ver 1.1.4) downloaded from (http://www.broadinstitute.org/cancer/cga/mutect_download) can generate vcf output with --vcf parameter. hope it helps

ADD COMMENTlink written 5.4 years ago by henryvuong720

Hi,

I am trying to run the latest mutect ver 1.1.4 to generate the output in vcf format but i am running with error. The command am using is 

java -Xmx14g -jar /scratch/GT/softwares/mutect/muTect-1.1.4.jar --analysis_type MuTect --reference_sequence /scratch/GT/vdas/test_exome/exome/hg19.fa -B:cosmic,VCF /scratch/GT/softwares/mutect/hg1/hg19_cosmic_v54_120711.vcf -B:dbsnp,VCF /scratch/GT/vdas/test_exome/exome/databases/dbsnp_137.hg19.vcf --input_file:normal /scratch/GT/vdas/pietro/exome_seq/results/N_S8980/N_S8980.realigned.recal.bam --input_file:tumor /scratch/GT/vdas/pietro/exome_seq/results/T_S7998/T_S7998.realigned.recal.bam --vcf /scratch/GT/vdas/pietro/exome_seq/results/mutect/mutect_S_313soma2.vcf --coverage_file /scratch/GT/vdas/pietro/exome_seq/results/mutect/coverage_S_313soma2.wig.txt

Can you please tell me where am getting it wrong for getting a desired output of Mutect in vcf format?

ADD REPLYlink written 4.0 years ago by vchris_ngs4.5k

yes it is working fine now with the MuTect (ver 1.1.4)

ADD REPLYlink written 4.0 years ago by vchris_ngs4.5k

Hi there,

--coverage_file: This corresponds to input or output parameter? Considering we are inputting two bam files (N/T pair), how did you generate one wig file? Simply by merging two bam files?

 

--intervals:  What input did you give for this ? If i want to run the analysis on genome-wide.
ADD REPLYlink written 3.9 years ago by Chirag Nepal2.1k
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: 697 users visited in the last hour