Hi all,
I am working with S. pombe doing 100bp single end Illumina HiSeq 2500 RNA-seq. Despite being told the trueseq stranded protocol was used one lot of my libraries appear unstranded (roughly 50/50 split between fwd and rev mapping reads in IGV). The second lot appears stranded as expected (>95% mapping bias to antisense strand in IGV). I'm trying to understand the best way to do read summarisation for these libraries using featurecounts for DGE analysis. I am fearful that the unstranded nature of my first lot of libraries will bias those datasets for genes which have known antisense transcripts. So what is the recommended course of action here? I'm not overly interested in antisense stuff so should i treat everything as unstranded for the featurecounts run?
Thank you in advance.
I am working with S. pombe doing 100bp single end Illumina HiSeq 2500 RNA-seq. Despite being told the trueseq stranded protocol was used one lot of my libraries appear unstranded (roughly 50/50 split between fwd and rev mapping reads in IGV). The second lot appears stranded as expected (>95% mapping bias to antisense strand in IGV). I'm trying to understand the best way to do read summarisation for these libraries using featurecounts for DGE analysis. I am fearful that the unstranded nature of my first lot of libraries will bias those datasets for genes which have known antisense transcripts. So what is the recommended course of action here? I'm not overly interested in antisense stuff so should i treat everything as unstranded for the featurecounts run?
Thank you in advance.
Comment