Question: gatk SelectVariants odd behavior
0
gravatar for jon.klonowski
5 weeks ago by
jon.klonowski10 wrote:

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 modified 21 days ago by Biostar ♦♦ 20 • written 5 weeks ago by jon.klonowski10

You may want to post on the GATK Support Forum.

ADD REPLYlink written 4 weeks ago by Kevin Blighe42k
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: 1753 users visited in the last hour