search

nf-core / smrnaseq

A small-RNA sequencing analysis pipeline
https://nf-co.re/smrnaseq
MIT License
74 stars 125 forks source link

Paired-end data? #475

Closed anastasiaprime closed 1 month ago

anastasiaprime commented 1 month ago

Hello

Does pipeline work with paired-end data? Or only SE?

bounlu commented 1 month ago

It currently works with SE data only. I suggest to merge PE reads with a tool like NGmerge or fastq-join and then run the pipeline with the merged reads.

apeltzer commented 1 month ago

I would even suggest to only use R1 and ditch R2 entirely. It's small rna and thus covered by your R1 in its entirety ... no need to merge anything there.

bounlu commented 1 month ago

It depends on your read length and the small RNA species you are looking for. If you have PE50 or even PE75 and you want to quantify tRNA, snRNA, snoRNA, etc., then you better use both reads. If you are only interested in miRNAs, then R1 would be sufficient.

anastasiaprime commented 1 month ago

Thank you @apeltzer @bounlu , I appreciate that! Indeed, I am processing microRNA data, but I have paired 50 bp reads from the lab. And it seems to me that there is no microRNA in the reverse (R2) reads at all, so I will work only with R1.