Trimmamotic out with 4 files
0
0
Entering edit mode
6 weeks ago
Saeed • 0

Dear all,

I use trimmomatic to QC for my RNA seq data and I got 4 outputs. I was wondering how can I put all these 4 files for alignment step?. I am going to use Hisat2.

RNA Hisat2 seq trimmomatics • 205 views
0
Entering edit mode

Most programs do not allow you to use unpaired fastq files that remain after trimming along with properly paired files in the same alignment run. You will need to align the files with singleton reads individually.

bbwrap.sh from BBMap can do this:

BBWrap is a simple wrapper that allows BBMap to be run multiple times without reloading the index each time. So, it can save some compute resources (particularly with a small number of reads and large reference), and is also handy for things like mapping paired and unpaired reads to the same reference, then outputting them in the same file.