Regarding the error message you got, there is some discussion at:http://seqanswers.com/forums/showthread.php?t=3551
If you moving out the reads (from sam file) corresponding to rRNA, you can remove them from your original fastq file (using read IDs) and then re-run the alignment.
-
I just looked at the coordinates of each SAM output, and removed lines that overlapped with rRNA locations, as I read in each line.
I reasoned that since my SAM file was sorted before editing, it should have remained sorted if I just moved certain lines.Leave a comment:
-
How did you removed the reads from SAM which correspond to rRNA?Leave a comment:
-
Editing SAM files
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
-
Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.
The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
...06-02-2026, 10:05 AM -
With the launch of new single-cell sequencing platforms in 2026, the field stands at an exciting inflection point. This article surveys the most impactful advances in the field and discusses how they’re reshaping research in cancer, immunology, and beyond.
Introduction
Single-cell sequencing technologies have undergone remarkable advances over the past decade, transitioning from low-throughput experimental approaches to highly scalable platforms capable of...05-22-2026, 06:42 AM
Leave a comment: