I have RNA seq data for total yeast RNA (not poly-A purified). 90% of my reads should be ribosomal RNA. When aligning with bwa and then looking at read depth at the rDNA using bedtools' genomeCoverageBed or IGV, I see little to no coverage of the rDNA - certainly far far less than I should be seeing. If I go into the original fastq file and just pull out sequences at random and blast them to the yeast genome, almost all of them map to the rDNA, as they should. I imagine that when bwa sees greater than a certain read depth that it just tosses those reads. (I can't think of another way to explain this.) Looking at the help pages for bwa and gsnap (I'd be happy to use either), I don't see an explanation of this nor an option to turn it off. Can anyone suggest what is going on and what I can do to fix it? I want to see all of the reads.
Thanks.
Eric
Thanks.
Eric
Comment