This is a test version of Biostars. For the public version, visit https://www.biostars.org.
SNPs not filtered from a GATK vcf file

Hi All,

I have a vcf file with SNP information from multiple samples made using GATK:

gatk --java-options "-Xmx4G" HaplotypeCaller -R ref.fa -I bams.list -L ch01 -O 01.vcf

Individual vcf's were made chromosome wise and then concatenated:

bcftools concat -o merge.vcf 01.vcf 02.vcf 03.vcf 04.vcf 05.vcf

I want to keep the SNPs only for the final vcf file so I did:

bcftools filter -i 'TYPE="snp"' merge.vcf > merge_SNP.vcf

But the output file still has INDELS. Then I tried using bcftools view for the same job:

bcftools view -v snps merge.vcf > merge_SNP.vcf

The output file again has variants other than SNPs. I am not sure what is going wrong. I will appreciate any suggestions. Thank you!

snp
• 1,509 views • link
updated by inedraylig • inedraylig • written by evelyn evelyn

You can also use SelectVariants module from GATK

• 0 views • link
written by svp svp

1 answer

The recommended way to filter indels would be to use --exclude types:

bcftools view --exclude-types indels merge.vcf > merge_SNP.vcf

assuming that you only have SNPs and INDELs in your vcf file. bcftools filters using the INFO field, so you can look at your vcf and see where the identity of the call (SNP/INDEL) appears and if it's indeed in the INFO field.

• 0 views • link
written by inedraylig • inedraylig

Log in to answer this question.

More posts like this
  • variant calling with platypus
    written by evelyn evelyn
    Hello Everyone, I used Platypus for joint variant calling using: python Platypus.py callVariants --bamFiles=bams.list --regions=ch01 --refFile=ref.fa output=out1.vcf I got vcf files for each chromosome and …
  • bcftools multiple samples
    written by evelyn evelyn
    Hello, I have downsampled sorted bam files using: samtools view -h -s 0.5 E.sorted.bam > E.5x.sorted.sam And again compressed the files: samtools -bS E.5x.sorted.sam > …
  • GATK multiple samples
    written by evelyn evelyn
    Dear all, I am trying to call SNPs for multiple sorted bam files altogether using GATK. I am using this code: INPUT_DIR=path REF_DIR=path INPUT=$(ls -1 …
  • how to merge gvcf using bcftools
    written by evelyn evelyn
    Hello Everyone, I created gvcf files using bcftools: bcftools mpileup -Ov --gvcf 0 -f ref.fa example.sorted.bam | bcftools call -m --gvcf 0 -o example.vcf I …
  • gvcf error bcftools
    written by evelyn evelyn
    Hello, I used `bwa` for alignment and I am using `gvcf` for calling variants with this code: bcftools mpileup -Ov --gvcf 5 -f ref.fa example.sorted.bam …
  • non-SNP variants in a vcf merge file
    written by evelyn evelyn
    Hello, I have used `bwa` for alignment and `platypus` for variant calling and filtered only `SNPs` using `bcftools`. I then merged the vcf files with …
  • problem with vcf file by gatk
    written by evelyn evelyn
    Hello, I am making `vcf` files using `gatk` with following codes: java -jar $PICARD_DIR/picard.jar ValidateSamFile \ I=example.sorted.bam \ MODE=SUMMARY java -jar $PICARD_DIR/picard.jar AddOrReplaceReadGroups \ I=example.sorted.bam …
  • vcf from sorted bam
    written by evelyn evelyn
    I am using `bcftools` to make `vcf` from `sorted bam` using: bcftools mpileup -Ov -f ref.fa sorted.bam | bcftools call -mv -o sample.vcf But I …
  • samtools SNP calling
    written by evelyn evelyn
    I am using `samtools` for SNP calling from `.sorted.bam` file: bcftools mpileup -f ref.fa bwa.sorted.bam | bcftools view -Ov - > bwa.bcftools.vcf VCF file gives …
  • Remove duplicate SNPs only based on SNP ID in bcftools
    written by janhuang.cn janhuang.cn
    I want to remove duplicate SNPs only according to SNP ID (e.g. rs123), from a vcf file. I used the below command, but it duplicate …