Strandedness of RNAseq data from SAM alignment flag
1
0
Entering edit mode
3 days ago
rs.yousaf • 0

Hi,

I am very new to RNAseq analysis. I have paired end RNAseq data from illumina Hiseq that I aligned using BWA MEM to a virus genome. This generated a BAM file that has alignment flags. I want to distinguish positive strand reads from negative strand as we have RNA from a virus that makes both strands. I have been using the following tags to filter negative and positive reads:

83 and 163 for positive reads 99 and 147 for negative reads

Can someone please confirm if this is correct? Also, how is strandedness determined during library prep. I really struggle to understand this.

Thank you!

SAMflag Strandanalysis RNAseq • 92 views
ADD COMMENT
0
Entering edit mode
2 days ago

Also, how is strandedness determined during library prep.

It depends on the kit. You really need to stop and understand your experiment before you analyze it.

ADD COMMENT

Login before adding your answer.

Traffic: 2169 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6