Entering edit mode
3.5 years ago
venura
▴
70
Hi,
I have this question (prob a silly question). We have a set of RNASEq data that came from BGI sequencing platform. The dataset is PE 150bp (and stranded).
I am going to use HISAT2 for alignment and wondering about the strandness setting to be used. This is the exp procedure they used:
Any thoughts on what setting I should use R/RF or F/FR . By the Image I feel the latter should be used. Just trying to make sure. Thanks in advance.
Run GUESSmyLT and you should have your answer.
Thank you for the link, Jacques!
Is this even a stranded library? I would ask the facility that made the libraries for advise.
I did. They said to follow the Illumina pipeline. :| When we placed the order we asked for stranded sequencing and PE 150
Following Illumina pipeline is one thing but as ATPoint said earlier this library does not appear to be stranded (see: HISAT2 rna-strandness option )
That is not an answer I would accept if I was you. You paid for it, so I personally would insist on details which means a protocol and details on the library prep. As of the picture you provide it does not seem that the library is stranded (at least based on the image) because (for me) it is not obvious whether the the first strand synthesis preserves the strand information. Ask for details, don't let them waive you off until you get the details you need for your analysis, you are the paying customer ;-)
Writing back to them now! You are correct, I am going to give them a hard time! :)
They replied and said it was a mistake :) and sent their original logs as proof. and asked me to replace the old image with this.
Please edit the original post and replace the image in it.
Looks like this is a standard dUTP library.
Just updated the original post.
So standard dUTP I guess: https://bioinformatics.stackexchange.com/questions/4074/hisat2-which-option-should-mention-for-strand-specific-library-read