After running some RNA-seq experiment, we have observed an over-represented sequence which do not correspond to any adapter used. It clearly seems to be a contaminant.
But my main question is why do I get this sequence only in one strand?
As I get the sequence in 2 400 000 reads, I don't get its reverse-complemented sequence. I don't know why, because if I get a contaminant sequence, I should get half of the reads with the sequence, and half with its complement, shouldn't I?
Have you tried blasting it?
Yes, it's 60S Rrna human, so a contamination is obvious. But what i can't explain is why is there only this: R2adapter - Sequence of 300 bases - R1adapter and none of this: R1adapter - Sequence of 300 bases - R2adapter The adapters must be set in a random way, no?
May be first you map your reads against rRNA and tRNA sequences (preferably via bowtie), and discard the mapable reads. If remaining reads still show over-represented sequences then there is something weird, it could also be the case of chimeric transcripts.
Actually, i've already done this mapping, and i know it's form this sequence: Chain 5, Structure Of The H. Sapiens 60s Rrna, 5070 bp. The problem is not what the sequence is, but why is it represented only in the 5'-3' strand, and not in the other way. I mean, if adapters are setting randomly around the rna fragment, why do i only observe R2adapter-Sequence of 300 bases-R1adapter, and not R1adapter-Sequence of 300 bases-R2adapter?
Which library prep method was used for the RNA-seq?
Script Seq, Illumina
Thanks for your answers/time!
So it's strand specific and rRNAs are supposed to be removed.
So if it is strand specific, it means that I will observe always get the R2adapter in 5' of my sequence and never in 3' ? So this: R1adapter-Sequence of 300 bases-R2adapter can't exist ?
Yes, from the FAQ of Script-seq:
You should not have RXadapter - Sequence of 300 bases - RYadapter in first place, as the first adapter will be sequenced by the primer anyway. Have you checked with the technician who did the experiment what went wrong? something seems very fishy on wet-lab level!