This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Issue with cDNA PCR artifact in RNA-Seq dataset

I have a large RNA-Seq dataset that I am working on. In half of my samples approximately 2-5% of reads are duplicates of this sequence "AAGCAGTGGTATCAACGCAGAGTACTTTTTTTTTTTTTTTTTTTTTTTTTTT". I understand this sequence is a primer from the cDNA step of RNA-Seq data preperation. This leads to a clear divide between my samples, such that in the MDS space, one axis divides samples by type and the other divides samples by whether or not they have this artifact.

How should I deal with this? Can I completely remove this sequence from my files? How should I go about that?

pcr rna-seq artifact cdna

1 answer

You can remove such sequences with multiple tools: BBDuk for example. However, I do not think that the divide will emerge just because of that. Rather, it may be low-input/complexity samples that happen to suffer from this artifact the most.

Log in to answer this question.