Display read depth
Hi: I would like to overlay RNA-seq data in the genome browser. Is this possible? If not, is there a mechanism to display read depth in the genome browser. Something like a track in the genome browser that would show coverage vs genome position? If so, would this be visualized in both the desktop and web version of ptools? Thanks.
Bedtools has a pretty slick utility that converts a BAM files to genome read coverage as described in http://bedtools.readthedocs.io/en/latest/content/tools/genomecov.html?highlight=genomecov (genomecov)
I could convert a BAM file of RNA-seq data to a http://genome.ucsc.edu/goldenPath/help/bedgraph.html (bedgraph) format
chr1 554304 554309 25
(chromA chromStartA chromEndA dataValueA, where dataValueA is read coverage)
and then stuff text into each line to convert to something GFF2 like
chr1 rna-seq cover 554304 554309 25
which presumably could be ingested into ptools and displayed as a Track. That got me thinking that the extra columns in the GFF2, although necessary for GFF2, really are not informative in this use case. That my way of trying to convince SRI to tweak ptools to accept bedgraph format files.