Dear all,
We performed an RNAseq experiment with 50bp single (unpaired) reads. The analysis was performed using CASAVA pipeline.
I would like to do some QC on the resulting bam-file. However, when I try using Picard (for example CollectRNASeqMetrics) it gives me the following error:
SAM validation error: ERROR: Record 1, Read name HWI-ST571_152:6:1204:18235:6093, First of pair flag should not be set for unpaired read.
I searched this forum, but I didn't really find out if and how I can change this FLAG so that Picard won't throw an error anymore and quits.
Many thanks
Lien
We performed an RNAseq experiment with 50bp single (unpaired) reads. The analysis was performed using CASAVA pipeline.
I would like to do some QC on the resulting bam-file. However, when I try using Picard (for example CollectRNASeqMetrics) it gives me the following error:
SAM validation error: ERROR: Record 1, Read name HWI-ST571_152:6:1204:18235:6093, First of pair flag should not be set for unpaired read.
I searched this forum, but I didn't really find out if and how I can change this FLAG so that Picard won't throw an error anymore and quits.
Many thanks
Lien
Comment