Biostar Beta. Not for public use.
Question: gatk SelectVariants odd behavior
0
Entering edit mode

gatk SelectVariants is not working properly for me...

I have a file of SNP ids (i.e chr1:861219:G:C) delimited with /n and a vcf where the snp ID is also formatted chr:pos:ref:alt (chr1:861219:G:C).

The list file of SNPs were generated by subsetting a portion of the vcf SNPs (based on 1000G frequency) when I try:

gatk SelectVariants --variant File_name.vcf.gz -O Output_file_name.vcf.gz --keep-ids SNPS_of_interest.txt

I get a file of just the headers.

and then when I try:

 gatk SelectVariants --variant File_name.vcf.gz -O Output_file_name.vcf.gz --exclude-ids SNPS_of_interest.txt

The output is the same as the input excluding a zless | wc -l of ~10

ADD COMMENTlink 10 months ago jon.klonowski • 10 • updated 9 months ago Biostar 20
Entering edit mode
0

You may want to post on the GATK Support Forum.

ADD REPLYlink 10 months ago
Kevin Blighe
43k

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.0