error: BAM to FASTQ

Question

[mailard]

Hi! I run

samtools bam2fq bx_mapped.bam

and my result is

[E::sam_parse1] no SQ lines present in the header [bam2fq_mainloop] Failed to read bam record.
[bam2fq_mainloop] Error writing to FASTx files.: No such file or directory [M::bam2fq_mainloop] discarded 0 singletons
[M::bam2fq_mainloop] processed 0 reads

How can I fix it?

Thank you for your answer!

[colindaven]

Try remapping your BAM (which aligner did you use?). It looks to be corrupted/incorrect as the header does not contain reference chromosome lines.

samtools view -h BAMFILE | less

[mailard]

I use histat2 Thanks for your help, I decided my problem!