Dear all,
I want to remove reads that map to ribosomal RNAs. I tried to modify the sam file (output of TopHat), and remove the reads that correspond to rRNAs. However, I could not run cufflinks on the modified sam files. Then I tried to sort them by converting them to bam (using samtools), but even then, it still did not work. I got the following error messages from cuffdiff.
[22:39:41] Inspecting maps and determining fragment length distributions.
SAM error on line 66185: CIGAR op has zero length
SAM error on line 94994: CIGAR op has zero length
SAM error on line 135374: CIGAR op has zero length
Has anyone seen this error? What does it mean?
What would you recommend to remove ribosomal reads?
thanks,
Grace
I want to remove reads that map to ribosomal RNAs. I tried to modify the sam file (output of TopHat), and remove the reads that correspond to rRNAs. However, I could not run cufflinks on the modified sam files. Then I tried to sort them by converting them to bam (using samtools), but even then, it still did not work. I got the following error messages from cuffdiff.
[22:39:41] Inspecting maps and determining fragment length distributions.
SAM error on line 66185: CIGAR op has zero length
SAM error on line 94994: CIGAR op has zero length
SAM error on line 135374: CIGAR op has zero length
Has anyone seen this error? What does it mean?
What would you recommend to remove ribosomal reads?
thanks,
Grace