I'm trying to analyze some old miRNAseq data and I'm having a hard time getting the reads trimmed to the expected length. I ran FastQC on the fastq files and despite all the reads being the same length (55bp), they all had <5% adapter content.
I've tried trimming with CutAdapt using the adapter sequences from Illumina, and with different settings for CutAdapt. However, none of it worked, the majority of the reads were still ~55 bp. I also checked a few of the full length reads in mirBase to see if the data actually needed to be trimmed, and none of them yielded any hits. So, I'm assuming they do, but I could be wrong since this is my first time starting from fastq files for RNAseq data..
How do I figure out what adapter/sequence needs to be trimmed from my data or if it needs to be trimmed at all before alignment?
Other relevant info:
- These were sequenced ~10 years ago.
- The library was made with either the Qiagen or Illumina kits for small/miRNA.
- It's paired end.
- Expected read length is ~22 bp.
Thank you!