Question: SortVCF in Picard
gravatar for Maheswari P
3.0 years ago by
Maheswari P0 wrote:

i am working on SORTVcf picard tool. java -jar picard.jar SortVcf I=filename.vcf o=output.vcf SEQUENCE_DICTIONARY=reference.dict I am getting error as Exception in thread "main" java.lang.IllegalStateException: Key found in VariantContext field INFO at chr1:10327 but this key isn't defined in the VCFHeader. We require all VCFs to have complete VCF headers by default. at htsjdk.variant.vcf.VCFEncoder.fieldIsMissingFromHeaderError( at htsjdk.variant.vcf.VCFEncoder.encode( at htsjdk.variant.vcf.VCFRecordCodec.encode( at htsjdk.variant.vcf.VCFRecordCodec.encode( at htsjdk.samtools.util.SortingCollection.spillToDisk( at htsjdk.samtools.util.SortingCollection.add( at picard.vcf.SortVcf.sortInputs( at picard.vcf.SortVcf.doWork( at picard.cmdline.CommandLineProgram.instanceMain( at picard.cmdline.PicardCommandLine.instanceMain( at picard.cmdline.PicardCommandLine.main(

I am working in 2.5.0 picard tool. Please any1 help me in this problem. I am beginner to this RNA seq analysis. Let know the early solution with reason of error.

rna-seq • 2.0k views
ADD COMMENTlink modified 23 months ago by Biostar ♦♦ 20 • written 3.0 years ago by Maheswari P0
gravatar for Pierre Lindenbaum
2.8 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum126k wrote:

The message is clear: your VCF is malformed , there is a field in the INFO column (eg 'AF', 'AC' ..) which definition is missing in the VCF header. You can add it with awk...

awk '/^#CHROM/ {printf("##INFO=<ID=AAC,Number=1,Type=Integer,Description=\"Alt Allele Count\">\n");} {print;}'  input.vcf > out.vcf
ADD COMMENTlink written 2.8 years ago by Pierre Lindenbaum126k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1338 users visited in the last hour