Quantcast
Channel: vcf — GATK-Forum
Viewing all articles
Browse latest Browse all 624

How to select samples that are polymorphic on a specific locus from a joint genotyped vcf file?

$
0
0

Hi, I am trying to select the samples that are polymorphic on a specific locus from a joint genotyped vcf file using SelectVariants tool and JEXL expressions with no success. The command I am trying to use is:
gatk SelectVariants \ --reference "/home/hg19.fa" \ --variant "/home/recalibrated_snps_recalibrated_indels_VQSR.vcf" \ --select 'POS == '11240190'' \ --select 'vc.getGenotype().isHet()' \ --output "/home/isHet.vcf"
Selecting on the position alone seems to work fine, but I want only the samples that are Het for this exact position, any suggestions please?
Thanks in advance
Nawar


Viewing all articles
Browse latest Browse all 624

Latest Images

Trending Articles



Latest Images