GATK licensing moves to direct-through-Broad model -- read about it on the GATK blog

Select indels with lenght smaller than 10.

WimSWimS Posts: 25Member

How can I select indels with lenght smaller than 10 bp from a vcf file?

I tried

java -jar GenomeAnalysisTK.jar -T SelectVariants -R ref.fa --variant INDEL.vcf -o INDEL_maxLenght10.vcf -select 'vc.getIndelLengths().0 < 10'

but the output still contains all the Indels, also the ones larger than 10 bp.

Best Answer


Sign In or Register to comment.