This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Get genes and coverage/depth per gene from sam/bam

Hi, Is there any tool (preferably one that could be run locally) that gets genes and coverage/depth per gene from sam/bam file (like sum of reads that map to given gene)? I mean coverage per gene, not position/locus, with no prior knowledge of what genes to look for (no bed file).

ngs sam bam coverage per gene depth per gene
• 2,262 views • link
updated by Charles Warden Charles Warden • written by prishly • prishly

And how is that tool supposed to delineate genes, then?

• 1 views • link
written by WouterDeCoster WouterDeCoster

There isn't, but this is a very good idea for a machine learning method.

• 1 views • link
written by h.mon h.mon

1 answer

I'm not sure what issues there could be with looking at coverage when the de novo annotations are from the same sample (compared to an independent set of annotations), but you could theoretically generate regions with derfinder (or do a reference-based de novo assembly, with a program like cufflinks).

Once you had your regions, you could use RSeQC to generate coverage plots.

• 1 views • link
written by Charles Warden Charles Warden

Log in to answer this question.

More posts like this