Hi, I want to get information about the amount and average length of Poly-A tails in my RNAseq reads. Is there any tool that performs this? Preferably in a time efficient manner, because I have around a billion 100bp SE-reads. If I happen to find a lot of Poly-A in my reads, how can I best trim them?
I looked at RNA-SeQC but it doesn't seem to give any statistics about Poly-A tails.
Stroehli
•
written
by
Stroehli •
1 answer
If a read is all or almost all poly-A then it probably won't align anyway. If you're OK with not going through the entirety of the dataset, fastQC will let you know how much A (or AAAAA, since I think it uses pentamers) enrichment you have in your reads. I've not personally had any issues with too much polyA in the reads, but you could always use any read trimmer and just tell it to trim AAAAAAA, or something like that, from the 3' end.
Log in to answer this question.
-
Are the reads in the GTEx BAM/FASTQ files trimmed?
written by Andrew •
-
Did I do this right? Testing efficiency of poly(A) tailing.
written by Morgan S. •
-
Would you remove rRNA reads in silico before testing for diff. expression?
written by mschmid
-
Question about rRNA in data
written by joerodger2017 •
-
How does the calculated inner-distance between reads change after trimming (in relation to insert s…
written by jomaco
-
Differences In Alignments When Run Twice (Using The Same Reads/Reference As Input)
written by Stroehli •
-
Simulation Of Rnaseq Read Mapping With Reads Having Different Polya Tail Lengths.
written by Gaima Kushwaha •
-
Removal Of Low Quality Reads With Trailing "B"
written by Bioscientist
-
Abi Solid Transcriptome Assembly
written by Damian Kao
-
Fastq Format Problem When Using Bwa Aln -Q Option (Read Trim)
written by Bioscientist
