Biostar Beta. Not for public use.
Usage of Mutect2
0
Entering edit mode
14 months ago

Hi,

I have been trying to use the Mutect2 of GATK toolkit, and noticed an example usage:

java -jar GenomeAnalysisTK.jar \
-T MuTect2 \
-R reference.fasta \
-I:tumor tumor.bam \
-I:normal normal.bam \
 [--dbsnp dbSNP.vcf] \
 [--cosmic COSMIC.vcf] \
 [-L targets.interval_list] \
 -o output.vcf

Here, I am assuming that the reference.fasta is the genome sequence in fasta format, and tumor.bam is the file containing details about the mapped reads. I am not sure what is normal.bam. Is it the bam file containing details of the mapped reads of a reference transcriptome with the reference genome?

ADD COMMENTlink
1
Entering edit mode
16 months ago
MatthewP • 110
China

Hello, normal here means normal tissue. People need to compare variants from tumor tissue and normal issue(example: PBMC) to find some variants new(somatic) in tumor tissue.

ADD COMMENTlink
0
Entering edit mode

Since I do not have the data for any normal tissue from the same human, is there a way to skip this step?

ADD REPLYlink
1
Entering edit mode

Maybe you should try germline pipeline which use HaplotypeCaller instead of Mutect2 .

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3