Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • #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


    • #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


      • #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. "
        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.


        @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


        • #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


          • #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. "
            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.


            @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


            • #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


              • #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


                • #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


                  • #24
                    Originally posted by Milestailsprowe View Post
                    What does it mean

                    Comment


                    • #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

                      • seqadmin
                        Current Approaches to Protein Sequencing
                        by seqadmin


                        Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
                        04-04-2024, 04:25 PM
                      • 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

                      ad_right_rmr

                      Collapse

                      News

                      Collapse

                      Topics Statistics Last Post
                      Started by seqadmin, 04-11-2024, 12:08 PM
                      0 responses
                      25 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 04-10-2024, 10:19 PM
                      0 responses
                      28 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 04-10-2024, 09:21 AM
                      0 responses
                      24 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 04-04-2024, 09:00 AM
                      0 responses
                      52 views
                      0 likes
                      Last Post seqadmin  
                      Working...
                      X