Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • loading bam file failure on IGV

    Hi everyone,

    I tried to find the reason for a long time. But couldn't solve my problem.
    When loading my sorted bam file (with bai file in the same folder) on IGV, I could see nothing.
    My bam file data looks like this:
    HWI-ST261:6:66:13462:17077#TAGCTT 83 chr01 148 255 95M = 94 -149 TGTAGTTCAAAAAATTAAAGTATGTATTCAATAGAGCTGAGGAAATTGTTCATGGAGGAATGATAGTGCTTATTACACCATTTTCAAGTTATGTA EE?D@EEBDE?DHHCHFHFFDB>HEE>EEEDDHDBHDDDHEFEFHEEHEHFHHHHHHHHHHHHHFEEEEFFHFFDDEHHHFHHHHHHHHHHFHHH XA:i:0 MD:Z:8T86 NM:i:0

    I tried with a bam file downloaded online. It worked when I loaded the file on IGV.
    IL7_1788:1:61:1161:1010 99 MAL1 17 0 6H48M = 548 584 ACCCTGAACCCTAAACCCTAAACCCTGAACCCTAAACCCTGAACCCTG>>>>>>>>>>>>>>>>>>>>>>>><<<<<666626666666+66/662 AS:i:48 MS:i:11

    My DNA-seq data is very big. My bam file was produced by galaxy. I used the steps as follows in galaxy:
    map with bowtie
    filter sam
    sam to bam

    Any suggestions will be very appreciated!

  • #2
    Of course it is always hard to troubleshoot a problem, but the first thing that stands out to me as a difference between your good data and bad data are the reference names. In the good one you have 'MAL1' while in the bad one the name is 'Chr01'. If the names are not matching your reference then you will not see anything.

    Are you using the same reference in both cases? Or different ones? If so, what are the names inside the reference file itself?
    Last edited by westerman; 12-22-2011, 08:00 AM. Reason: Fixed some English.

    Comment


    • #3
      Like westerman said, I think the problem might be the reference names. If this human data and you are using one of the pre-loaded human genomes on the IGV, then they have "chr1" and not "chr01" as the name of the chromosomes. So try changing the names of the reference in your file and see if that works.

      Comment


      • #4
        I will echo the suggestion that this might be a chromosome name problem. What assembly were these aligned against?

        You can correct this by creating an "alias" file, as decribed here (at the bottom of the page)



        Jim

        Comment


        • #5
          I followed your suggestions. You're correct. It's a chromosome name problem. I'm so happy that I solved this problem finally. Thanks a lot.

          Comment


          • #6
            Originally posted by ykingh View Post
            I followed your suggestions. You're correct. It's a chromosome name problem. I'm so happy that I solved this problem finally. Thanks a lot.
            By the way, this is, by far, the most common failure mode for IGV that I see. BAM file has one set of names for the chromosomes and the reference gff or gtf has another set. The result? No alignments are depicted. The solution? There are several, but all of them involve delving into the data files and either changing them or creating a table file and putting it in the right place for IGV to read it.

            At the very least I would propose that IGV should display a visible warning when a .bam file is loaded that contains no sequences mapped to any of the current reference sequences. Some along the lines of "the reference contains the sequences "chr1, chr2 ..." while the alignment file contains hits to the sequences "1, 2, ..." -- since none of these sequence names are the same, no alignments can be displayed."

            --
            Phillip

            Comment


            • #7
              I agree with the warning suggestion, I will do that. Its on our faq, and its not just an IGV issue (samtools, the GATK, UCSC, etc won't work either), but nevertheless a warning would be helpful. Fundamentally it means you are viewing your BAMs against a different fasta file than was used for the alignment step, so some caution is in order.

              -- Jim

              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
              49 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