easyRNASeq

[calmac3]

This dev package looks very promising to simplify read summaries. I'm sorry if I missed this in the documentation, but I am wondering how the package handles paired-end reads. I have a dataset with a mixture of paired- and single-end reads. How does this impact the "readLength" line?

Thanks

[calmac3]

I suppose I will answer my own question. At this point, without input from the developers, do not try to mix paired-and single-end reads.

[regyre]

Hi Calmac3,

The best place to get an answer about Bioconductor packages is the bioconductor mailing list. I honestly would not have expected posts about easyRNASeq here on SEQanswers.

To answer your question: if you use version 1.2.5 of easyRNASeq, there's no support for PE. If you use 1.3.14, you do not need to provide the readLength argument as this information will be extracted directly from your data. Concerning the PE data, they are still considered as SE as I'm working on a proper way to define/count fragments and not reads. I hadn't thought of mixed SE-PE setup but will keep that in mind. If you have other questions, or suggestions, comments, please post on the Bioc mailing list. At list you'll get a quicker answer!

Nico