I align forward and reverse fastq reads to a fasta reference using bwa, then sort the alignment file:
bwa mem -M -t 16 ref_galr2a_AA.fa E01_S49_L001_R1_001.fastq.gz E01_S49_L001_R2_001.fastq.gz | samtools sort > test.bam
Now say I want to convert back this bam to fastq files; I am used to doing:
bedtools bamtofastq -i test.bam -fq forward.fastq -fq2 reverse.fastq
Then put the result in a gunzip
It creates file
On macOS, I can typically double-click a file like this, and I'll get the fastq file inside (it uses Archive Utility by default).
However, with this one I get:
Unable to expand "forward.fastq.gz" into "folder".
(Error 79 - Inappropriate file type or format.)
This does not happen if I unzip/put back in a gzip the original fastq directly.
i.e. say I unzip E01_S49_L001_R1_001.fastq.gz by double-clicking it (Archive Utility), then run:
It generates back E01_S49_L001_R1_001.fastq.gz, which I am able to open with Archive Utility.
So that tells me bedtools bamtofastq is not innocent? What is going on?