Annotate BAM with gene strand quickly
0
0
Entering edit mode
3 months ago
mb • 0

I would like to annotate RNA-seq alignments with the strand on the genome they come from. I can do a brute force method using GTF file and for every read look up the strand information and store it as a tag. However I would like to do in a fast and efficient way.

RNA-seq • 275 views
ADD COMMENT
0
Entering edit mode

The bitwise flag in the BAM file (2nd field) already has strand information. I guess I don't understand why you would need to look it up.

ADD REPLY
0
Entering edit mode

No the flag only tells you either forward or reverse complement of the 'read' is mapped to the genome. If I knew the directionality of the reads in the first place I could use that information. However my reads are split half half between strands of the input cDNA so regardless of gene strand half maps to positive strand and half to negative. However every gene has a true strand which I need to know.

ADD REPLY

Login before adding your answer.

Traffic: 1799 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