Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Visualization tool to compare distribution profiles of NGS reads?

    Hello,

    Does any one know of a (free) software that is able to normalize aligned NGS data of different samples (e.g. in .bam format) to enable comparison of their read distribution profile across a reference sequence? That is, I want to compare not only the distribution profile of reads mapped onto the reference but also compare their relative abundance to one another, simultaneously.

    I've tried Integrated genome viewer, but would like to try other software as well.

    Thanks for any information!

    Ken

  • #2
    i try artemis and tablet but i failed to view in IGV even though I already sorted and index the BAM file.
    any tips can you share with me to view in IGV because so many times i tried but still can't view from IGV.

    Comment


    • #3
      Tablet only loads one BAM file at a time, so to do comparisons with Tablet you'd need to have one Tablet window open per BAM file, and manually line up the views. Not ideal.

      I would use Artemis with the BamView plugin. This can load a reference sequence with multiple BAM files, which you can color code. You can show a coverage plot as well as the raw reads. The user interface is a bit cumbersome, but it works.

      Comment


      • #4
        If your reference sequence is an assembled genome then you can do this in SeqMonk (which comes from our group)

        Comment


        • #5
          Seq monk

          Hi Simon,

          Can we plot RPKM/FPKM values for visualization in Seqmonk?

          Comment


          • #6
            Originally posted by honey View Post
            Can we plot RPKM/FPKM values for visualization in Seqmonk?
            Kind of. In the current release you can calculate values which are adjusted to the size of the largest dataset rather than per million reads. These are equivalent to RPKM but scaled up (or down) by a constant factor. In the next release there will be an option to correct per million reads to get RPKM directly.

            The reason we don't do RPKM at the moment is that if you log transform your data (which makes a lot of sense when viewing most types of expression data), then you often end up with a lot of negative values which make the plots more difficult to interpret than if everything is on a positive scale.

            Comment


            • #7
              Hi rururara

              I've been viewing the .tdf file with the 'count' function from igvtools.
              did you try to generate the .tdf coverage file with igvtools, or just tried to view along with the bam file? Sometimes if the reference is too large it won't display it until you zoom in quite abit.
              Let me know what you did specifically and i'll get back to you with more details.

              thanks everyone else

              Comment


              • #8
                If you create a coverage (tdf) file with igvtools you can normalize it when viewing by changing a preference. Specifically, select View > Preferences, click the "Tracks" tab, and select the last checkbox "Normalize Coverage Data". This will scale each value by 1,000,000 / total mapped count.

                I realize this is not a very intutitve way to turn this feature on, its on my list to fix.

                Also on my list is implementing RPKM and FPKM, both in igvtools and "on-the-fly" in igv. I understand the concept but was here on seqanswers hoping to to find a definition that acounts for alternative isoforms. Is the normal thing to take the union of all exons for all transcripts, or what?

                Jim

                Comment


                • #9
                  Jim,

                  Not sure but these might be relevant to your question:

                  Accurate quantification of transcriptome from RNA-Seq data by effective length normalization


                  Prediction of alternative isoforms from exon expression levels in RNA-Seq experiment (they used IGV in this one)

                  Comment


                  • #10
                    Igv

                    Hi Jim,
                    I was wodnering have you implemented changes in IGV so as to plot FPKM/RPKM normailzed values?
                    Thanks.

                    Comment

                    Latest Articles

                    Collapse

                    • seqadmin
                      Strategies for Sequencing Challenging Samples
                      by seqadmin


                      Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                      03-22-2024, 06:39 AM
                    • seqadmin
                      Techniques and Challenges in Conservation Genomics
                      by seqadmin



                      The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                      Avian Conservation
                      Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                      03-08-2024, 10:41 AM

                    ad_right_rmr

                    Collapse

                    News

                    Collapse

                    Topics Statistics Last Post
                    Started by seqadmin, Yesterday, 06:37 PM
                    0 responses
                    10 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, Yesterday, 06:07 PM
                    0 responses
                    9 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 03-22-2024, 10:03 AM
                    0 responses
                    50 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 03-21-2024, 07:32 AM
                    0 responses
                    67 views
                    0 likes
                    Last Post seqadmin  
                    Working...
                    X