Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • qinbird
    Junior Member
    • Apr 2010
    • 1

    Tophat: how to define primary alignment

    From the sam file generate by Tophat, if there are multiple mapped reads, one of the alignments will be labeled as primary alignment.
    does anybody know how Tophat define primary alignment?
    is there any reference mentioned it?
    thanks!

    i checked one read, there are four alignments.
    all of them have 2 mismatch, but only one are labeled as primary alignment.
  • fongchun
    Member
    • May 2011
    • 55

    #2
    Did you get ever figure this out? I would be interested in knowing this answer.

    Thanks,

    Comment

    • dvanic
      Member
      • Jan 2012
      • 61

      #3
      I am also quite interested in this! Still...

      Comment

      • dvanic
        Member
        • Jan 2012
        • 61

        #4
        Because some of my "non-primary" alignments look quite good in IGV...

        Comment

        • apredeus
          Senior Member
          • Jul 2012
          • 151

          #5
          There seems to be confusion in what's meant by "primary" in the field.

          E.g. SAM flag interpreter http://picard.sourceforge.net/explain-flags.html clearly uses "non-primary" as an equivalent for spliced

          Comment

          • Brian Bushnell
            Super Moderator
            • Jan 2014
            • 2709

            #6
            Originally posted by apredeus View Post
            There seems to be confusion in what's meant by "primary" in the field.

            E.g. SAM flag interpreter http://picard.sourceforge.net/explain-flags.html clearly uses "non-primary" as an equivalent for spliced
            I don't see the word 'spliced' anywhere on that page. The difference between primary and non-primary is that only one primary alignment is allowed. It was pretty clear, at least until 'supplemental alignments' and made things murky. The intent of supplemental appears to be for representing chimeric alignments, and thus it could be used for splice events, while secondary shouldn't be used that way.

            Comment

            • apredeus
              Senior Member
              • Jul 2012
              • 151

              #7
              I don't see the word 'spliced' anywhere on that page.
              Flags of 256/272 in splice aligners' output usually means gapped (spliced) alignment to + and - strand respectively. That's what they call "not primary alignment" on that page.

              The difference between primary and non-primary is that only one primary alignment is allowed.
              Not sure what that means? I figured that if there's one best alignment score the read would be considered as mapping to one genomic position, while if there are few scores that are equally high would be a multimapper. Now, most RNA-seq aligners do allow multimappers up to a certain number of locations, e.g. Tophat's default is 15.

              If you call the read that has the best score primary, it's not clear what's the use of having a term like that since others would be discarded anyway?

              Comment

              • Brian Bushnell
                Super Moderator
                • Jan 2014
                • 2709

                #8
                Originally posted by apredeus View Post
                If you call the read that has the best score primary, it's not clear what's the use of having a term like that since others would be discarded anyway?
                Depends on the application; sometimes they are discarded, sometimes not - they can be used fractionally when calculating coverage, for example. The most important case is when there is no single best-scoring location. Then one location can be assigned 'primary' and the rest 'secondary'. This allows a tool (or person) aware of secondary alignments to use that information later, or for simplicity (for example, when converting sam to fastq), all the secondary sites can just be ignored.

                Comment

                • apredeus
                  Senior Member
                  • Jul 2012
                  • 151

                  #9
                  You're right, I was actually totally wrong to assume that flag to mean spliced as opposed to non-spliced alignment. It's "secondary" as defined by some other criteria.

                  Thank you

                  Comment

                  • dpryan
                    Devon Ryan
                    • Jul 2011
                    • 3478

                    #10
                    Originally posted by apredeus View Post
                    Flags of 256/272 in splice aligners' output usually means gapped (spliced) alignment to + and - strand respectively. That's what they call "not primary alignment" on that page.
                    That's not splicing, it's a chimeric/fusion/non-linear alignment (there are a bunch of terms for this). These days, those should have the 0x800 bit set in the flag.

                    Comment

                    Latest Articles

                    Collapse

                    • seqadmin
                      Pathogen Surveillance with Advanced Genomic Tools
                      by seqadmin




                      The COVID-19 pandemic highlighted the need for proactive pathogen surveillance systems. As ongoing threats like avian influenza and newly emerging infections continue to pose risks, researchers are working to improve how quickly and accurately pathogens can be identified and tracked. In a recent SEQanswers webinar, two experts discussed how next-generation sequencing (NGS) and machine learning are shaping efforts to monitor viral variation and trace the origins of infectious...
                      03-24-2025, 11:48 AM
                    • seqadmin
                      New Genomics Tools and Methods Shared at AGBT 2025
                      by seqadmin


                      This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

                      The Headliner
                      The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
                      03-03-2025, 01:39 PM

                    ad_right_rmr

                    Collapse

                    News

                    Collapse

                    Topics Statistics Last Post
                    Started by seqadmin, 03-20-2025, 05:03 AM
                    0 responses
                    49 views
                    0 reactions
                    Last Post seqadmin  
                    Started by seqadmin, 03-19-2025, 07:27 AM
                    0 responses
                    57 views
                    0 reactions
                    Last Post seqadmin  
                    Started by seqadmin, 03-18-2025, 12:50 PM
                    0 responses
                    50 views
                    0 reactions
                    Last Post seqadmin  
                    Started by seqadmin, 03-03-2025, 01:15 PM
                    0 responses
                    200 views
                    0 reactions
                    Last Post seqadmin  
                    Working...