We've moved!
This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!

Error in VariantAnnotator

hi
Because my reference is in scaffold level and GATK not designed to analyse such referenes, thus my variant caller project done using samtools. Now i want to add QD(qualityofdepth) to vcf files using VariantAnnotator but i get error as following:

ERROR MESSAGE: Invalid command line: Argument annotation has a bad value: Annotation QD was not found; please check that you have specified the annotation name correctly

my used command:

java -jar GenomeAnalysisTK.jar -T VariantAnnotator -R /home/mylinux/refdic/DROMREF.fasta -I /home/mylinux/Drom1/realigned_D1.bam -V /home/mylinux/Drom1/variants/SNPSwhitGQ.vcf -L /home/mylinux/Drom1/variants/SNPSwhitGQ.vcf -A QD -o /home/mylinux/Drom1/variants/GATKANN.vcf

how can i add QD to vcf file?

Answers

  • rezakhalkhali29rezakhalkhali29 IranMember

    i fixed above mentioned error by using -A QualByDepth but i get another error:

    ERROR MESSAGE: Bad input: We encountered a non-standard non-IUPAC base in the provided reference: '13'

    i used that reference for local realigning in GATK and it worked well but now it give error. what happend?

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie admin

    That typically happens when your reference file was edited or saved on a Windows machine because of character encoding differences. See previous threads on this topic for help with resolution.

Sign In or Register to comment.