[bam_header_read] EOF marker is absent. The input is probably truncated.

Question: [bam_header_read] EOF marker is absent. The input is probably truncated.

3.0 years ago by

Mic • 70

Australia

Mic • 70 wrote:

Hello, I did the following steps:

1. Prinseq 
2. Sickle
3. BWA for Illumina and 
4. SAM to BAM
5. MarkDuplicates 

After MakrDuplicates, I got the following error:

[bam_header_read] EOF marker is absent. The input is probably truncated.
[bam_header_read] invalid BAM binary header (this is not a BAM file).
[bam_index_core] Invalid BAM header.[bam_index_build2] fail to index the BAM file.
Fatal error: Exit code 1 ()

Picked up _JAVA_OPTIONS: -Xmx2048m -Xms256m
Ignoring SAM validation error: ERROR: Record 579, Read name NS500334:41:HLLLKBGXX:3:21601:2035:8235, MAPQ should be 0 for unmapped read.
Exception in thread "main" htsjdk.samtools.SAMException: Value was put into PairInfoMap more than once.  370: C_BATCH2:NS500334:41:HLLLKBGXX:1:21102:18509:10014
	at htsjdk.samtools.CoordinateSortedPairInfoMap.ensureSequenceLoaded(CoordinateSortedPairInfoMap.java:132)
	at htsjdk.samtools.CoordinateSortedPairInfoMap.remove(CoordinateSortedPairInfoMap.java:86)
	at picard.sam.markduplicates.util.DiskBasedReadEndsForMarkDuplicatesMap.remove(DiskBasedReadEndsForMarkDuplicatesMap.java:61)
	at picard.sam.markduplicates.MarkDuplicates.buildSortedReadEndLists(MarkDuplicates.java:285)
	at picard.sam.markduplicates.MarkDuplicates.doWork(MarkDuplicates.java:114)
	at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:206)
	at picard.cmdline.PicardCommandLine.instanceMain(PicardCommandLine.java:95)
	at picard.cmdline.PicardCommandLine.main(PicardCommandLine.java:105)
[bam_header_read] bgzf_check_EOF: Invalid argument

I was able to fix it with samtools view -h Galaxy62-\[SAM-to-BAM_on_data_10_and_data_54__converted_BAM\].bam | grep -v null | samtools view -bS - > Galaxy62.bam

I just wonder whether it is possible to do it anyhow to fix it in Galaxy or is it possible to fix the FASTQ files after PrinSeq and Sickle?

Thank you in advance.

Mic

picard bwa markduplicates samtools • 4.0k views

ADD COMMENT • link •

modified 3.0 years ago by Jennifer Hillman Jackson ♦ 25k • written 3.0 years ago by Mic • 70

3.0 years ago by

Jennifer Hillman Jackson ♦ 25k

United States

Jennifer Hillman Jackson ♦ 25k wrote:

Hello,

Are you working in Galaxy? This forum focuses on that usage. If running tools line-command, help from the 3rd party tools (either manual or dedicated help lists) and several other common bioinformatics forums could be a better troubleshooting resource (http://biostars.org, http://seqanswers.com, http://scicomp.stackexchange.com/).

If you loaded to Galaxy mid-analysis, then double checking that upload was successful is a good place to start when a tool runs into a malformed input: https://wiki.galaxyproject.org/Support#Loading_data

Should a header be removed during a step with Galaxy (would not be expected when performing an operation on a BAM dataset, but can occur with a SAM dataset), it can be replaced with the tool Picard: ReplaceSamHeader replace header in a SAM/BAM dataset.

Thanks, Jen, Galaxy team

ADD COMMENT • link written 3.0 years ago by Jennifer Hillman Jackson ♦ 25k

Similar posts • Search »