This is a test version of Biostars. For the public version, visit https://www.biostars.org.
A Question About Samtools And Calculating Depth

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
• 4,653 views • link
updated by Ashutosh Pandey Ashutosh Pandey • written by Sameet • Sameet

1 answer

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.

• 0 views • link
written by Ashutosh Pandey Ashutosh Pandey

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

• 0 views • link
written by Sameet • Sameet

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

• 0 views • link
updated by Ram Ram • written by Ashutosh Pandey Ashutosh Pandey

Thank you @ashutoshmits. This was very helpful.

• 0 views • link
written by Sameet • Sameet

Log in to answer this question.

More posts like this