Annotating big VCF file hg38
Hi , Please I need to annotate a vcf file ( phased sequencing data of about 1000 persons ). Its a hg 38 built , I tried with this file common_all_20180418.vcf.gz and this command :
bcftools annotate -a common_all_20180418.vcf.gz -c ID mydata.phased.call.allchr.vcf.gz --output-type z -o mydata_dbSNP151.vcf.gz
and i get this error message
[W::hts_idx_load2] The index file is older than the data file: mydata.phased.call.allchr.vcf.gz.tbi
[W::hts_idx_load2] The index file is older than the data file: mydata.phased.call.allchr.vcf.gz.tbi
[E::bgzf_flush] File write failed (wrong size)
[E::vcf_parse_format] Number of columns at chr14:77174450 does not match the number of samples (60 vs 929)
[E::bgzf_flush] File write failed (wrong size)
Do you have any idea how to solve this ???!!
Thank you
1 answer
Log in to answer this question.
More posts like this
-
SAMtools error closing standard output: -1
written by theadrijasaha •
-
bgzip error 4
written by yb •
-
Error with Tabix file
written by bioinformatics.girl •
-
Use RSEM and Bowtie2 to align paired-end sequences
written by melissachua90 •
-
File write failed (wrong size) when filtering a bam file with interset
written by ManuelDB
-
Failed to open plain.vcf: unknown file type
written by obidobi •
-
Vg gamsort temporary folder
written by arunkumarramrsh •
-
suppress pysam index warnings.
written by b10hazard •
-
Samtools sort error with a bam file
written by Vasu
-
Error with samtools - indexing alignments: Writing to standard output failed: broken pipe - python
written by coyot001 •
