Dear Ann and Jeremy:
We have this discussion long time ago, and I am sorry that I brought it
up here
again. I am just thinking that as Ann said, can we add this tool which convert
bam into wig file into galaxy? Or make a workflow to generate a wig
file from a
bam file generate by tophat? In
Hi Ying,
You're in luck because I've been working with genome browsers lately, so I
think I can help you address your problem. What you're looking for is a
visualization of a coverage histogram for the BAM reads produced by Tophat,
yes?
It turns out that some genome browsers provide this
I can answer IGV questions, sadly I'm still coming up to speed on
Galaxy.
I've lost track of the original question, but IGV computes a coverage
histogram on the fly, a bit like the Galaxy Track Browser, but you
have to be zoomed in. However, you can also precompute a coverage
Ying and all;
Nice discussion and lots of great ideas. One other approach is to generate
bigWig files. You can serve these directly from Galaxy and they have
the same advantages as BAM files; they are compressed and individual
regions that be fetched on demand by UCSC.
This tool in the community
HI,
The option you need in IGV tools is count. You set a window size and this
gives you a tdf file from your sorted bam (or sam) file which is nice and quick
to view on IGV.
Best Wishes,
David.
__
Dr David A. Matthews
Senior Lecturer in Virology
Room E49
Ying,
You can create a coverage file from a BAM file using this workflow:
http://main.g2.bx.psu.edu/u/jeremy/w/create-coverage-dataset-from-bam-dataset
While this doesn't generate a proper wiggle file, it can be a useful file for
statistical analysis.
I and others have suggested ways to
Hi,
You can get an equivalent visualisation from the IGV viewer by the Broad
Institute - its under IGV tools and generates a tdf file from bam or sam files.
This also gives a quick and easy way of looking at depth at any particular site
and is very accessible.
Cheers
David
On 21 Feb 2011,
-user] get wig file after tophat
Hi:
I am using tophat in galaxy to analyze my paired-end RNA-seq data and find out
that after the tophat analysis, we can not get the wig file from it anymore
which is used to be able to. Do you have any idea of how to still be able to
get the wig file after tophat
8 matches
Mail list logo