A question on GetPileupSummaries (GATK)
0
0
Entering edit mode
5.1 years ago
Laven9 • 0

I am confusing what "common_biallelic.vcf.gz" is? Am I right if I run it like this:

 gatk SelectVariants \
 -R Homo_sapiens_assembly38.fasta \
 -V gnomad.vcf.gz \
 --restrict-alleles-to biallelic \
 -O common_biallelic.vcf.gz

Or should I change the -V for something else like my normal.vcf?

GetPileupSummaries • 2.4k views
ADD COMMENT
0
Entering edit mode

Hey, I don't think you should use your normal.vcf. Here dbSNP site has common_all_xxx.vcf.gz , also I think GATK bundle will provide common variant file.

ADD REPLY
0
Entering edit mode

Thank you, MatthewP ! But as I am new in the field, could you please tell me how to interpret the file common_all_xxx.vcf.gz? If it is right that in my opinion, it should be a file containing all the common snps in a large group of people, so that gatk4 could use it to exclude those snps in the Mutect2.output.file?

BTW, sorry to use the answer section, my "add reply" is out of work.

ADD REPLY

Login before adding your answer.

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