Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • blancha
    Senior Member
    • May 2013
    • 367

    #16
    @GenoMax
    IGV is just about the only software program that I know of smart enough not to care about the presence or absence of "chr" in the name of the chromosome, at least for the human and mice genomes. I would congratulate the IGV programmer who programmed IGV not to care whether the "chr" is there or not.

    Apparently, the human race is smart enough to sequence the human genome, not smart enough to decide whether chromosome 1 should be written "chr1" or "1".

    Comment

    • Milestailsprowe
      Member
      • Jun 2015
      • 31

      #17
      Originally posted by blancha View Post
      In IGV, just pick the mm10 genome instead.
      It's the same genome as GRCm38.
      I think you'll have to reload the BAM file after you've selected the correct genome.

      Then, zoom into a location where you know you will have coverage.
      For RNA-Seq for example, you could pick a housekeeping gene, like GAPDH.
      Just type GAPDH in the search box, and click on the Go button.
      Didnt work

      Originally posted by GenoMax View Post
      @blancha: mm10 may not work if the one included in IGV is UCSC version which has the "chr" in front of all chromosome numbers.

      @Milestailsprowe: If above does not work, create a new "genome" by pointing to the iGenomes (/path_to/WholeGenomesFasta/genome.fa file) and use the corresponding GTF file from (/path_to/Annotations/Genes/genes.gtf). Open your BAM file in IGV.


      Where do I put the .gtf file?

      Thank you for all your help. I'm a graduate student trying to figure this out and my professor never has time sadly

      Comment

      • blancha
        Senior Member
        • May 2013
        • 367

        #18
        @GenoMax
        Your comment got me thinking, and I've finally found out how come IGV is the only program that can handle smoothly both UCSC and Ensembl chromosome nomenclature. IGV uses a chromosome name alias file.
        Simple, but so useful.

        "Note: Certain well-known aliases are built into IGV and do not require an alias file. These include mappings that involve adding or removing the prefix "chr" to the name, for example 1 -> chr1 and chr1 -> 1. Also, NCBI identifiers that start with "gi|" and follow the pattern illustrated in the example above are automatically mapped. "
        https://www.broadinstitute.org/igv/LoadData

        @Milestailsprowe
        It really shouldn't be necessary to load your own genome fasta file, or your own GTF file. If you absolutely insist on doing it, the Gene file would be the GTF file.

        I'm wondering if you aren't using an ancient version of IGV. The chromosome alias file only appeared in IGV 1.53, which would explain why you can't view your data, if the chromosome nomenclature is different. @Genomax's solution to load your own genome file would then be correct, although a better solution would be to upgrade IGV.

        Check which version of IGV you are using. If it's a prehistoric version, just upgrade to 2.3 before doing any more troubleshooting. Then select mm10, reload the BAM file, and view a region that you know has high coverage.

        I don't recognize the interface you show in your screenshot in recent versions of IGV, and I believe the interface should be the same across all operating systems. I've only seen this interface in very old versions of IGV.

        I could be completely wrong, of course. It wouldn't be the first time.
        Last edited by blancha; 11-02-2015, 09:55 PM.

        Comment

        • GenoMax
          Senior Member
          • Feb 2008
          • 7142

          #19
          Originally posted by Milestailsprowe View Post
          Where do I put the .gtf file?

          Thank you for all your help. I'm a graduate student trying to figure this out and my professor never has time sadly
          Point to the GTF file in the option box that says "Gene" file.

          You can also read in the GTF file along with your BAM later on.

          Comment

          • Milestailsprowe
            Member
            • Jun 2015
            • 31

            #20
            Originally posted by blancha View Post
            @GenoMax
            Your comment got me thinking, and I've finally found out how come IGV is the only program that can handle smoothly both UCSC and Ensembl chromosome nomenclature. IGV uses a chromosome name alias file.
            Simple, but so useful.

            "Note: Certain well-known aliases are built into IGV and do not require an alias file. These include mappings that involve adding or removing the prefix "chr" to the name, for example 1 -> chr1 and chr1 -> 1. Also, NCBI identifiers that start with "gi|" and follow the pattern illustrated in the example above are automatically mapped. "
            https://www.broadinstitute.org/igv/LoadData

            @Milestailsprowe
            It really shouldn't be necessary to load your own genome fasta file, or your own GTF file. If you absolutely insist on doing it, the Gene file would be the GTF file.

            I'm wondering if you aren't using an ancient version of IGV. The chromosome alias file only appeared in IGV 1.53, which would explain why you can't view your data, if the chromosome nomenclature is different. @Genomax's solution to load your own genome file would then be correct, although a better solution would be to upgrade IGV.

            Check which version of IGV you are using. If it's a prehistoric version, just upgrade to 2.3 before doing any more troubleshooting. Then select mm10, reload the BAM file, and view a region that you know has high coverage.

            I don't recognize the interface you show in your screenshot in recent versions of IGV, and I believe the interface should be the same across all operating systems. I've only seen this interface in very old versions of IGV.

            I could be completely wrong, of course. It wouldn't be the first time.
            I'm usuing version 2.3.14 (16)

            Can you please give me a run down on how you would do it from start the finish. NOthing is showing up

            Where would you download your orginal file? Please. dont understand how I can have file multiple gigs big and no data show up

            Comment

            • GenoMax
              Senior Member
              • Feb 2008
              • 7142

              #21
              Can you run the samtools flagstat/idxstats commands on your BAM file and post the output? That would give us an idea of what the status is as far as alignments are concerned.

              Did you try my suggestion with the iGenomes files/IGV?

              Comment

              • Milestailsprowe
                Member
                • Jun 2015
                • 31

                #22
                Originally posted by GenoMax View Post
                Can you run the samtools flagstat/idxstats commands on your BAM file and post the output? That would give us an idea of what the status is as far as alignments are concerned.

                Did you try my suggestion with the iGenomes files/IGV?
                What does it mean

                Comment

                • Milestailsprowe
                  Member
                  • Jun 2015
                  • 31

                  #23
                  Oh. I have one file that I got to work in goldenhelix. I just compare it to the the one I showed you and there is a big difference in the %. What does that mean

                  Comment

                  • GenoMax
                    Senior Member
                    • Feb 2008
                    • 7142

                    #24
                    Originally posted by Milestailsprowe View Post
                    What does it mean
                    https://www.biostars.org/p/12475/

                    Comment

                    • GenoMax
                      Senior Member
                      • Feb 2008
                      • 7142

                      #25
                      Originally posted by GenoMax View Post
                      Looks like one of your files has good # of reads mapped where as the other one has far fewer. You should be able to look at both of them in IGV at the same time (if you are finally able to figure out how to use IGV).

                      Comment

                      Latest Articles

                      Collapse

                      • SEQadmin2
                        Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
                        by SEQadmin2



                        Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

                        Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
                        ...
                        10-14-2026, 09:55 AM
                      • SEQadmin2
                        New Genomics Technologies Take Aim at Long-Standing Limits
                        by SEQadmin2


                        Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

                        We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
                        ...
                        09-28-2026, 10:25 AM

                      ad_right_rmr

                      Collapse

                      News

                      Collapse

                      Topics Statistics Last Post
                      Started by SEQadmin2, Yesterday, 11:13 AM
                      0 responses
                      16 views
                      0 reactions
                      Last Post SEQadmin2  
                      Started by SEQadmin2, Yesterday, 10:48 AM
                      0 responses
                      17 views
                      0 reactions
                      Last Post SEQadmin2  
                      Started by SEQadmin2, 10-15-2026, 10:52 AM
                      0 responses
                      35 views
                      0 reactions
                      Last Post SEQadmin2  
                      Started by SEQadmin2, 10-14-2026, 05:39 AM
                      0 responses
                      42 views
                      0 reactions
                      Last Post SEQadmin2  
                      Working...