Hello,
I have used bowtie2 to map some newly sequenced reads to the sheep reference genome.
After that with samtools I converted the sam file to bam file. I'm interested for the reads that are aligned in the mitochondria, so I extracted the particular region from the bam fie. So far everything seems ok, except the fact that the new bam file has "invalid BAM binary header". How can I keep the header to the new bam file?
Thank you very much in advance
Vasilis
Can you show the commands used? How did you extract the region of interest?
I used:
and after that I converted again into bam.
Just do
Yes, it works!
Thank you very much