This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to tell if RNA-Seq Data has already been trimmed

I was wondering if there was an easy way to tell if this RNA-Seq data has already been trimmed or not? I got this data second hand so I'm not entirely sure what was done to it. It looks like the adapters were removed from each read? I'm not entirely sure though. I will attach some data to look at. Does anyone have experience with this type thing?

@HISEQ:249:C9MM3ANXX:8:1101:1925:2235 1:N:0:CACCGG AGAGCACATAATATCATCGTTTACATTTNNNNNNNNNNNNTAAATGTATATTGCACTTTACTGCGCACGAGAATCAATCAATCTGGACCTACATCAGTGTCACAAGGACTCGAAGAATTTTGACG + BBB@BFGGFGGGGGGGGGGGGGGGGGGG############===BFFFGCGGGGGDGGGGGGGGG<ggggfeggggggggggggg:cf0dgggggfgdfeggggggggggggdgccfgfg;cfggf @hiseq:249:c9mm3anxx:8:1101:2258:2233="" 1:n:0:caccgg="" gcttgttgtagtgagctagacgagacgnnnnnnnnnnnnnncgctcgaaaatatcatttacaaagctgttcattatgctcatcgccttcgacgaaattccggtgtcaggatggacctgcttgaga="" +="" bbbbcfggggggggggggggggggggg##############0="=FGGGGBGGGGGGGGGGGGGGGFGGGDGGGGGGGGGGGGGGGGGGGBCBGGGGEGGGGGGGGGGGEGGGGEGGEGGGGGGGD" @hiseq:249:c9mm3anxx:8:1101:2349:2235="" 1:n:0:caccgg="" tgatcacagagctcttcatgctctccgcnnnnnnnnnnnngtcggaggagtatcgctcctgtggtcccctagagctgctcttagctcgccgtgccggcgatgccaagctgtccatagacagattg="" +="" bb@bbgggggggggggggggfggggggg############="==CFGGGGGDE1FGGGGGFGDGGCFBFGGBGC@FGGGGGGGGEGGEGGGGGGGGEC9G.CGGGGGGGGGGGGGGGGGB@EGGGC" @hiseq:249:c9mm3anxx:8:1101:2602:2232="" 1:n:0:caccgg="" ttctggtgcagggagtcaatgttgtagnnnnnnnnnnnnnnnttgttggctctatacagaagctgactgtgcttgggatttcgcaagagggcctccacaatcgtcttgctcctcgcccgcattgc="" +="" bbbccgggggggggggggggggggggg###############="==EGGFGGGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG" @hiseq:249:c9mm3anxx:8:1101:2536:2234="" 1:n:0:caccgg="" gtgacatgccctccaagctatctggaatnnnnnnnnnnnnttagcagcgcacacatgcatgtttgtatatgcttgtatgtatgtatagcaaatatataattttagaacgaagaggagatcggaag="" +="" cccccgggggggggggggggggggggcg############="==FGGGGGGGGGGGGGGGGGFG">FGGGGGGGGGGGGGGGGGGGGGGGGFGGGGGGGGGGFGGGGGGGGBGGGGGDGGEGGGGG> @HISEQ:249:C9MM3ANXX:8:1101:2961:2232 1:N:0:CACCGG TATGGACCTTCGGTCTGGGCCAGTTCGNNNNNNNNNNNNNNNGCGCCATAGTAAATCGTTTCGAATATCTGCTGATTCAAGAGACCGGCTTCTTCGCTTTCGTACGGGAAACGCATGAGAATCAA + BCCBBGGGGGGGGGGGGGGGGGGGGGG###############===FGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGDGGEEGGEGGG @HISEQ:249:C9MM3ANXX:8:1101:3751:2240 1:N:0:CACCGG CATAACACCCGTATTGTACTTTTCCACGTGCNNNNNNNNNAAGTTGAGAAATACGTTCCAGATCACAAAGAACATACGAATCGGTGGCAAGTCGTGTGTGCCGAATAGGGAAAATAGGTAAATGG + CBCBBGGGGGGGGGGGGGGGGGGGGGGGGGG#########=<:C@GGGGGFGGGGEGGGGGGGGGGGGGBGGGGGGGGGGGGGGGGGGGGGEFDFGGGGGGGGGGGGGGGGGGFGGGBEGGG@GG @HISEQ:249:C9MM3ANXX:8:1101:5145:2235 1:N:0:CACCGG CTTGCGGATGGTCTTCACAGTAAGGTTCANNNNNNNNNNNCAGTTAATGGTGCACCCCGTGCACTTGTAGATCTCTGGGCCTTCAAATGCGAACGGATCGTTGGGATCCACGGTGGACTTCAAAA

rna-seq alignment genome

1 answer

Try running the sequences through FastQC. It's a great quality control program that will inform you about the presence of adapter or primer contamination. It also has a lot of other metrics, such as examining the quality of base-calling with respect to the length of the read.

Log in to answer this question.