paired reads have different names BWA MEM after samtools bam > fastq
0
0
Entering edit mode
6 weeks ago

I have used bwa mem to align to a host genome and output the unmapped reads. I have then sorted this resulting BAM and split into pairs fastq files. Whya fter sorting the BAM file am I still getting paired reads have different names error from bwa mem and can I solve this without resorting the fastq file itself?

alignment • 230 views
ADD COMMENT
0
Entering edit mode

R1 can align, R2 cannot, the new fastq files would miss the R1 name for that read and boom -- out of sync. You probably want to resync/repair it with something like repair.sh from BBMap suite.

ADD REPLY
0
Entering edit mode

So there is no way to fix this during the conversion?

I have read that such repair tools can lead to data loss, is this not the case?

ADD REPLY

Login before adding your answer.

Traffic: 2144 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6