Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • coverage plot

    Can anyone recommend a good tool to display a coverage plot over a chromosome? The images produced by UCSC aren't of great quality nor can I customize it. I'm curious what most people use to generate high quality coverage plots with.

  • #2
    What's your desired input format for the script?

    Comment


    • #3
      to start with, the length of the chromosome and the coverage information, perhaps in wig format.

      For fancier stuff, say customizing color, size of the image (in pixels), how the +/- strand information is displayed, etc.

      Comment


      • #4
        Maybe the alignment view of IGV is enough for want you want, it has a coverage track:
        If you are unable to find something or have a question about our new website, please email [email protected]. For other inquiries related to the Broad Institute, the necessary contact information can be found here.


        For a completely different idea you could also have a look at (the R package for) Hilbert curves:

        Comment


        • #5
          I've done coverage plots of small (circular) genomes using Biopython and GenomeDiagram, things like mitochondria or viruses - should work fine for bacteria too. However, this isn't an off the shelf point and click solution.

          Comment


          • #6
            One of the features for the next release of SeqMonk will be a whole genome view. This will display the current quantitation for a single sample over the whole genome so you can configure it to show pretty much anything you want.

            Comment


            • #7
              Savant can do that as well. If've got a BAM file of your aligments you can convert it into a BAM Coverage files. That allows you to display coverage along a chromosome easily. I attached a picture of a whole human chromosome 1 coverage plot.

              Hope that helps,
              Best regards,
              Peter
              Attached Files

              Comment


              • #8
                Originally posted by svl View Post
                Maybe the alignment view of IGV is enough for want you want, it has a coverage track:
                If you are unable to find something or have a question about our new website, please email [email protected]. For other inquiries related to the Broad Institute, the necessary contact information can be found here.


                For a completely different idea you could also have a look at (the R package for) Hilbert curves:
                http://www.ebi.ac.uk/huber-srv/hilbert/
                I think IGV has the issue of displaying alignment only if you zoom-in enough.. not sure if that is true for the coverage track as well!
                --
                bioinfosm

                Comment


                • #9
                  IGV does a great job of displaying coverage information, but you do need to zoom in enough. I want something more on the chromosome level to get an overview of my coverage.

                  Thank you for all the suggestions. I'll have to spend some time looking at these.

                  Comment


                  • #10
                    hi, simonandrews

                    Do you have any experience in trying SeqMonk in linux? I use SSH to connect our linux server. I just downloaded SeqMonk to my directory and tried to launch it by using

                    java -Xmx1500m -classpath $CLASSPATH:. uk.ac.bbsrc.babraham.SeqMonk.SeqMonkApplication
                    but it didn't pop up any window.. I was wondering what the correct way is to use SeqMonk in Linux.

                    Thanks

                    -c

                    Comment


                    • #11
                      cliff - Make sure you use -X for X11 forwarding with ssh

                      Comment


                      • #12
                        Originally posted by golharam View Post
                        cliff - Make sure you use -X for X11 forwarding with ssh
                        hi, golharam

                        How to use -X for X11 forwading with SSH? I am not familiar with X11. I turned on "Tunnel X11 Connections" in ssh settings but it didn't work..

                        Comment


                        • #13
                          Are you using Windows? You'll need to install an X Server like ReflectionX, eXceed, or Cygwin. Googling for "X Forwarding on Windows" turned up http://www.cs.caltech.edu/courses/cs.../xwindows.html

                          Comment


                          • #14
                            Originally posted by cliff View Post
                            hi, simonandrews

                            Do you have any experience in trying SeqMonk in linux?
                            If you're running Linux locally on your machine then what you tried should have worked (and we use that here very successfully).

                            As others have pointed out, if you want to run SeqMonk on linux from a remote system over SSH you'll need to make sure that your ssh session is started with an X tunnel available (either ssh -X or set the appropriate flags in PuttY) and that you have a local X server running on your machine. This is the same for all graphical linux programs.

                            Having said that, I'm not sure I'd encourage people to run SeqMonk this way. The nature of the program means that it does big complex updates to its display all the time, which means that it has to send an awful lot of display data over X to work remotely. We know of one site who do it this way, but they run an nx server to improve the display performance.

                            Really, SeqMonk is better run as a local native application. It's designed to work on modest hardware so you don't need a particularly powerful PC to run it (the default configuration is optimised for a dual core machine with 2GB of RAM). If you've not tried running it directly on your desktop I'd try that before going to the hassle of getting a decent remote connection working.

                            Comment


                            • #15
                              Hi everyone,

                              IGV can also visualize the coverage of a whole chromosome at once without having to zoom in:
                              You need to create an extra .tdf-file (a binary tiled data file) using igvtools with the command "count" (see http://www.broadinstitute.org/igv/igvtools).

                              For example
                              igvtools count -z 7 -w 25 -e 0 alignments.bam alignments.coverage.tdf hg19

                              It supports the formats .sam, .bam, .aligned, .sorted.txt, and .bed.

                              Then you can load the .tdf-file into IGV and you get a pretty nice coverage plot.
                              Last edited by ForeignMan; 09-30-2010, 12:00 AM.

                              Comment

                              Latest Articles

                              Collapse

                              • seqadmin
                                Exploring the Dynamics of the Tumor Microenvironment
                                by seqadmin




                                The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
                                07-08-2024, 03:19 PM
                              • seqadmin
                                Exploring Human Diversity Through Large-Scale Omics
                                by seqadmin


                                In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
                                06-25-2024, 06:43 AM

                              ad_right_rmr

                              Collapse

                              News

                              Collapse

                              Topics Statistics Last Post
                              Started by seqadmin, 07-10-2024, 07:30 AM
                              0 responses
                              24 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 07-03-2024, 09:45 AM
                              0 responses
                              201 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 07-03-2024, 08:54 AM
                              0 responses
                              211 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 07-02-2024, 03:00 PM
                              0 responses
                              192 views
                              0 likes
                              Last Post seqadmin  
                              Working...
                              X