This is a test version of Biostars. For the public version, visit https://www.biostars.org.
FastQC polyA and polyG trimming in RNA-seq

Hi!

I have raw Fastq data with 151bp paired-end reads, which were trimmed using Trim Galore (v0.6.10). After trimming, I observed overrepresented polyG sequences (purple line) in the FastQC results for Read 2, along with a small presence of polyA sequences (green line). Should I remove both the polyG and polyA sequences, only the polyG sequences, or leave both as they are? enter image description here

polya qc trimming rna-seq fastq
• 2,593 views • link
written by joonhong kwon • joonhong kwon

1 answer

Yes. Trim both before performing any downstream analysis. trimgalore uses cutadapt under the hood, so: https://cutadapt.readthedocs.io/en/v1.8.1/recipes.html#trimming-poly-a-tails

• 0 views • link
written by Mark Mark

just make sure the adapter overlap is long enough I'd say over 10 bases

Out of curiosity I checked that 6 base overlap, which typically suffices for artificial adapter sequences would not probably be too short for Gs,

6Gs in a row occur over 6 thousands times even jus the shortest chromosome 22 of the human genome

• 0 views • link
written by Istvan Albert Istvan Albert

Log in to answer this question.

More posts like this