But I download my bamfile from this link: "ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/data" I think that problem is from other thing. I work with some variant somtool to create sorted file for using in HAPCUT, but error occurred every time.Is it possible that samtool do not work correctly? or my error is from other thing? thanks,
Error: Eof Marker Is Absent When Processing A Bam File
Hi
I work with HAPCUT software to make Haplotypes using a sorted bam file and a vcf file containing variants. When I run this command,
./extract_hairs --VCF variantfile.VCF --bam alignedreads.sorted.bam --maxIS 600 > fragment_matrix_file
an error says that
[bam_header_read] EOF marker is absent
Also before that, when I tried to sort bam file with samtools this error occurred. but samtools continued to make sorted bam and then I used it in HAPCUT.
Thanks for any advice on this problem.
2 answers
You have an invalid bamfile. Recreate it.
• 161 views
•
link
written
by
Istvan Albert
If samtools raises an error then it means that your bamfile is incorrect. There are many ways a download can go bad, give us the link to the actual bamfile and perhaps someone can test it for you on their system.
• 0 views
•
link
written
by
Istvan Albert
The actual link that i downloaded my bamfile is: " ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/data/HG00096/alignment/HG00096.chrom11.ILLUMINA.bwa.GBR.low_coverage.20111114.bam ".
I work on chromosome 11, if this error is from my file please tell me another link for download bamfile.
If hard drive was full when made bam file it would stop but not give any indication so maybe the person that made the bam file did not notice that during the operation the hard drive became full. Try to get the sam file and convert yourself or use a different bam file to check it is the bam file which is not complete. I recently created a bam file but my hard drive was full so had this issue.
Even though the file will be sorted using samtools it is not the complete bam file so your results will not be based upon the full data set.
• 144 views
•
link
updated
by
Ram
•
written
by
rob234king
•
written
by
rob234king
Log in to answer this question.
More posts like this
-
samtools mpileup fail to create bcf
written by mthm •
-
getting error when indexing bam files
written by Sara
-
Errors with samtools depth on merged bam file
written by rbronste
-
Samtools error sorting
written by Msadeghibiology91 •
-
Problem on indexing of sorted BAM file
written by Gaurab2017 •
-
bam-readcount error bam_header_read
written by fangqiong.ling •
-
sam to bam conversion
written by sukesh1411 •
-
BAM header file error
written by Bioblazer
-
Using samtools and Picard tools with cram files
written by Chris Gene
-
Filtering nth length from mapped SAM/BAM file.
written by unique379
HAPCUT is interesting. How long are the haplotypes it generates from paired-end sequencing data?
If you get this error. You can maybe ask yourself if the bam file you created was created in a perfect manner. If I get this error, I'll always look at the sam file, and convert it to a bam file again. Don't run things in parallel, since this can result in errors. After this conversion, sort your bam file and then perform indexing again. Normally you won't run into this problem afterwards.
A good pipeline converts alignment stream to BAM right away, with or without multithreading does not matter. If the file is corrupt then because of memory shortages or corruption during the alignment.