A Question About Samtools And Calculating Depth
1
0
Entering edit mode
11.6 years ago
Sameet ▴ 300

I have some RNA-Seq data. I want to know if there is a way to use samtools to calculate strand-specific depth from the accepted_hits.bam file generated by tophat.

samtools depth-of-coverage • 3.8k views
ADD COMMENT
0
Entering edit mode
11.6 years ago

Not sure what exactly you mean. Can you be little more descriptive. So far whatever I understand I think you should use Cufflinks for the next step. You can also make use of SAM flags (Check out - http://picard.sourceforge.net/explain-flags.html.) that tell you which strand, forward or reverse a read belongs to ? You can use this information to split your bam file into two - forward and reverse and calculate coverage.

ADD COMMENT
0
Entering edit mode

I do have the cufflinks output for these. But i only wanted to see if i can visualize the alignment itself.

ADD REPLY
0
Entering edit mode

Read this post. If you need to visualize it on genome browser, you can create strand specific wiggle files as described in the post.

ADD REPLY
0
Entering edit mode

Thank you @ashutoshmits. This was very helpful.

ADD REPLY

Login before adding your answer.

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